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Preface to the second edition 


The objectives outlined in the preface to the first edition have remained 
unchanged in preparing this edition as they have continued to be the basis of 
my teaching programme. This edition is therefore not radically different 
from the first, which to my pleasure and relief was well received by those 
who obtained and used the book. 

I have taken the opportunity to correct errors that occurred, have 
improved some diagrams and added others, and brought all the material on 
cavitation together into Chapter 3: I hope that this gives a more connected 
account of this very important topic. I have added some updated material in 
places, have added some references, and hope that by this means the reader 
can pursue some topics in more depth after reading this introduction. 

The worked examples that were included in the text have been retained, 
and extra exercises have been added where students have commented on the 
need for further clarification. A major change has been the addition of sets 
of problems for solution by the reader. These are given at the end of all 
chapters but four, five and ten. These are based in most cases on the 
questions set over the years in the Finals in the course on Turbomachinery 
at Loughborough University of Technology, and I am grateful for the 
permission granted by the University authorities to use them. While the 
problems are placed at the end of each chapter, the solutions are collected 
together at the end of the book. It is hoped that readers will attempt the 
problems first and then turn to the end for help. 

I hope that this edition is free from error and ambiguity, and as an earnest 
seeker after truth will be grateful for comments and suggestions. 

I must acknowledge the invaluable help of Mrs Janet Redman for her 
translation of my sketches and of Mrs Gail Kirton who typed the new 
chapters. Finally, my thanks to my dear wife who has been patient and 
helpful as always. 



Preface to the first edition 


This text outlines the fluid and thermodynamic principles that apply to all 
classes of turbomachines, and the material has been presented in a unified 
way. The approach has been used with successive groups of final year 
mechanical engineering students, who have helped with the development of 
the ideas outlined. As with these students, the reader is assumed to have a 
basic understanding of fluid mechanics and thermodynamics. However, the 
early chapters combine the relevant material with some new concepts, and 
provide basic reading references. 

Two related objectives have defined the scope of the treatment. The first 
is to provide a general treatment of the common forms of turbomachine, 
covering basic fluid dynamics and thermodynamics of flow through passages 
and over surfaces, with a brief derivation of the fundamental governing 
equations. The second objective is to apply this material to the various 
machines in enough detail to allow the major design and performance 
factors to be appreciated. Both objectives have been met by grouping the 
machines by flow path rather than by application, thus allowing an appreci- 
ation of points of similarity or difference in approach. No attempt has been 
made to cover detailed points of design or stressing, though the cited 
references and the body of information from which they have been taken 
give this sort of information. 

The first four chapters introduce the fundamental relations, and the suc- 
ceeding chapters deal with applications to the various flow paths. The last 
chapter covers the effects of cavitation, solids suspensions, gas content and 
pumped storage systems, and includes a short discussion of the control of 
output. These topics have been included to highlight the difficulties encoun- 
tered when the machine is not dealing with a clean Newtonian fluid, or in 
systems where problems are posed that can only be solved by compromise. 
Chapter 5 discusses all the conventional centrifugal machines, covering in a 
uniform manner the problems faced w'ith liquids and gases: since high 
pressure rise machines have a number of stages, the ways in which fluid 
is guided from stage to stage are introduced. Thrust load problems are 
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described and the common solutions adopted are outlined. The discussion of 
axial machines has been divided between two chapters, as the technologies 
of pumps, fans and water turbines are similar but differ from those used in 
compressible machines. Radial flow turbines form the subject matter of 
Chapter 8, and the common designs in use in industry and in turbochargers 
are discussed. 

Worked examples have been included in all chapters but the last. They 
are intended to provide illustration of the main points of the text, and to 
give a feel for both the shape of the velocity triangles and the sizes of the 
velocity vectors that normally apply. They are of necessity simplified, and 
must not be regarded as representing current practice in all respects. No 
problems for student solution have been provided. Teachers normally prefer 
to devise their own material, and may obtain copies of examination questions 
set by other institutions if they wish. 

As a matter of course the SI system of units has been used throughout, 
except in some diagrams. To assist the reader, a list of symbols used in the 
early chapters, together with a statement of the conventional dimensions 
used, follows the Preface. As far as possible the British Standard on symbols 
has been followed but, where current and hallowed practice dictates the use 
of certain symbols, these have been used; it is hoped that where the same 
symbol appears to have different meanings the context makes the usage 
clear. 

The material presented forms the core of a lecture course of about 46 
hours, and the author hopes that in the inevitable distillation no ambiguities 
have occurred. He will be grateful for comments and suggestions, as he is 
still an earnest ‘seeker after truth’. 

Finally, it is necessary to offer some words of thanks, especially to Mrs 
Redman, who ensured that the diagrams were clear, to Mrs Smith and Mrs 
McKnight, who helped with the typing, and finally to my dear wife, who was 
so patient and who ensured that the proof-reading was done properly. 




Symbols 

used: 


their meaning and 


dimensions 

a 

acoustic velocity 

m s 1 

b 

passage height 

m 

Cv. 

lift coefficient (Table 4.1) 


C D 

drag coefficient (Table 4.1) 


C P 

pressure rise coefficient (equation 4.15) 


Cp 

specific heat at constant pressure 

kJ kg~ 1 KT 1 

Cv 

specific heat at constant volume 

kJkg“‘ K ' 

D 

diameter 

m 

D 

drag force on an aerofoil 

N 

F a 

force acting in the axial direction on a foil 
or blade 

N 

F , 

force acting in the tangential direction on 
a foil or blade 

N 

g 

acceleration due to gravity 

ms" 2 

gH 

specific energy 

J ke ' 1 

h 

specific enthalpy 

J kg' 1 

H 

head 

m of liquid 

K 

lattice coefficient (equation 4.11) 


k 

an alternative to y (= C p /C v ) 


k. 

dimensionless specific speed 


L 

lift force on an aerofoil 

N 

M 

pitching moment acting on a foil 

N m 

M n 

Mach number (= Via) 


m 

mass flow rate 

kgs' 1 

N 

rotational speed 

rev min' 1 

NPSE 

net positive suction energy 

kJ kg" 1 

NPSE a 

net positive suction energy available 

kJ kg' 1 

NPSEr 

net positive suction energy required 

kJ kg' 1 

NPSH 

net positive suction head 

m of liquid 

N, 

specific speed 


o 

opening or throat in a turbine cascade 

m 
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P 

Po 

Pv 

P 

Q 

R 

R 

R c 

^em 

5 

t 

r 

T 

To 

T 

u 

V 

V, 

V n 

Vis or 

(^,scn) 

^R 

V'u 

w 

w u 

Yn, 

z 

a 

P 


6 


>1 

>lss 
n - rs 
'/tt 

e 

K 

P 


P 


pressure 

stagnation pressure 
vapour pressure 
power 

volumetric flow rate 
reaction (Section 1.3) 
specific gas constant 
Reynolds number 
model Reynolds number 
suction specific speed 
blade thickness 

blade passage minimum width or throat 
temperature (absolute) 
stagnation temperature (absolute) 
torque 

peripheral velocity 
absolute velocity 

axial component of absolute velocity 
normal component of absolute velocity 

isentropic velocity (equation 1.34) 
radial component of absolute velocity 
peripheral component of absolute velocity 
relative velocity 

peripheral component of relative velocity 

loss coefficients (equation 4.27) 

blade number or position 

angle made by absolute velocity 

angle made by relative velocity 

ratio of specific heats 

stagger angle 

deviation angle 

fluid deflection 

loss coefficient (equation 4.13) 
efficiency 

static to static efficiency 
total to static efficiency 
total to total efficiency 
camber angle 
elastic modulus 
absolute viscosity 
kinematic viscosity 

Markov's loss coefficient (equation 4.26) 
density 


N rrT" 2 
N m~ 2 
Nm' 2 
J s 1 = W 
m’s" -1 

kJ kg^ 1 K _I 


m 

m 


K 

K 

N 

m 

m 

m 

m 


m 


m 

m 

m 

m 

m 


s 

s 

s 

s 

s 


- 1 
-1 
-1 
-1 
- 1 


degrees 

degrees 

degrees 

degrees 

degrees 


degrees 

kgm~‘ s~ : 

kg m _1 s -1 
^ _ 1 
nr s 

kg m -3 
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Thoma’s cavitation parameter 
velocity ratio (equation 4.29) 
flow coefficient ( V.Ju ) 
specific energy coefficient 
//• rr /2fr? s (equation 4.30) 

Howell’s work done factor 
(u angular velocity rads -1 

Subscripts 1, 2 etc. indicate the point of reference. 

For a complete definition of blade terminology please refer also to Fig. 
4.2 and Table 4. 1 . 



Fundamental principles 


/ 


1.1 Introduction 

An important class of fluid machine has, as its characteristic, the transfer of 
energy between a continuous stream of fluid and an element rotating about 
a fixed axis. Such a machine is classed as a turbomachine: fans, pumps, 
compressors and turbines come into this group. Discussion is limited in this 
book to those machines where the fluid is at all times totally enclosed by the 
machine elements, so that it is controlled by passage walls. This restriction 
excludes the Pelton turbine and wind turbines. 

The machines will be categorized by flow path and by function, as 
indicated by the simple line diagrams in Fig. 1.1 of the typical machines to 
be covered. The ideal performance laws are introduced first: the discussion 
centres on the Euler equation and its applications, it being assumed that 
basic fluid mechanics and the principles of vector diagrams are understood. 
The incompressible cases are treated first, and then attention is paid to the 
problems posed by compressible considerations. Shock wave theory and 
basic gas dynamics are also taken to be understood by the reader, who is 
referred to basic texts like those by Shapiro (1953) and Rogers and Mayhew 
(1967). 


1 .2 Euler equation 

An outward flow radial machine is illustrated in Fig. 1.2. Fluid approaches 
along the suction pipe, is picked up and operated upon by the rotor and is 
discharged into the casing at a higher level of energy. The rotor has 
imparted both a velocity and a radial position change to the fluid, which 
together result in momentum changes and resultant forces on the rotor. The 
resulting axial and radial forces on the rotating system are treated later: 
present concern centres on the changes experienced by the fluid. 

In the pump in Fig. 1.2 a typical stream surface is examined which 
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Figure 1.1 Typical flow paths in machines: (a) centrifugal or centripetal; (b) mixed 
flow; (c) bulb or bowl; (d) axial. 



Figure 1.2 Typical radial machine. 


intersects the inlet edge at 1 and the outlet edge at 2. Since momentum 
changes in the tangential direction give rise to a torque and thus to work, 
moment of momentum equations for elemental areas of flow at the points of 
entry and exit will be written down. The normal fluid velocities are V nl and 
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V n2 - ^ elemental areas of flow dn, and d« 2 are examined, the moments of 
momentum entering the rotor at 1 and 2 are given by 

dM, = (p^'n\da { )V ui R i 

AM 2 = ( p \ pod^) I u2^2 

Thus the total moments of momentum are 

Mi— entering plane (1) 

M 2 = — § pV^V^RzAai leaving plane (2) 

The fluid torque is the net effect given by 
T = Mi + M 2 

T— § />V n \ Vui/^dui — J/'I'n2l / u2^2'-l < r: (1-1) 

It is assumed that V U R is a constant across each surface, and it is noted that 
J pV n Aa is the mass flow rate m. Then equation (1.1) becomes 

T = m(Y,jiRi - V u2 R 2 ) (1.2) 

The rate of doing work is wT. and since <»R is the rotor peripheral velocity u 
at radius R , equation (1.2) can be transformed to give work done per unit 
mass: 

gH = u-\\,i - u 2 V u2 (1.3) 

This is one form of the Euler equation. 

To distinguish this gH the suffix E (Euler) will be used, and the two forms 
of the Euler equation used are: 

for pumps gH E = u 2 V u2 - "i' u i (1.4) 

for turbines gH E = n t V ul - u 2 \\ l2 (1.5) 

If gH is the specific energy change experienced by the fluid and // h is the 
hydraulic efficiency, then 


for pumps 

gH 

>1 h = u 

gH E 

(1-6) 

for turbines 

gH E 

>lh = —TV 

gH 

(1.7) 


1.3 Reaction 

This concept is much used in axial flow machines as a measure of the 
relative proportions of energy transfer obtained by static and dynamic 
pressure change. It is often known as the degree of reaction, or more simply 
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as reaction. The conventional definition is that reaction is given by. 


energy change due to, or resulting from, 
static pressure change in the rotor 
total energy change for a stage 


or, in simple enthalpy terms, 


static enthalpy change in rotor 
stage static enthalpy change 


( 1 . 8 ) 


1.4 Application to a centrifugal machine 

A simple centrifugal pump is illustrated in Fig. 1.3. Liquid passes into the 
rotor from the suction pipe, is acted upon by the rotor whose channels are 
wholly in the radial plane, and passes out into the volute casing which 
collects the flow and passes it into the discharge pipe. 

The velocity triangles of Fig. 1.4 assume that the fluid enters and leaves 
the impeller at blade angles /?, and /? 2 , and that the heights F R] and K R2 are 
obtained from relations like V R = QlnDb. Applying the Euler equation, 

= U 2^u2 ~ Ml^ul (1-9) 

Also it can readily be shown from the triangles that 

gH E = \{{V\ - V\) + («; - u\) + (W\- W 2 2 )} (1.10) 

On the right-hand side (RHS) of equation (1.10), the first bracket is the 
change in fluid absolute kinetic energy, the second is effectively the energy 
change due to the impeller rotation, and the third is the change in relative 



Figure 1.3 Simple centrifugal pump. 
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Figure 1.4 (a) Inlet and (b) outlet velocity triangles for the pump in Fig. 1.3. 


kinetic energy. Equation (1.10) is thus a statement that the total energy 
change is the sum of velocity energy change (the first bracket) and the static 
equivalent energy change (the sum of the second and third brackets). 

Consider now the ideal case where = 0 (called the zero inlet whirl 
case). The inlet triangle is now right-angled and, depending on the blading 
layout, the outlet triangle can take one of the three forms shown in Fig. 1.5. 
The outlet triangles are based on the same peripheral and radial velocities, 
and they demonstrate how the absolute velocity V 2 and its peripheral 
component increase with jU. This increase, if the fluid is compressible, could 
lead to high outlet Mach numbers, and in pumps would lead to casing 
pressure recovery problems, so that the forward curved geometry tends only 
to be used in some high performance fans. 

It is instructive also to study the effect of /? : on the impeller energy 
change. So, considering zero inlet whirl, equation (1.4) becomes 

gH E = U 2 V u2 

— Z/i(W2 l^cot fi 2 ) 

Since F R2 = QIA 2 (A 2 = nD 2 b in Fig. 1.3), 

gH E = ( 1 ■ COt fin 'j It 2 

\ A 2 u 2 J 

or 

gH E = u\ - k 2 Qu 2 

for a given machine. 

Equation 1.13 is plotted in Fig. 1.6 and this figure illustrates the effect of 
p 2 upon the energy rise. It will be realized that this only relates to the ideal 
case, as the usual pump characteristic departs considerably from the straight 
line owing to friction and other effects. 


( 1 . 11 ) 

( 1 - 12 ) 

(1.13) 
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Outlet 

velocity 

triangles 


Figure 1.5 Effect of outlet angle on the outlet triangles for a centrifugal pump. 



Figure 1.6 Influence of outlet angle on the ideal centrifugal pump characteristic. 


1.5 Application to axial pomps and turbines 

, 5.1 Axial pump or fan 

. , hese idealized cases, flow is assumed to be passu,, through the machine 
arallel to the machine centre of rotation at all iadn. the amuilus dimensions 

remain constant as shown in Fig. 1.7. , 

If one plane of flow is chosen, say at radius /%. the peupherai vchvm , ,s 



Figure 1.7 Axial how 
(GV = guide vane). 


pump or compressoi stage and the ideal veloeiiv uiangk 



8 Fundamental principles 


coR and is constant through the machine. This, and the assumption that V a is 
constant, allows the triangles in Fig. 1.7 to be drawn. 

The Euler equation, equation (1.4), reduces to 

gH E = u(V u2 - V ul ) = uAV u (1.14) 


or, using the extended form of equation (1.10), 

gH E = \[{V\ - Vf) + {W\ - W\)\ 
The reaction, equation (1.8), can be expressed as 

UV ul + F u2 ) 


R 


1 - 


(1.15) 


(1.16) 


A convenient way of drawing the velocity triangles is to draw them together 
on a common base (Fig. 1.8) or on a common height, which is not so usual. 
The common base will be used throughout, and allows AV U and changes in 
velocity size and direction to be seen at a glance. 

To illustrate how the shape of the velocity triangles is related to the 
degree of reaction, Fig. 1.9 has been drawn for three cases (a) axial inlet 
velocity, (b) axial outlet velocity and (c) 50% reaction. Using equation 
(1.16), the reaction in case (a) is between 0.5 and 1, depending on the size 
of V u2 . In case (b) R is always greater than 1. The "symmetrical’ case (c) 
(50% reaction) yields R as 0.5, and has been much used in compressors 
since K, = W 2 and V 2 = WV The triangles are symmetrical on a common 
base, making them easy to draw and understand. 


1.5.2 Axial turbine stage 


An axial turbine stage consists of a stator row, usually called a nozzle ring, 
directing fluid into the rotor row. Again the assumption is that fluid is 



Figure 1.8 Axial velocity triangles based on a common base for an axial pump stage. 
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Figure 1.9 Effect of reaction choice on the velocit\ triangles for an axial compressor 
stage. 


moving parallel to the axis of rotation, and equations (1.5) and (1.10) 
become 


gH V: = u(V u] - V u2 ) (1.17) 

or 

= $[(Vt - V\) + (W; - Kt)] (1.18) 

Figure 1.10 illustrates blade layout, velocity diagrams and the common base 
velocity triangles, which may be considered typical. 
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Common base diagram 

Figure 1.10 Velocity triangles for an axial turbine stage (NGV = nozzle guide vane). 



(1.19) 


Alternative operating modes 1 1 




Figure 1.11 Effect of reaction on an axial turbine stage velocity triangle: (a) zero 
reaction. R = 0 ( />', = /?•>. IV, = W-.): (b) 50% reaction. R = 0.5 (IV, = W. VV- = 
V,). 


Two common reaction cases are sketched in Fig. 1.11, (a) zero and (b) 
50% reaction. Case (a) is the so-called ‘impulse’ layout, where the rotor 
blade passages have inlet and outlet angles the same, and are of constant 
area. Many steam and gas turbine designs use this layout. Case (b) is the 
50% reaction layout, so often used in compressors. Many turbine designs 
use low reactions of 10-20% design, as then the ‘leaving loss’ due to V' : is 
minimized. 


1 .6 Alternative operating modes 

Before going on to discuss cavitation and thermodynamic limitations on 
machine performance, a comment will be made on the possible modes in 
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Figure 1.12 Possible operating modes for a radial machine rotating in the positive 
(normal) direction. 


which a machine may be used. For the centrifugal pump sketched in Fig 
1.3, two directions of rotation may be applied: the normal for pumping, 
identified as positive, and the reverse. Possible characteristics are sketched 
in Figs 1.12 and 1.13. These possible modes are briefly discussed in Table 
1 . 1 . 


1.7 Compressible flow theory 

1.7.1 General application to a machine 

The preceding sections assumed that the density was constant throughout 
the machine concerned. If the density does change considerably it is 
necessary to relate enthalpy changes to the equations already discussed. 
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Figure 1.13 Possible operating modes for a radial machine rotating in the negative 
(reverse) direction. 


Table 1.1 Possible operating modes for a radial machine 


Mode 

Flow 

Change 
in gi f 

Rotation 

Comment 

A 

+ 

+ 

4 - 

Normal pumping mode 

B 

+ 

- 

+ 

Energy dissipation, outward flow' turbine. line 





B' in Fig. 1.12 

C 

- 

+ 

+ 

Energy dissipation, rotor resisting back flow 

D 

+ 

“b 

- 

Pump being driven wrong way round 

E 

+ 

- 

- 

Energy dissipation 

F 

— 

+ 

— 

Turbine operation 


Suppose that a compressor is considered to be the control volume (Fig. 
1.14) and that the fluid system moves through the control volume in time dr 
as sketched. If the heat added in the time interval is d Q and the work 
extracted is dVF. then application of the first law of thermodynamics to the 
system at time t and time t 4- dr in Fig. 1.14 yields the equation (following 
Shapiro, 1953) 

dQ - dH 7 = (f B — f A ) 4- (rne 2 — mey) + (p 2 >nv 2 - ppnvy) (1.20) 

On the RFIS of equation (1.20). the first bracket is the change in internal 
energy of the fluid in the control volume, the second is the difference 
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Control volume 




Time t Time f+df 

Figure 1.14 Compressible flow theory example. 


between the internal energy at exit from and inlet to the control volume, 
and the third is the net work done by the system on the surrounding fluid as 
m passes through the control volume. If flow is steady, the first bracket is 
zero and the second and third brackets can be rewritten in terms of 
enthalpy, velocity and potential energy as follows: 

dQ - dW = m[(h 2 - h { ) + (l /? - V])I2 + g(Z 2 - Z,)] (1.21) 

This can be expressed in terms of work rate and mass flow rate. The 
suppression of heat exchange results, if stagnation enthalpy is used, in the 
equation 

W = w(/j(, 2 - (1.22) 

or 

W = mC p (T iP _ - T m ) (1.23) 

In a turbomachine, 

W = m{gH) = (u 2 Vu2 - Mi^ui)^ 

Thus 


C p (T a2 — T(n) — (ii 2 V u2 ~ w iT u i) (1-24) 

or 

(hoi ~ hoi) = (u 2 V u2 ~ u { V ul ) (1.25) 

Thus stagnation enthalpy change relates to the velocity triangles. The way 
this relates to the compression and expansion processes will now be 
discussed. 
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Figure 1.15 h-s diagram for compression between two planes in a machine. 


1.7.2 Compression process 


Figure 1.15 illustrates the compression from state 1 to state 2; both static 
and stagnation conditions are depicted. Actual compression is shown as 
accompanied by an increase in entropy from S, to S 2 . The isentropic 
compression is shown as increasing stagnation enthalpv from point 01 to 
point 02s. Isentropic efficiency statements may be written as total to total 


hti2\ fhu 

~ /«,„ 


(1.26) 


or as static to static 


>1 ss = 


/7 2 S — /?, 

h 2 - h i 


(1.27) 


The choice of efficiency depends on the system in which the compressor is 
placed. //-(• i being appropriate for situations where stagnation conditions at 
inlet and outlet are proper indexes of performance. The equations apply to 
the whole process of compression but, because of the divergence of the 
constant pressure lines (Fig. 1.16), the sum of the stage isentropic rises 
exceeds the overall rise. Thus, if i] r is the small stage efficiency and ;; c the 
overall efficiency. 


>h > 


(1.28) 
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Pa 



Figure 1.16 Preheat effect in a multistage compressor. 


The preheat factor is given by 

r/ p _ £ stage rise 
// c overall rise 

The small stage or polytropic efficiency is given by 

_ d/t is _ vdp RTdp 
' lp ~ ~dh ~ C p dr _ C p pdT 

Thus 

d T (y - 1) dp 

T yn p p 

since 

Cp = 


(1.29) 


When integrated for the compressor this yields, if ;/ p is the same for all 
stages. 


h = /pgV’ 

Tl \pi) 


(1.30) 


If a; p = 1 for the ideal compression process. 


L* = (P2Y~~ lV7 
T\ \pj 
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R 

Figure 1.17 Plot of the variation of polytropic and overall efficiencies with pressure 
ratio in compression. 


Thus 


(;-!)/•/ _ 


'h = 


(PiJP t) 
(pMf~ lv ^ - i 


(1.31) 


Figure 1.17 is a plot of >/ c against the pressure ratio at various i/ 

Figure 1.18 illustrates a simple centrifugal compressor. Figure 1.19 is the 
h~s diagram relating to the compression flow path, which has been con- 


, = 

hoi — hiu 

. From this figure. 



(tlMsss 

~ hm 


>h 

' ' h, u - 

hm 

U) 


l> 4sss - 

- h i 


>h 

s h, - 

hi 

(133) 


1.7.3 Expansion process 

In compressible machines some energy transformation takes place in a 
nozzle before the fluid passes through the turbine rotor. If heat losses in a 
nozzle (Fig. 1.20) are neglected, the stagnation enthalpy remains constant. 
The isentropic heat drop is h m - h 2s = A/t iscn . If this could be utilized, the 
outlet velocity would be 

V K = V(2Ah lsc J (134) 

Since the actual Ah = /?,,, - h 2 , a commonly used nozzle efficiency mav be 
written as 

(’PI ~ (t; 


H nozzle 


(1.35) 
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Two other definitions are quoted by Kearton (1958): 

hi - h 2 


and 


(1.36) 


(^01 - h 2 ) - 7con(^01 - h 1 ) 




hi - /t2s 


(1.37) 
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Figure 1.20 Simple expansion process. 


where /; c ,, n is the efficiency with which the inlet kinetic energy is converted. 

If a complete stage is considered, the h-s diagram in Fig. 1.21 results. 
Also shown in this figure are the rotor inlet and velocity triangles, and a 
sketch of a typical annulus. The It-s diagram for the nozzle is based on 
h {) i = /to:, and the rotor h-s is based on /legislative = /fiorciuiivc' since the 
velocity vectors W\ and allow total enthalpy relative to the rotor blades 
to be calculated. It can be shown by using equation (1.25) that for an axial 
stage 

>h)2R _ "ti?R - 9 ~ 9 - d 


as assumed. Using the h-s diagram. 


«(H ~ “U3ss 


(1.38) 


Vrs = 


/’(II ~ h 03 
*01 — /l.lss 


(1.39) 


In general. >/ TS is a more realistic statement for machines in systems where 
the outlet velocity energy is not utilized, and i] TT for applications like the 
pass-out steam turbine where the energy is utilized. 

Another efficiency much used in basic steam turbine theory is the so- 
called diagram efficiency: 


work from velocity triangles 
energy available to rotor blades 


(1.40) 


The h-s diagram in Fig. 1.21 is drawn for a reaction machine layout, but it 
is necessary to comment upon the distinction between ‘impulse’ and ‘zero 
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Nozzle 



Figure 1. 21 h-s diagram and velocity triangles for an axial turbine stage. 


reaction' systems. By zero reaction it is understood that the whole of the 
stage enthalpy drop takes place over the nozzles and none in the rotor, so 
that the rotor h-s diagram results in the shape of Fig. 1.22(a). If the true 
impulse concept is applied. p 2 = p 2 . The result is Fig. 1.22(b), and an 
enthalpy increase from h 2 to /i 3 follows, so the impulse concept is strictly one 
of negative reaction. 

Figure 1.23 illustrates the 'reheat' effect which follows from the divergence 
of the pressure lines. The reheat factor is given by 


R H 



nr 


(1.41) 


where r/ T is turbine efficiency and ;; p is 'small stage’ efficiency. 

Figure 1.24 is a plot of ij t against pressure ratio for various /; p . obtained 
by an argument similar to that for the compressor. 

Figure 1.25 illustrates a radial inflow turbine, and Fig. 1.26 shows the h-s 
diagram and the ideal velocity triangles. For the turbine. 


'1 tt — 


^01 ~~ ft 03 
hf)l — h 03ss 


(1.42) 
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Figure 1.22 Distinction between the ‘impulse’ and ‘zero reaction’ concepts. 



Figure 1.23 Reheat factor. 



Figure 1.24 Effect on overall efficiency and polytropic efficiency of expansion ratio. 


htn hiii 
hoi ~ ^3ss 


0.43) 


It will be noted that it was assumed that h 0] = /r 02 , but that since peripheral 
velocities u 2 and u 2 are dissimilar, the rotor portion of the diagram takes the 
shape shown. 
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Figure 1.25 Radial inflow turbine. 



Figure 1.26 h-s diagram and velocity diagrams for a radial inflow turbine. 


1 .8 Shock wave effects 

The reader is referred to Shapiro (1953) and Ferri (1949) for shock wave 
theory, as the object here is simply to comment that when shock waves form 
they create fluid dynamic difficulties. 

In compressors, the conditions at inlet to a row typically cause shocks to 
form, and these have an effect on the flow into the passages in addition to 
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Figure 1.27 Prandtl-Meyer shock-wave effect. 


the flow instabilities caused by the wakes of the preceding blade rows in a 
multistage machine. 

In turbines, flow near the trailing edges of nozzles is of particular 
importance. If expansion ratios are selected that give rise to shock con- 
ditions, the flow is deflected by a Prandtl-Meyer expansion (Fig. 1.27), the 
misdirection caused being several degrees. Horlock (1966), for example, 
shows that a deflection of 12° resulted from a nozzle with an expansion ratio 
over the critical, with a Mach number of unity in the throat expanding to a 
back pressure half that at the throat. Clearly, a substantial correction is 
needed in design calculations. 


1.9 Illustrative examples 

1 .9. 1 Radial outflow machine (pump ! 

The pump sketched in Fig. 1.28 is driven at 1470 rev min -l and delivers 
1001s -1 with a specific energy change of 400 J kg -1 . Sketch the inlet and 
outlet triangles, assuming a hydraulic efficiency and zero inlet whirl. 

The shapes of the triangles are shown in Fig. 1.28. The calculations are as 
follows: 


'6 


1470 x 7i x 0.2 
30 x 2 


15.39 ms 1 


u 2 = 28.48 ms ! 

The radial normal velocities F R , and F R2 are found using the flow rate: 

F R i = 0.1/tt x 0.2 x 0.03 = 5.31ms -1 
F R2 = 0.1/tt x 0.37 x 0.03 = 2.87 ms -1 
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4/ 2 ( = 28-48 




u x (= 15 - 39 ) 


Figure 1.28 Radial outflow machine. 


The Euler specific energy rise is given by 


400 

085 


28.48 V 


ii : 


Thus 

From the inlet triangle. 


V u2 — 16.52 m s 1 

/?, = tan -1 5.31/15.39 
= 19.04" 


Similarly, from the outlet triangle, 

fi 2 = tan -1 2.87/(28.48 - 16.5) 


13.47° 


1.9.2 Axial pump and turbine 

The axial machine sketched in Fig. 1.29 is driven at 45 rads -1 . If the energy 
change is 120 J kg -1 , sketch the velocity triangles for both pumping and 
turbining modes of operation assuming V a = 12 m s -1 . Ignore efficiency, and 
assume zero inlet whirl for the pump and zero outlet whirl for the turbine. 
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Figure 1.29 Axi;il machine. 


Work on the mean diameter. 

Tihm n 

For the pump. 


23.63 ms 1 


Therefore 


120 = 23.63 V u2 


F U 2 = 5.08 ms 1 

The velocity triangles thus take the shapes shown in Fig. 1.30 with 

/I i = 26.9° 

P 2 = 32.9° 

For a turbine, the triangles follow as shown in Fig. 1.31 with 

Pi = 32.9° 
p 2 = 26.9° 


1 .9.3 Compressible flow problem 

A simple air turbine of the axial type has a nozzle angle of 20° referred to 
the peripheral direction. The nozzle efficiency is 90%. and the temperature 
drop over the nozzle is 125 K. Construct the velocity triangles if the rotor 
outlet angle is 30°, and suggest the air power available. Assume a rotor 
tangential velocity of 250ms -1 . no flow losses through the rotor, a flow 
rate of 4 kgs -1 , and zero outlet whirl. (Assume also that C p = 
1.005 kJ kg -1 K - '.) 

Using equation (1.34), and introducing the nozzle efficiency from equation 
(1.35), the nozzle outlet velocity is given by 

F, = V(2 x 1.005 x 10’ x 125 x 0.9) = 475.53 ms - ' 
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Outlet triangle 

Figure 1.30 Velocity triangles. 



Outlet triangle 


Figure 1.31 Velocity triangles for a turbine. 
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Figure 1.32 Velocity triangles for compressible flow problem. 


In the inlet velocity triangle (Fig. 1.32) it is now possible to calculate V 7 U| . 
and hence 

//, = 59.23° 

VF, = 189.29 m s -1 

With zero outlet whirl the outlet velocity triangle shape is as shown in Fig. 
1.32 and 

U3 = 288.68 m s~ 1 

The air power is given by 

muV u = 4 X 250 X 446.85 
= 4.468 x KFVF 


1.10 Exercises 

1.1 A radial outflow puntp has an impeller with an outside diameter of 
305 mm. an inside diameter of 75mm, and passage height of 49mm. If 
the blade inlet angle is 45° to the tangent and the outlet angle is 30° to 
the tangent, find for a rotational speed of 1500 revolutions per minute 
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the flow rate, pressure rise, torque and power needed, ignoring losses 
and assuming zero inlet whirl. Assume water is pumped. 

1.2 An axial flow pump has a rotor tip diameter 1.5 m and hub diameter 
0.75m. If the rotational speed in 45 rads -1 , the axial velocity through 
the rotor is 10ms -1 , sketch the mean diameter velocity triangles and 
estimate the Euler specific energy rise. Assume ideal flow, zero inlet 
whirl, the outlet blade angle referred to the tangential direction is 25° 
and that water is the fluid pumped. 

1.3 If the machine in 1.2 is used as a fan passing air of density 1.2kgm -3 
suggest, basing calculations on the mean section, the volume rate of 
flow and the idea! pressure rise. 

1.4 An axial flow fan has a tip diameter of 2m, a hub diameter of 0.8m, 
and rotates at 1450 revolutions per minute. For the condition of zero 
inlet whirl estimate the velocity diagrams at the tip section, if the inlet 
absolute velocity is 55 ms -1 the air has a density of 1.2 kg m -3 and 
losses are ignored. Estimate also the fluid power, if Ap is 5kNm -2 . 

1.5 A small water turbine is of the axial type, and is installed in a river 
barrage where the level drop is 8 m. The rotor tip diameter is 1.5 m, the 
hub diameter is 0.62m, the rotational speed is 45 rads -1 , the axial 
velocity through the rotor is 10ms -1 , the hydraulic efficiency is 92% 
and the machine mechanical losses plus electrical losses reduce the 
turbine output by 5%. Estimate the mean radius blade angles assuming 
ideal flow conditions, and the power output in MW. 
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Principles and practice 
of scaling laws 


2.1 Introduction 

This section covers the fundamentals of similarity, scaling, and the problems 
involved when models are used to predict full size machine performance. 
The fundamental hydrodynamic principles leading to Euler's equation are 
introduced and applied to incompressible and compressible machines, and 
the problems associated with their use are discussed. 

Three machine flow paths are possible - radial, mixed and axial flow, as 
indicated in Fig. 1.1 - and the flow direction and rotation is related to the 
machine action, either as energy input or as extraction. 


2.2 Performance laws 

In a simple approach, a turbomachine can be considered as a black box. 
Figure 2.1 describes a functional unit. A shaft transmitting torque is the 
external work output input P. flow is passing through at a volume rate Q, 
the fluid is experiencing an energy change gH , and the rotor speed is to. 
How the energy gH is achieved, and its relation to Q and P, were discussed 
in Chapter 1. The quantities listed are all measurable externally, and can be 
used to indicate machine performance in a way relevant for the user. 

If a pumping machine is considered, the input power P is a function of 
flow rate Q. specific energy rise gH. machine size (characteristic dimension 
D). rotational speed to. and liquid properties //, p and k (modulus). Using 
the principles of dimensional analysis, sets of non-dimensional groups may 
be obtained, one such being 

P _ Q gH poo!) 2 per D 2 
pc/'D ' _ tolf' crD " p k 

(1) (2) (3) (4) (5) 


( 2 . 1 ) 
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Figure 2.1 



The following may be noted for equation (2.1). 

• Grouping 1 is a power coefficient. 

• Grouping 2 is often known as a flow coefficient 4> (or, since Q VD 2 , 
and coD u, rj> = V/u, which is a velocity coefficient). 

• Since arD 2 u 2 , grouping 3 can he written as gH/u 2 , and without g is 
known as y. the head coefficient. 

• Grouping 4 is a Reynolds number based on a typical machine dimension. 

• Grouping 5 is a form of Mach number irrelevant in pumps, fans and 
water turbines, but not in other turbomachines. 

One use of the first three groupings is to present the typical characteristics, 
shown in Fig. 2.2, in non-dimensional form as in Fig. 2.3. 

They may also be used to predict the probable dynamically similar per- 
formance of the same pump running at different speeds, or of a pump in the 
same ‘family’, using the ‘scaling’ laws 

P 

pco 3 D 5 

Q 

wD 3 

gH 

o/D 2 


constant 

constant (2.2) 

constant 


Figure 2.4 illustrates the point. 
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Figure 2.2 Constant speed characteristics of a pump. 



Figure 2.3 Non-dimensional presentation of the complete performance charac- 
teristics of a pump. 


Equation (2.1) can thus be seen as a basic performance law. which should 
be valid for all fluids, compressible and incompressible, providing account is 
taken of changes in fluid properties. For compressible fluids it is conven- 
tional to express quantities in terms of inlet conditions, and instead of 
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Figure 2.4 Application of scaling laws to pump performance translation. 


equation (2.1) there is an alternative form, much used in the gas turbine 
industry in particular: 


P 02 _ j- 
Pei 


ml) 

-VRTm' 


m\Z(RT Q \) 

D 2 Pm 


,R c ,M n . .. 


(2.3) 


Conventionally the gas constant and D are not quoted explicitly (the gas and 
size being specified separately), so equation (2.3) is often written as: 


Pm _ y 
Pm 


m 

-V7V 


mV T 0 , 
Pot 


,R c .Af n . . . 


(2.4) 


and performance plots are presented in Fig. 2.5 for compressors and 
turbines. 


2.3 Concept of specific speed 


The classical view of the problem of characterizing the performance of a 
pump without involving its dimensions directly is discussed very clearly by, 
among others, Addison (1955). He defines a pump as being of standardized 
‘size' when it delivers energy at the rate of one horsepower when generating 
a head of one foot. This imaginary wheel has a speed termed the specific 
speed. He shows that the specific speed N s is given by 


jV s = K 


NVQ 


(2.5) 


where K contains the fluid density and the horsepower equivalent, N is in 
revolutions per minute (rpm), Q in gallons per minute (gpm), and H in feet 
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Figure 2.5 Non-dimensional plots of compressor and turbine performance; (a) com- 
pressor; (b) turbine. 


of fluid. It is customary to suppress K. and for many years the accepted form 
of specific speed has been 


Ay = 


NVQ 

H ?a 


( 2 . 6 ) 


Though called specific speed, the dimensions of vary with the units 
involved. N is usually in rpm. but Q can be in imperial or US gpm, litres per 
second or cubic metres per second, and H in feet or metres. Caution 
therefore is needed when A/ s is used. 

The equivalent expression relating to hydraulic turbines takes the form 


N s = 


NVP 

/T 4 


(2.7) 


Here P may be in kilowatts or horsepower, so that care is needed in reading 
literature. Neither of equations (2.6) or (2.7) is non-dimensional, and care 
has to be taken to determine the dimensions used when using them. This 
fact, and the introduction of the SI system, has given rise to the so-called 
characteristic number 

_ (oVQ 

5 " (gHf A 


( 2 . 8 ) 




Radial Mixed flow Axial 


Figure 2.6 Variation of overall efficiency with non-dimensional specific speed (after 
the well known Worthington plot). 


where co is in radians per second, Q tn cubic metres per second and gH is 
the specific energy change. This gives a non-dimensional result, as does the 
Q number (where Q and gH are in the foot, pound, second system) due to 
Csanady (1964). 

Figure 2.6 indicates how this concept may be used as an index to flow' path 
shape and type of machine in a range of applications. This figure is a well 
known chart for pumps. It illustrates how flow path changes as the specific 
speed increases. A low number means that the flow rate is low relative to 
energy change, and a radial impeller is needed; a large value suggests a high 
flow relative to energy change, and an axial rotor is required. It can readily 
be seen that specific speed could be a basis for presenting dimensions, as 
discussed by Balje (1981), for example. 

2.4 Scale effects in incompressible units 

2.4.1 Hydraulic machines 

The use of models for predicting full size water turbine performance is well 
established, but the same technique for pumps is relatively recent, having 
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Figure 2.7 Non-dimensional presentation of model and full size performance for the 
Eggborough pump (adapted from Nixon and Cairney. 1972a). 


been made necessary by the increasing size of cooling water pumps for 
thermal power stations. Significant differences in efficiency between model 
and full size occur, suggesting departures front strict dynamic similarity. It is 
argued by both turbine and pump authorities that losses differ, and that 
work/capacity curves differ too. All are agreed that exact mathematical 
treatment is not possible, and attempt empirical formulation from ex- 
perimental data. The references cited, which show wide variations, may be 
consulted as an extension of this section. The discussion here is based on the 
reports by Nixon (1965), and Nixon and coauthors (1972a. b). which are 
admirable examples of a logical approach. 

As these authors suggest, departures from the scaling laws are due to: 

• geometrical dissimilarities due to tolerances 

• clearance variations 

• surface finish 

• hydrodynamic effects 

• testing errors 

• installation effects. 

These were examined for the prototype and four models, three of which 
were in aluminium and the other in fibreglass, but being one-eighth scale for 
comparison with the aluminium surface produced. Figure 2.7 illustrates the 
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variations observed from the prototype machine when the three aluminium 
models were tested. 

Geometrical dissimilarities and their effects were examined. For example 
blade shape, as determined by blade angle, was compared. A variation of 
±0.5° about a mean value which was 1.5° less than design was found for the 
full size machine, and in the models the scatter was of the same order, apart 
from the glass-fibre impeller in which it was much worse. It must be com- 
mented that much higher variations are quite common in commercial cast 
impellers. Since the small aluminium models were shell mouldings the area 
variations were small, but the larger one was floor moulded, giving at least a 
10% variation in area from passage to passage. Again, larger variations have 
been observed, particularly in cast impellers. The volute throat is a most 
important area, and correction is needed if variations occur. 

Clearances tend to be of the same order in models as on the prototype, so 
leakage losses tend to be larger, and model surface roughness cannot be 
supersmooth in scale without large cost. This is therefore a significant 
problem because it affects boundary layers, as discussed in the following 
section. 

Turning to hydrodynamic problems, Nixon and Cairney (1972). Oster- 


walder (1978). 
relation: 

and Osterwalder and Ettig (1977). 

suggest the 

following 

II 

H 

II 

> 

1 

<5 m + 0, + 

<5d + 



total 

mechanical leakage 

disc 

skin 

inertia 

loss 

loss loss 

friction friction 

loss 

(2.9) 


Here, <5 N1 and <5] are unaffected by the Reynolds number and <5, is usually 
assumed to remain the same. c5 M is considered to vary as speed, in contrast 
to the other hydrodynamic losses which tend to follow an N 3 law, and 
reduces with reducing speed at a lesser rate, thus being proportionally more 
important at low speeds. Nixon and Cairney (1972) present a method of 
finding <5 M , and suggest that prediction from low speed tests be limited to 
differential head readings. 

The estimation of disc friction loss has been a subject for argument, as the 
classical work was done on plain thin discs rotating in a close fitting closed 
casing. Nixon used work by Necce and Daily (1960) and Watabe (1958) for 
"smooth' and ‘rough’ discs, and showed an error from measured data of 
about 10%. Sutton (1968) studied this problem, particularly the effect of 
leakage flow through wear rings and its relation to disc friction. Osterwalder 
(1978) commented that there is little current data of general applicability, 
but Kurokawa and Toyokura (1976) and Wilson and Goulburn (1976) 
extended the database. 

The same situation is found in attempting to correlate <) F . Both Nixon and 
Osterwalder suggest the applicability of Nikuradse and Colebrook data with 
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Table 2.1 A selection of model scale formulae (as quoted for example by Nixon 
(1965). Subscript m denotes ‘model’ 


Moody (1942): 


Moody: 


Anderson: 


Pfleiderer: 


1 - 9 
1 - '7 m 


Hutton: 


Ackeret: 



1 ~ i ( H A" m ( D A" 

i - » /m \h) \D I 


1 - // _ 0.94 - Q~" 
1 - », m _ 0.94 - 


, R sm y M (2 m y- 2S 


valid between 1/12 < R cm /R c < 20 



1 ~ >1 
1 “9m 




a limiting Reynolds number criterion for 'transition'. Nixon proposes a 
relatively simple approach which is applicable and gives reasonable ac- 
curacy. Osterwalder surveys the published material and does not propose a 
general correlation, but indicates that a computer-based in-house study is 
relevant. 

Nixon (1965). illustrating the deficiencies of a number of formulae in 
predicting pump efficiency from a model (some of which are in Table 2.1), 
demonstrated that there was little correlation. Anderson (1977) preferred a 
‘design’ approach using a diagram based on data from many pumps and 
turbines to estimate probable efficiency. 

2.4.2 Fans and blowers 

The prediction of fan performance from models was studied at the National 
Engineering Laboratory (NEL), UK, resulting in a report by Dalgleish and 
Whitaker (1971) in which work on three small fans was used to predict the 
performance of fans at a scale of 1.5:1 and at a scale of about 2.5: 1. The 
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report underlined Nixon’s comments on tolerancing and proposed a formula 
for //a similar to the pump equations just discussed: 

(1 ~ >/a) p = Q 3 + 0 7 (g^i) (2.10) 

^/A/m 

where >] A is the air efficiency and R c = proD 2 lp for a range 0.8 X 10 6 < R c < 
6.5 X 10 6 . They comment that clearances are important, and that the best 
efficiency moved to higher flow coefficients than suggested by scaling as size 
increased. It is of interest to note that improvement in surface roughness 
does not give better efficiency except for small high-speed fans. 


2.5 Scale effects in compressible machines 

Equation (2.1) was quoted as applying to compressible machines, but limits 
of application apply, as for incompressible machines. If single-stage ma- 
chines are considered, the effect of compressibility may be neglected for low 
Mach numbers (below about 0.5); the divergence caused increases with 
Mach number. Work with compressors using refrigerant and other 'heavy' 
gases indicates that the effect of the adiabatic exponent may be neglected 
for few stages. However, as the overall pressure ratio and hence the number 
of stages increase, density change particularly is important. As an example, 
if a pressure ratio of 8: 1 is chosen, with air ( k = 1.4) the density ratio is 
4.41:1; if k is 1.05 this becomes 7.23:1. This is clearly important in a 
multistage compressor where the stages are designed assuming a constant 
flow coefficient for air, but the machine is used for a heavier gas ( k = 1.05) 
even though use of the gas reduces power demand for dynamically similar 
conditions since acoustic velocity goes down. The only technique proved 
satisfactory is to blend gases to give the right k, as Csanadv (1964) showed. 

The scaling problems in pumps, discussed in the previous section, occur in 
compressors. For example, Henssler and Bhinder (19“7) studied the in- 
fluence of size on a family of small turbocharger compressors; these had 
60mm, 82 mm and 94 mm impeller diameters, and other dimensions were 
scaled, but surface finishes had the same surface roughness. They show how 
performance varied with flow coefficient and Reynolds number based on 
peripheral speed. The authors comment that these changes are not predicted 
by the similarity laws, but do not attempt to suggest the correlation that 
different laws apply to different families of machines. 

A contribution by Miller (1977) surveyed earlier correlations of scaling 
predictions for axial compressors. His approach was based on the need to 
conserve rig power consumption by testing at reduced pressure levels and 
dealing with the associated problems of correcting efficiency pressure ratio 
and flow to normal operating level, and was also concerned with scaling 
model tests to full size performance. Miller examined a number of Reynolds 
number correction approaches and concluded that although careful testing 
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would yield effective prediction for one compressor design, this could not be 
applied to another, and that pressure level effects appear to be more 
pronounced than scale effects. 

Others have contributed to the discussion as part of their study; for 
example, McKenzie (1980) shows a scatter in his predictions. No completely 
satisfactory general prediction appears to be available at the moment, 
although individual companies and groups use their own approaches and 
satisfy their own needs. 


2.6 Illustrative examples 


2.6.1 Similarity laws applied to a water turbine 


The turbines in a river barrage hydroelectric plant are designed to give 
55 MW each when the level difference is 25 m and they are running at 
94.7 rev min -1 . The designed overall efficiency is 93%. and the runner 
diameter is 6 m. A model with a runner diameter of 300 mm is to be tested 
under the same level difference. Suggest the probable rotational speed, flow 
rate, efficiency and power produced when the model is operating in dy- 
namically similar conditions. 


The full size flow rate is first calculated: 


55 x 10 h 
0.93 


25 g x 10 3 x Q 


Therefore 

Q = 241 mV 

Applying the scaling laws, equations (2.2), 

Q 

= constant 


gH 


= constant 


N 2 D 2 

and substituting the known data, it is found that 

£?modci = 0.602 ms -1 
^modci = 1894 rev min -1 

The model efficiency must now be found using one of the equations in Table 
2.1. The well known turbine equation due to Hutton will be used: 

. 0 .3 + 

1 - 9 model V 94.7 x 6 
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from which the model efficiency is found to be 

>7model ~ 89.77% 

The model power developed, assuming the same mechanical efficiency, is 
therefore 

P model = 25 g x 10 3 x 0.602 X 0.8977 
= 132.54 kW 

2.6.2 Compressor performance prediction problem 

A compressor for hydrogen duty is to deliver 18 kgs -1 while increasing the 
pressure from 1.01 x 10 5 Nm -2 to 16.5Nm -2 . The inlet temperature is 
expected to be 300 K, and the rotational speed 2900 rev min -1 . 

For development purposes a half-scale machine is to be tested using air as 
the medium, with inlet conditions 10 5 Nm~ 2 and 288 K. Suggest the model 
mass flow rate, delivery pressure and rotational speed for dynamical 
similarity. 


The full size pressure ratio is 16.5/1.01 = 16.34:1. The model outlet pres- 
sure will be 16.34 x 10 s Nm' 2 . Following equation (2.3), for similarity, 


mV(RT m ) 

PmD 2 


= constant 


and 


ND 

V(RT m ) 


= constant 


For hydrogen, R = 4.124 kJ kg 'K *; for air, R= 0.287 kJ kg ‘K '.Thus 

18V(4.124 x 10 3 x 300) mV(0.287 x 10 3 x 288) 

1.01 x 10 5 x D 2 ~ 10 5 x (0.5D) 2 


m = 17.24 kgs 1 


and 


2900 P N x 0.5P 

V(4.124 x 10 3 x 300) ~ V(0.287 x 10 3 x 288) 

N = 1499.15 rev min -1 


2.7 Exercises 

2.1 A centrifugal pump rotates at 185 rads -1 and at best efficiency has a 
pressure rise of 4.5 x 10 5 Nm -2 when pumping water at the rate of 
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0.28 m 3 s -1 . Predict the corresponding best efficiency flow rate and 
pressure rise when rotating at 80% of the design speed. If the efficiency 
is 85% in both cases estimate the input power required. 

2.2 If a pump geometrically similar to that in Exercise 2.1 is increased in 
diameter to give the same pressure rise of 4.5 x HPNm -2 when 
rotating at 148 rad s^ 1 suggest the scale ratio, and the flow rate and 
power input if the efficiency is 85%. 

2.3 A multistage pump is to lift water at a rate of 0.03 m' s~' from a mine 
820m deep. If the motor speed is 2900 rpm, estimate the number of 
stages needed if the specific speed in 20 (rpm nr s~ '.m). (Remember 
that the specific speed applies to one stage only in a multistage 
machine.) 

2.4 The specified duty for a hydraulic turbine is 37.5 MW at 90 rpm under a 
head of 18 m with an efficiency of 93%. The factory test bay can give a 
head of 6 m and absorb 45 kW. Suggest the necessary scale ratio for a 
model, its rotational speed and the flow rate. Assume that the model 
efficiency is 93%. 

2.5 A one-third scale model of a ventilating fan is to be constructed. The 
machine to be modelled is to run at 1450 rpm, with a duty flow rate of 
5nv's~' and pressure rise of 450Nm~ : and overall efficiency of 78% 
when the intake conditions are 5°C and 0.989 x HflNirT 2 . The model 
fan will draw air in at 20°C and l(TNm -2 . Predict the model perfor- 
mance and rotational speed when performing under dynamically similar 
conditions. 

2.6 A water turbine is to develop 100 MW when running at 93.7 rpm under 
a head of 20m with an efficiency of 92%. The runner diameter is to be 
7m. If a water model is to be tested and the maximum test head and 
flow available are 5m and 0.9m 3 s~’ estimate the scale ratio, the model 
runner diameter and rotational speed and the power to be absorbed. 
Discuss the possibilities of using air models to studv such a machine. 

2.7 A new design of compressor is to be tested. It is to deliver helium with a 
pressure ratio of 12:1 when rotating at 3000 rpm with a mass flow rate 
of 12000kgh~'. The inlet conditions are lObar and 300K. For in- 
house testing air is to be used with intake conditions of 1 bar and 300 K. 
Suggest the pressure ratio, mass flow rate, and rotational speed to give 
dynamically similar conditions and discuss the validity of your approach 
and the use of air to model other gases. 
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2.8 A gas turbine tested with air passed 15 kgs -1 when the expansion ratio 
was 1.6:1, the rotational speed was 6000 rpm, and the inlet static con- 
ditions were 10 X 10 4 Nm -2 and 288 K. The same machine when used 
for its designed purpose was supplied with combustion products 
at the inlet conditions of 1.4 x 10 6 Nm -2 and 660 K. If the expansion 
ratio was also 1.6:1, estimate for dynamically similar conditions the 
rotational speed and mass flow rate. (For combustion products R = 
2291 J kg -i K _1 .) 


3 


Cavitation 


3.1 Introduction 

Cavities form in a flowing liquid in an area where the local pressure ap- 
proaches vapour pressure; the entrained and dissolved gas form bubbles 
which grow and then collapse. In a venturi section of a pipe. Fig. 3.1, 
bubbles form near the throat and are swept downstream to an area where 
the pressure is higher and where they then collapse. There are a number of 
contributory factors: liquids absorb gas into solution (Henry's law) and the 
capacity to retain gas in solutions is a function of the local pressure, thus gas 
is liberated as pressure falls. Gas and vapour nuclei act as centres of growth 
and liberated gas helps bubble growth. (Henry's law' suggests that w'ater will 
absorb air through a free surface at atmospheric pressure at the rate of 2% 
of the water volume; hydrocarbon fluids have a greater capacity to absorb 
gas.) Bubble growth and collapse is rapid, as Figs 3.2 and 3.3 indicate. It has 
been observed that bubbles can burst and reform several times, indicating 
that the process is unstable. This gives rise to flow instability and in ma- 
chines the result is noise, vibration and surface damage. 

Since cavities form where the pressure is low. cavitation is a problem in 
the suction zone of pumps and in the area of discharge from a water turbine. 
Figure 3.4 illustrates the effect on the characteristics of a pump of reduction 
in the suction pressure leading to cavitation, and Fig. 3.5 illustrates the 
effect of cavitation on the behaviour of a turbine. 

In a centrifugal pump. Fig. 3.0. the rotational effect of the rotor is 
increasingly imposed on the liquid as it approaches the impeller blades. This 
results in an increasing absolute velocity with a reduction in pressure which 
continues until the leading edges of the impeller blades are reached. 
Usually, the pressure then increases on the pressure side of the blades but 
reduces further on the suction surface with bubbles forming if the vapour 
pressure is approached. Flow' in the blade passages is affected by the 
formation and collapse of the cavities and this gives rise to a reduction in 
performance, as sketched in Fig. 3.4. and also to noise and vibration with, in 
addition, damage to surfaces (section 3.3). 
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Figure 3.2 The life history of a bubble when growing and bursting in regions at two 
different pressure heads. 


3.2 Net positive suction energy (NPSE) or head (NPSH) 

Hydraulic machinery engineers have used an expression for a long time 
which compares the local pressure with the vapour pressure and which states 
the local pressure margin above the vapour pressure level. In the era before 
the SI system it was always known as the net positive suction head (NPSH) 
and must now be known as the net positive suction energy (NPSE). Two 
definitions are used, that which relates to the pressure presented to the 
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size 



Figure 3.3 The life history of growth and collapse for a single cavity over 6 micro- 
seconds (based on data in Knapp and Daily. 1970). 



Figure 3.4 The effect of net positive suction head on a pump characteristic. 


machine by the system (NPSEJ and that generated by the dynamic action of 
the machine (NPSE R ). 

3.2.1 NPSE available (NPSE.,) or NPSET, 

There are two main suction systems applied to pumps: the drowned system 
and the suction lift system (Fig. 3.7). The lowest pressure P s in the suction 
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Figure 3.5 The effect of cavitation inception on a water turbine 
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Figure 3.6 The variation of liquid pressure in the suction region of a pump 
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Figure 3.7 Sketches of typical suction systems: (a) drowned system; (b) suction lift 
system. 


line will occur at the suction flange. NPSE a is then defined by the equation 

NPSE A = PJp-PJp (3.1a) 

where P y is the total pressure of the fluid and P v the vapour pressure at the 
temperature prevailing. In the older system 

NPSH a = total head - vapour pressure head (3.1b) 

For the system in Fig. 3.7(a) 

NPSE a = P!p + gH - flow losses in the suction line — PJp (3.2) 
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Figure 3.8 Variation of vapour pressure of water with temperature. 


and. for the system in Fig. 3.7(b), 

NPSE a = P/p — gH — flow losses in the suction line — PJp (3.3) 
where P is absolute pressure. 

The flow losses must be calculated for the flow rate that the pump is to 
deliver. An example of a calculation for the system in Fig. 3.7(b) now 
follows. 

A water pump has its centreline 5 m above the free surface in a sump 
which is at atmospheric pressure. The water temperature is 20 °C. If the flow 
losses are calculated to be lOJkg -1 estimate the NPSE a . From data sheets 
and Fig. 3.8, P v is 2.354 x 10 3 . Thus, 

in 5 2 354 

NPSE a = ^3 - 5g - 1° - -~r = 38.596 J kg- 1 (NPSH a = 3.93m) 

As an example of an NPSE a calculation for a drowned suction, consider a 
tank that has a free surface 2.3 m above the pump centreline, containing n- 
butane which is held at its boiling point of 37.8 °C. Friction losses are 
computed to be 5 Jkg -1 at the flow condition under study - estimate NPSE a 
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if the vapour pressure for n-butane at 37.8 °C is 3.59 bar, and the relative 
density is 0.56. 

Since the liquid is "boiling’ the free surface is at 3.59 bar, therefore 


NPSE ; , = 


3.59 x 10 5 „ _ 

0.56 x 10 3 + 


s — 


3.59 x 10 5 
0.56 x 10 3 


1 9.525 J kg“' 


or, NPSH a is 1.99m of butane. 

If the suction pressure is obtained from a pressure gauge mounted in the 
suction flange instead of being estimated from the system data, it has to be 
remembered that the gauge reads static pressure, so the total pressure is 
obtained by adding the kinetic term. 


P = l w + Yl 

P P 2 


(3.4) 


3.2.2 NPSE required (NPSE R ) 

This is the statement of the margin above vapour pressure that the pump 
can cope with by its dynamic action. It is clearly a function of the impeller 
geometry and is very difficult to estimate with any accuracy. The classic text 
by Pfleiderer (1961) contains equations that will given an estimate, but in 
general the approach has been to use empirical data to give NPSE R . The 
most common approach has been to use the Thoma number, first developed 
for water turbines: 


_ NPSE R 
Machine gH 

or 

c t = NPSH R /machine head rise 

Figures 3.9 and 3.10 are plots of o against characteristic number k s , based on 
well known diagrams to be found in the classic texts. This approach gives 
conservative values of NPSE R , and yields values which are based on design 
flow and energy rise. 

3.2.3 Critical or limiting NPSE 

The two statements of NPSE are based on different criteria, so their relation 
to one another must be discussed. If a pump is considered, NPSE., and 
NPSE R can be plotted against flow rate (Fig. 3.11). NPSE a will fall with 
increasing flow rate since the flow losses vary as (flow' rate) 2 ; NPSE R varies 
in the manner shown. The point at which the curves cross is called the 
critical flow, as at higher flow rates the pump cannot sustain the suction 
pressure presented to it, and at lower rates the pump is supplied with 


(3.5a) 

(3.5b) 
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Figure 3.9 Variation of Thoma's cavitation parameter with characteristic number for 
pumps. 


much larger suction pressures than it needs and cavitation should not be a 
problem. It can be argued that at flow rates higher than the critical value the 
pump will cavitate badly, and at lower flow rates it should behave satisfac- 
torily. While this is a general rule, studies over the last 15 years have 
supported the view put forward by Bush et al. (1976) that at low flow rates 
recirculation in the suction line and between impeller and casing can give 
rise to cavitation-induced damage and instability (section 3.4). 

In general, pumps are tested for their cavitation performance by reducing 
the NPSE whilst maintaining the pump speed and flow constant. The result- 
ing curve (Fig. 3.12) shows a level energy rise until at low suction pressure it 
falls away. The point at which the performance is stated to be unacceptable 
is when the energy rise falls by x% (usually 3% for standard pumps as laid 
down by the American Hydraulic Institute). The corresponding NPSE is 
called the critical value. 

Hydraulic turbines, when model tested, are usually tested at the design 
full size head drop and the point at which cavitation is observed noted. In 
turbines the effect is not so marked as cavitation is downstream and air is 
injected into the draft tube to reduce pulsations. Damage is. of course. 
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Figure 3.10 Variation of Thoma s cavitation parameter with characteristic number 
for Francis anti Kaplan water turbines. 



Figure 3.11 NPSLi available and NPSE required curves presented on a base of flow 
rate to illustrate the critical Mow rate where the curves cross. 
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Figure 3.12 The usual way of testing and presenting pump cavitation behaviour 
(The test is performed of constant speed and constant design' flow.) 


important in all machines so effort is expended to reduce cavitation where 
possible. 

As an example of NPSE R calculation consider a pump with a design duty 
of 35nv'h _l when delivering against a specific energy rise of 409 Jkg~' and 
rotating at 3000 rpm. The system NPSE is computed to be 30Jkg _1 . 
Comment on the risk of cavitation. 


30007T 1/ 35 \ 



from Fig. 3.9 o = 0.05. thus 

NPSE r = 0.05 x 490 = 24.5 Jkg~' 
or, NPSH r = 2.5 m of liquid. 

NPSE a is greater than NPSE r by a small margin. This indicates that if cold 
water is the liquid that there may be little difficulty, if the suction system 
does not ingest air. If the fluid is volatile and is bubbly, there will not be 
enough margin and some thinking is needed about the levels of tank and 
pump and the suction system. 


3.3 Cavitation damage 

Surface damage is a serious problem in pumps and turbines. Many studies 
have been made over the years of the causes and of the ways in which 
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Figure 3.13 The jet collapse model for damage due to cavitation as illustrated by 
Lush (1987). 


damage can be reduced. Lush ( 1987a. b) summarized findings about the way 
the bubble collapse gives rise to very high local pressures on the surfaces of 
the machine channels. The existence of a microjet is proposed which, as the 
cavity collapses, bursts across the void of the bubble to hit the opposite 
surface at very high velocity (Fig. 3.13). The impact results in very high 
stresses which are equal to or higher than the ultimate strength of the 
material in many cases. There is a debate about the mechanism that leads to 
surface damage, but there appears to be agreement that the material loss is 
due to a combination of mechanical, chemical and electrolytic actions. A 
suggested mechanism, based on the computed stresses being of the order of 
2000 bar in some cases, is that work hardening with attendant temperature 
rise occurs in the material at the surface. This can lead to small changes in 
chemical composition and a spongy subsurface that leads to cracking of the 
hardened surface with subsequent penetration of the material. Eventually, 
the molecular bonding breaks down and erosion and corrosion occurs. 
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Table 3.1 Relative losses of material under comparable conditions (obtained in a 
venturi test device using water at 20 °C) 


Material Relative 

volumetric loss 


Stellite 1 

Cast stainless steel: 12.88% Cr. 0.17% Ni, 0.43% Mn, 0.38% Si 7 

Stainless steel 18:8 Cr:Ni 5 

Monel 16 

0.33% carbon steel 37 

14% Cr stainless steel (forged or drawn) 98 

Manganese bronze 118 

Gun metal 230 

Cast iron (as cast without skin) 374 

Typical cast aluminium alloy 1176 


Typical damage is described by Tillner et al. (1990) and Pearsall (1978) who 
show ways that damage can be produced in the laboratory for test and 
development purposes. The extent of the damage suffered depends on the 
fluid, the materials and the hydrodynamic system, and it has been found that 
even with advanced material loss the machine has developed the duty 
required and damage has only been found during routine maintenance. In 
pumps repair is usually by replacement, but in hydroelectric plant it is a 
routine procedure to deposit metal in damaged areas and then to return the 
surface to the high finish required. Table 3.1 summarizes the resistance of 
common materials used in hydraulic machines. 


3.4 Hydrodynamic effects 

The effect of bubble formation and collapse have been discussed briefly in 
the introduction to this chapter. The effect of reducing suction pressure was 
illustrated in Fig. 3.4, with the flow range being progressively restricted as 
suction pressure falls. Pearsall (1973) discussed this and proposed ways of 
reducing the effects by increasing the suction diameter. While this is a valid 
approach care is needed to ensure that the suction diameter does not get too 
large as then recirculation can be set up at flow rates close to the design flow 
and cause instability; reference may be made to Bush et al. (1976) and to 
Grist (1986). The studies by Pearsall and his team at NEL (1978) demon- 
strated a close relation between noise and damage and they proposed a 
noise criterion for cavitation prediction. Wolde et al. (1974) described the 
problems to be solved in noise studies and demonstrated the validity of the 
technique. 
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The effects on water turbines have been discussed by Knapp and Daily 
(1970) and by Klein (1974) who with others at the same conference illus- 
trated the effect of cavitation on the performance envelope. 


3.5 Thermodynamic effects on pump cavitation 

It is well know n that the cavitation performance of a pump varies with the 
fluid condition. For example, when pumping water the required NPSE is 
highest when passing cold water and decreases as the temperature rises. A 
similar effect is noted when other liquids are being pumped. The change is 
too large to be explained by the so-called Reynolds number effects. An 
empirical approach to this problem is known as the B-factor or /1-factor 
method outlined by Knapp and Daily (1970) and Stahl and Stcpannof 
(1956): the technique correlates suction hydraulic behaviour with vapour 
volume in the cavitating region. 

When examining bubbles in cold water it has always been assumed with 
some justification that all energy terms involving the vapour in the cavities 
are negligible compared with those of the surrounding liquid. When the 
temperature increases such an assumption may not be applied since the 
latent heat required to supply vapour to the cavity can no longer be ignored, 
and neither can the energy exchange during the expansion and contraction 
of the bubble. If they are neglected and the NPSE calculated using cold data 
the resulting value will be very conservative. As a consequence, boiler feed 
pumps in the 1950s began to be uneconomic in size as both pressures and 
temperatures rose with turbine steam conditions; eventually designers 
developed techniques to design for lower cavitation numbers than those for 
cold duties. 

The heat required by vaporization must come from the liquid surrounding 
the cavity, thus causing a drop in temperature and vapour pressure in the 
immediate vicinity of the bubble. This has the effect of reducing the bubble 
size from that which would apply in cold liquid, thus reducing the effect on 
the flow. This reasoning has led to the approach of Stahl and Stepannof 
(1956). Figure 3. 14 shows the conventional NPSE variation curves for a cold 
and a hot liquid. The two critical points where the 3% head drop applies are 
B and C. Point B is known for the cold test, and the NPSE reduction can be 
found to estimate point C. 

The v olume ratio B or fi is defined by Stahl and Stepannof as 

_ volume of vapour 
volume of liquid 

they showed that using the Claperyron-Clausius equation 

S = A(NPSE)C p t(-y-) (3.7) 
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Figure 3.14 The variation on NPSE required of liquid temperature. 


or 


B’ = 


B 

ANPSE 



(3.8) 


Figure 3.15, taken from Stahl and Stepannof (1956). plots B' for a number 
of fluids based on refinery pumps of the double suction design and a 3% fall 
in gH. Also shown are lines of NPSE adjustment in the relation 


— NPSE adjustment 

^corrected O w-V) 

oH 

1 pump 

The method is based on the assumptions that the cavities are uniformly 
distributed across the flow cross-section and that there is the 3% drop 
criterion. There are many other approaches but these indicate that there is a 
considerable difference in opinion (see for example Hutton and Furness. 
1974). 


3.6 Inducer/pump combinations 

In some applications the suction pressure is too low for the pump and it 
cavitates. Where it is not possible to change the system design there is a 
solution, an inducer, which is an axial flow stage that is placed just before 
the impeller (Fig. 3.16). The device is a high specific-speed machine and is 
required to give a small pressure rise that raises the pressure at inlet to the 
main impeller to a level that it can sustain without cavitation. It has long 
passages so that bubbles that form burst before leaving the inducer, so 
that any damage is confined to the inducer blades. The blade shapes are 
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Figure 3.15 Plot of thermodynamic adjustment factors (from Stahl and Stepannof, 
1956). 

archimedean spirals with a very acute inlet angle referred to the tangential 
direction. 

The improvement in NPSE performance that an inducer can produce is 
illustrated in Fig. 3.17. Work has shown that the improvement is only in the 
flow range close to the design point, and away from that zone, because the 
incidence angles are higher, the impeller inducer combination can give 
worse NPSE values than the impeller alone. Contributions by Pearsall 
(1978), Susada et al. (1974), Lakshminarayana (1982), Turton (1984), 
Tillner et al. (1990) and Oschner (1988), may be consulted for more informa- 
tion on both the design and behaviour of inducers when applied to pumps. 
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Figure 3.16 Sketch of a pump/inducer design. 


3.7 Concluding comments 

The undesirable effects of cavitation, flow instability, noise, and surface 
damage have been described, and basic guidelines for avoidance have been 
introduced. For detailed discussion more advanced material may be found in 
the references cited and in the textbooks listed. 


3.8 Exercises 

3.1 A pump is sited 5.5 m above a water sump, and the friction energy loss 
at design flow is computed to be 8.5 kg -1 . If the water temperature is 
25 °C compute the NPSE A at sea level and at a level where the atmos- 
pheric pressure is 0.85 x 10 5 Nm~ 2 . 

3.2 A water pump is supplied by a tank with a free surface 4m above the 
pump centreline. At the design flow the suction system loss is 50 J kg -1 . 
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Figure 3.17 Characteristic of a pump fitted with an inducer. 


If the water is at 25 °C determine the NPSE A when the free-surface 
pressure is ltPNm -1 and 1.3 x lCPNm -2 . 

3.3 If in the system described in Exercise 3.2 a limiting NPSE of 100 J kg 
is required by the pump and the losses rise to 100 J kg -1 , calculate the 
necessary pressure on the free surface. 
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3.4 A system to deliver «-butane has a tank with a free surface at 5.5 n 
above the pump. Calculate for a suction system loss of 50Jkg _1 the 
NPSE a when the suction tank free surface is at 3.5 bar gauge, and wher 
it is at vapour pressure. Assume fluid density is 0.56 kg m -3 , vapou’ 
pressure is 3.59 bar absolute, and atmospheric pressure is 0.95 > 
10 5 N m~ 2 . 

3.5 If the tank level in Exercise 3.4 falls to 4 m, and the free surface is a 
vapour pressure determine the maximum suction system loss that can bt 
allowed if the NPSE A is not to fall below 10 J kg“ 1 . 

3.6 A pump is to deliver water at the rate of 80m 3 h -1 against a systen 
resistance of 50 m when rotating at 1450 rpm. The losses in the suctioi 
line are estimated to be 25 J kg -1 at the specified flow rate. Assuming ; 
suction lift suggest the maximum level difference from sump level t< 
pump that can be sustained. Assume the water is at 20 °C and th' 
margin between NPSE R and NPSE A is to be a minimum of 10 J kg -1 . 



4 


Principles of axial 
flow machines 


4.1 Introduction 

Axial flow pumps, fans, compressors and turbines have, as a common 
feature, blades which interact with the fluid involved. Typical machine flow 
paths are shown in Fig. 4.1, and in the two-dimensional approach that will 
be discussed it is assumed that the stream surfaces are symmetrical with the 
axis of rotation. Both the isolated aerofoil and cascade data are described, 
with the limitation just outlined. The three-dimensional problems are then 
introduced, and the concepts of radial equilibrium, the actuator disc ap- 
proach. and Howell’s work done factor are related to the common secondary 
flow problems. 

The basic equations and shapes of the ideal velocity triangles and their 
relation with the Euler equation and the concept of reaction were discussed 
in section 1.2. Figure 4.2 illustrates a set or cascade of blade profiles, and 
relates blade angles to fluid angles by introducing incidence and deviation. 
The figure also demonstrates the difference between the arrangement of 
blades used in a compressor, which increases the fluid energy, and in a 
turbine, which extracts energy. Also illustrated are the other common terms 
used to describe rows of blades, and the conventional definitions of the 
coefficients in general use are defined in Table 4.1. 


4.2 Wing theory 

It is assumed that individual blades in a machine interact with the fluid as if 
they were wings, and the presence of other blades is ignored. This allows lift 
and drag data obtained from wind tunnel research on wing sections to be 
used. 

The principles relating life and drag forces with the velocity directions for 
a static blade will first be discussed using Fig. 4.3. It is assumed that the lift 
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Figure 4.1 Simple axial machine stages. 


Table 4.1 Definition of terms used with cascades of blades (to be read in conjunction 
with Fig. 4.2) 


Term 

Static 


Moving 

Fluid inlet angle 

a. 


IU 

Fluid outlet angle 

a 2 


h 

Blade inlet angle 

aj 


K 

Blade outlet angle 

a 2 


P'2 

Blade camber angle 

0[= (ai 

- «))] 

[=(/? i - im 

Stagger angle 

v or ^ 



Deflection 

4= <«i 

- a 2 )] 

[= (lh - lh)\ 

Incidence angle 

i 



Deviation angle 

(> 



Chord 

c 



Distance to point of maximum camber 

a 



Spacing or pitch 

s 



Blade height 

h 



Space : chord ratio 

s/c 



Solidity 

c/s 



Aspect ratio 

hie 




Lift coefficient C L = Ll\pV 2 (projected area) 
Drag coefficient C D = DlkpV 2 (projected area) 
Pitching moment M = LX 


The projected area for a blade is usually (c x unit height). The reference direction used here is 
the axial or normal, but often, particularly in pumps, fans and water turbines, the tangential 
direction is used. 
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(e) 

Figure 4.2 (a) Typical cascades of blades: (b) static cascade of blades; (c) forces on a 
blade: (d) definition of pitching moment: (e) definition of angle of attack. 


force /. acts at right angles to, and drag D along, the mean direction. Thus 
the tangential force on the blade may be written as 

F { = Lsin/f m + Deos/f m 


and the axial force as 


F, = L cos fi m ZD sin ji m 
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Figure 4.3 Relative velocity triangles and forces acting on a static blade. 


Since the lift coefficient C L = Ll\pV 2 (projected area) and the drag coef- 
ficient C D = Dt\pV 2 (projected area), and the projected area is the blade 
chord times the unit height, 

n K; /C L sin/? m + C D cos/? m \ /4 1N 

Ft = />yc( -jrj- ) (4.1) 

or 

F t = pVis{cotfi x - cot/?,) (4-2) 

From equations (4.1) and (4.2) it can be shown that 

C L = 2j/c(cot/f| - cot/i 2 )sin/) m - C D cot/i m (4.3) 

Since C D /C|. is approximately 0.05 for many blade sections, equation (4.3) 
may be reduced to a commonly used equation by suppressing the term 
containing C D : 

C L = 2s/c(cot/f] - cot/? 2 )sin/? ni (4.4) 

Since in machines blades move as well as being static, a blade moving across 
the flow field at a velocity U is shown in Fig. 4.4. Also shown are the 
velocity triangles at inlet and outlet, it being assumed that the fluid relative 
velocity directions are the same as the blade angles at inlet and outlet. Since 
the blade is moving and doing work, work done per second = F t U. Thus, 
since conventionally a work or specific energy coefficient is much used 
(defined as ip = gH/U 2 ), the work done equation may be rearranged to read 

V = C L - cosec /? m ( 1 + C D /C L cot/? m ) (4.5) 

2U s 

In many gas turbine textbooks and literature the axial direction is used as 
reference for angles, and so equation (4.5) may be rewritten in the form 

V = ^r;C L -secy? m (l + C D /C L tan/J m ) 

2U s 


(4.6) 
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Figure 4.4 Moving blade, the relative velocity triangles, and the forces acting. 


Since C D /C L may be neglected, equations (4.5) and (4.6) may be written 


V., c 

V = TT, Cl. - cosec fi„ 
2 U s 


(4.7) 


or 

H' = C L ~ sec fi m (4.8) 

1U s 

Thus, by using these equations and the ideal equations outlined in section 
1.2, it is possible to construct the velocity diagrams, obtain the fluid angles 
and determine the blade angles. By selecting a blade profile using the 
appropriate lift coefficient, the probable solidity, blade numbers and stagger 
angle may be computed. Before considering approaches to design in this 
way the data available for blade profiles must be discussed. 


4.3 Isolated aerofoil data 

The simplest section that may be used is the flat plate, for which Weinig 
(1935) produced a relation between angle of attack and lift: 

Ci, = Insina (4.9) 

This compares quite well for low angles of attack with experimental data 
presented by Wallis (1961; Fig. 4.5). Also shown is some information for the 
cambered circular arc plate. Simple plates of this sort are easy to manu- 
facture. and some low duty fans use them without marked loss of perfor- 
mance (Fig. 4.6) compared with an aerofoil of identical camber line shape 
but provided with a profile. In most cases strength considerations dictate 
some profiling and a number of sections have been used. Some of them were 
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Figure 4.5 Lift and drag characteristics of flat and cambered plates (tic = 0.02. 
R cC ~ 3 x 10 5 ) (adapted from Wallis. 1961). 


originally wing profiles, tested in wind tunnels at low speeds, which had a 
flat pressure side: for example, profile (c) in Fig. 4.7, the performance for 
which is presented in Fig. 4.8. Also shown in Fig. 4.7 are later profiles 
developed for compressors, which have been designed to give particular 
characteristics. Many families of profiles available are detailed in standard 
texts like Riegels (1961), and Abbot and Doenhoff (1959) (which concen- 
trates on the NACA profiles). Undergraduate texts deal with the hydro- 
dynamic theory which underlies the ’design' approach to profiles. Carter 
(1961) discusses this approach, and Fig. 4.9 from this source illustrates the 
effect of choosing a ‘concave’ velocity distribution (a) on the variation of 
drag with incidence, and contrasts this with a ‘convex’ distribution (b). As 
can be seen, there is a marked difference in behaviour. Carter commented 
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Figure 4.6 Eckert's results for an axial fan (R c = 3.25 x 1() 5 . same camber and 
rotational speed) showing comparison between the performance of cambered plate 
and profiled blades. 



Figure 4.7 Selection of profiles in use in axial machines: (a) axial compressor 
(British); (b) axial compressor (American); (c) fan profile; (d) curved plate; (e) 
transonic section; (f) supersonic section. 
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Figure 4.8 Characteristics of the Clark Y profile {tic = 10%, R cC = 8 x 10 6 ). 

that option (a) would give a better low speed performance, and option (b) 
will be better for high speed applications. In practical terms the choice of 
camber line geometry and the basic shape disposed round the camber line 
determine the blade performance. Some typical camber lines are illustrated 
in Fig. 4.10. It may be commented that if the point of maximum camber is 
towards the leading edge (LE) (well forward), the largest surface curvature 
and hence surface velocities will occur near the leading edge on the upper 
surface, producing a generally concave velocity distribution; a more rounded 
velocity distribution results if maximum camber is further towards the trail- 
ing edge (TE). It follows that profiles in the first category have a good low 
speed performance, a high stalling C L and a wide range of operation but are 
susceptible to cavitation in pumps, and those in the second have a lower 
range of operation and lower stalling C L but better cavitation behaviour. 

In practice the simplest camber line is a circular arc, but if the maximum 
camber is placed away from mid-chord a parabolic arc is needed. Carter 
suggested 40% of chord for maximum camber in fans and pumps, and 50% 
for compressors. The maximum thickness should be as small as practicable. 
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Figure 4.9 Effect of velocity distribution on the profile drag coefficient: (a) ‘concave’ 
velocity distributions, (b) ‘convex' velocity distributions (Carter (1961) courtesy of 
the Institution of Mechanical Engineers). 


and from data referred to by Carter it appears that thicknesses up to 12.5% 
of chord do not affect performance, and thicknesses above this dictated 
by strength considerations will result in a drastic fall-off in performance. 
Varying the position of maximum thickness has a similar effect to the 
movement of maximum camber, and for many sections the maximum thick- 
ness is located at 30% of chord from leading edge, as in the C4 section. The 
NACA 65 series all have the maximum thickness at 40% of chord. Table 4.2 
summarizes the variation in thicknesses for some commonly used profiles. It 
will be seen that the leading edge and trailing edge radii tend to be small, 
although the leading edge for low speed use tends to be fairly large to allow 
a tolerance of incidence changes typical of fans and pump usage, and in 
compressors the leading edges tend to be sharper. 


4.4 Cascade data 

When blades are in close proximity or ‘cascade’, the individual blade behav- 
iour is affected (Fig. 4.11) since the passages formed by adjacent blades 
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Figure 4.10 Common camber-line shapes. 


dominate the flow. A one-dimensional theoretical correction for this effect is 
to insert a correction factor in the Euler equation so that it reads, for a 
pump. 

gH — Ch(V u2 u 2 ~ V ul Mi) (4.10) 

This factor is shown in Fig. 4.12, and is a theoretical statement of the effect 
of blade spacing and blade angle. Weinig (1935) studied the two-dimensional 
problem by deriving relations for thin aerofoils which approximate to flat 
plates, and produced a lattice coefficient K used to correct the flat plate 
equation, equation (4.9): 

C L = 27r/fsina 

Figure 4.13 illustrates the effect on K of blade angle and spacing. 


(4.11) 


Cascade data 71 


Table 4.2 Thickness profiles of NA CA 6510, C4 and 
T6 sections 


% chord 

tl 2 

NACA 65 series 

tl 2 
C4 

t! 2 

T6 

0 

0 

0 

0 

0.5 

0.772 



0.75 

0.932 



1.25 

1.169 

1.65 

1.17 

2.5 

1.574 

2.27 

1.54 

5.0 

2.177 

3.08 

1.99 

7.5 

2.674 

3.62 

2.37 

10 

3.04 

4.02 

2.74 

15 

3.666 

4.55 

3.4 

20 

4.143 

4.83 

3.95 

25 

4.503 



30 

4.76 

5.0 

4.72 

35 

4.924 



40 

4.996 

4.89 

5.0 

45 

4.963 



50 

4.812 

4.57 

4.67 

55 

4.53 



60 

4.146 

4.05 

3.70 

65 

3.682 



70 

3.156 

3.37 

2.51 

75 

2.584 



80 

1.987 

2.54 

1.42 

85 

1.385 



90 

0.81 

1.60 

0.85 

95 

0.306 

1.06 

0.72 

100 

0 

0 

0 

LE radius 

0.687% of chord 

12% max. t 

12% max. r 

TE radius 

'sharp' 

6% max. t 

6% max. t 


Performance of cascades of blades was obtained using cascade tunnels 
sketched in Fig. 4.14 and by determining pressure and velocity changes at 
the mid-height of the centre blade and its associated passages. The tunnels 
had variable geometry to produce variations in incidence, and wall boundary 
layer control to ensure a two-dimensional flow field for the section under 
study. Figure 4.15 illustrates the typical performance of a cascade of NACA 
65 series blades. The design lift coefficient C L0 is related to camber for this 
foil shape in the way shown in the figure. 
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Figure 4.15 Data for an NACA 65 (27) 10 foil in cascade = 45°, sic = 1.0). 
(Note that 65 (27) 10 means a 65 series foil. C L() = 2.7. t.'c = 10%; C L o is design lift 
coefficient, related to camber angle 6 as sketched.) 


A pressure rise coefficient is often used: 


r _ P 2 - Pi 

p \pvl 

and a tangential force coefficient using Y as tangential force: 


( 4 . 15 ) 
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(quoted as y/ T by Zweifel (1945)). From these, 

C P = C f cot/? m - M (4.17) 

2 P v a 

Howell (1945) introduced a diffuser efficiency, defined as 


and it can be shown that 


Id 


Pi - Pi 

\piy\ - v\) 


>1d ~ 1 


2C D 

Cl cos 2/i m 


(4.18) 


(4.19) 


Howell also showed (assuming C L /C D is constant, which is approximately 
true for many profiles over a range of incidence) that /? m optimum is 45° and 
that 

'/□max =1-7^ (4.20) 

Cl. 

This suggests that >j D is maximum for a compressor row where /? m is 45°. 
Howell illustrates how t/ D varies with /?„, and the lift: drag ratio, and de- 
monstrates the small effect of variation in C L /C D for a conventional cascade. 
Figure 4.16 gives cascade performance plots presented non-dimensionally, 



Figure 4.16 Cascade data for a typical English cascade (after Howell (1945), cour- 
tesy of the Institution of Mechanical Engineers). 
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Figure 4.17 Variation of losses with flow coefficient for an axial flow compressor. 


and allows profile loss to be assessed over a range of flows and corresponding 
incidence angles. The other losses vary with flow rate, as Fig. 4.17 illustrates 
for a typical compressor stage. Cascade testing has revealed strong two- and 
three-dimensional flow patterns, the principal effects being illustrated in Fig. 
4.18. Boundary layer and tip clearance effects have a strong effect on the 
overall efficiency of a blade row and, owing to the secondary and wake flows 
affecting blades in succeeding rows, the performance of the whole machine. 
An example of this interaction is shown in Fig. 4.19, where time averaged 
axial velocity profiles are illustrated for four stages in a compressor. The 
effect is to cause a variation in work capacity of a blade along its length, and 
on the loading of successive rows. A much used correction in compressors is 
Howell’s work done factor Q, defined as used in 

C p AT = QuAV u (4.21) 

This factor is shown in Fig. 4.20. 

The flows in Fig. 4.18 show that fluid angles are likely to differ from 
actual blade angle, so that there is always a deviation angle to estimate 
which is affected by blade geometry. Howell proposed a formula for the 
deviation angle 3: 

S = m6V(s/c) (4.22) 

The well known plots of m against stagger angle for circular arc and para- 
bolic camber lines due to Carter (1948) are shown in Fig. 4.21. Carter’s 
discussion of blade profile behaviour has been referred to in section 4.2, and 
some comments relevant to cascade behaviour which extend the isolated foil 
now need to be made. Figure 4.9 shows a comparison of the behaviour of 
two types of velocity profile in cascade. Carter further discusses the effect of 
these profiles on loss (which is small) and shows that the effect of stagger on 
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Figure 4.18 Secondary flow effects in a compressor cascade. 


blade velocity distribution is small. Figure 4.22 is a reminder that as stagger 
increases, relative velocities must increase for the same flow rate and camber 
angle in compressor and turbine cascades. The basic difference between 
compressor (diverging) and turbine (converging) cascade passages is also 
illustrated. 

Figure 4.23 shows a zero reaction (impulse) blade cascade (b), illustrating 
the constant area passage and the blade shapes which follow. Also shown is 



78 Principles of axial flow machines 



0.6 0.8 1.0 1.2 1.4 

VIV 

v ' v a mean 


Figure 4.19 Velocity profile changes in a four-staae compressor (adapted from 
Howell, 1945). 



Figure 4.20 Work done factor proposed by Howell (after Horlock, 1958). 
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Figure 4.21 Deviation rules proposed by Howell and plotted by Carter (1948) (cour- 
tesy of the Institution of Mechanical Engineers). 
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Figure 4.22 Simplified presentation of the effect of stagger on cascade layout and the 
velocity triangles. 


a passage of a conventional reaction design (a), with a gradual contraction 
to the throat or minimum area. In Fig. 4.23(d) typical profile coefficients for 
impulse and reaction blades are illustrated, indicating the loss penalty 
imposed by the profile shapes, and the restriction in incidence range in- 
volved when impulse profiles are used. Also shown in Fig. 4.23(d) is the 
small variation in gas angle with incidence, demonstrating how the throat 
geometry has a great influence on outlet angle, a distinct difference from 
compressor cascades. A simple empirical rule much used by steam turbine 
designers is. following Fig. 4.24 

ci 2 = cos -1 o!s (4.23) 

where o is the cascade opening or throat. This is a good approximation to 
actual gas angles for high Mach numbers, but Ainley and Matheison (1957) 
comment that at low Mach numbers (0 < M n ^ 0.5) 

a 2 = [/(cos -1 ols)} + Asle (in degrees) (4.24) 

where /(cos -1 o/s) = —11.15 + 1 . 154 cos -1 ofs. and e is the mean radius of 
the flank leading up to the opening o. They also modified equation (4.23) by 
replacing o with throat area and s with annulus area just downstream from 
the trailing edges. 

The basic relations used in compressor and fan technology have been 
introduced in previous sections, and are in coherent groups, since in com- 
pressors both stator and rotor blades are treated as behaving in the same 
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Figure 4.23 (a) Typical shape of cascade or reaction’ (gas turbine section). Reaction 
blades may be standard section placed on camber line, like T6; (b) typical shape of 
'impulse’ cascade; (c) base profile of T6; (d) comparison of ‘impulse’ and ‘reaction’ 
cascade performance (after Horlock, 1966). 


way after allowance for rotor motion is made. In turbines, however, owing 
to the very different histories of steam and gas turbines development, a 
considerable range of relations exist, which need to be treated before con- 
sidering turbine design (section 8.3). 

A large range of equations has been used to express irreversibility in 
practical terms. A simple velocity ratio 


K = 


V 

V 

v lsen 


(4.25) 
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Figure 4.24 Empirical rule used by steam turbine designer. 


is much used, where F isen is the velocity that could be expected if expansion 
is isentropic to the prevailing back pressure. Kearton (1958) quotes three 
nozzle efficiencies, as already defined in section 1.7. Horlock (1966) quotes 
many different relations, among them Markov’s (1958) loss coefficients: 



(4.26) 


He also introduces stagnation pressure loss coefficients 


Pt 12 - P 2 

and shows that V N . for example, is given by 


Pi)2 rcl ffiQrd 
P(>3rel — Pa 


Y n ~ Cm l 


(4.27) 

(4.28) 


Clearly such indices are only valid if used when flow conditions give rough 
equality of M„. 

The velocity ratio o K is much used in steam turbine design for fixing the 
optimum efficiency point: 


u 

"^[2(^I()| — ^03ss)] 

Gas turbine engineers tend to use y/: 

7rr 

For R = 0, y/ = l/2crf s ; for R = 50%, y/ = 1 la\. 


(4.29) 


(4.30) 
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Figure 4.25 Forces on a fluid element in a rotating field of flow. 


For an isentropic turbine Vrr is unity, but rj TS < 1, and Horlock (1966) 
gives ‘design’ charts that may be used in situations where often // T s is a true 
reflection of turbine effectiveness. 


4.5 Radial equilibrium theories 

Theory and experimental information was based on two-dimensional test 
data, and corrections for three-dimensional flow' have been outlined in 
previous sections, but one effect that has not been discussed is the effect of 
the strong centrifugal forces exerted by blades. The centrifugal field distorts 
the flow velocity profiles considerably, and particles tend to move outwards 
rather than passing along cylindrical stream surfaces as classically assumed, 
particularly in low hub: tip ratio designs. An approach known as the radial 
equilibrium method, widely used for three-dimensional design calculations, 
assumes that flow is in radial equilibrium before and after a blade row, and 
radial adjustment takes place through the row. 

Consider a small element of fluid (Fig. 4.25) rotating in a pressure field at 
constant tangential velocity V u . If the angular velocity is o>, the centrifugal 
force is (pr dr dd)w 2 r on the element. Since V u = cor, the centrifugal force is 
pV^drdO on the element. 

The pressure force on the element, ignoring products of small quantities, 
is given by 

pressure force ~ r dp 66 

If the two forces are the only ones acting (viscous and other effects neg- 
lected), the particle will move at constant radius if: 

dp i'ii 

dr r 


(4.31) 
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or 


dp 

P 



(4.32) 


If the concept of stagnation pressure p () = p + \pV 2 is introduced, it can be 
shown that 


1 dpo = v d(V a ) | V u d (rV u ) 
p dr l ' dr r dr 

and if does not vary with radius, the usual radial equilibrium concept 
follows: 

d(fj d(rV u ) 

dr 2 dr 


or 


d(Vg) , 1 d(rV „) _ Q 
dr r dr 


(4.33) 


Although many approaches used in axial compressors are available as 
solutions to equation (4.33) the two of interest in pumps are the free and 
forced or ‘solid' vortex solutions. If the ideal flow situation is considered, 
Euler's equation will apply, which with equation (4.32) gives after inte- 
gration: 

V u r = constant (4.34) 

This is the free vortex law, and substitution in equation (4.33) results in 
V. d = constant. 

In a solid or forced vortex, 


— = constant (cl 
r 


(4.35) 


and equation (4.33) becomes 

d(V a ) 2 1 d(cV) 

dr r 2 dr 


a solution of which is 

V. t = V(constant — 2 c : r 2 ) (4.36) 

The free vortex solution is simple and, since V, is constant across the swept 
annulus, allows easy comparison of velocity triangles. The sections are laid 
out easily. As Horlock (1958) shows, the change in blade angles with radius 
is considerable, and Fig. 4.26 is a sketch of such a blade showing the ‘twist’ 
in the blade length. 
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Figure 4.26 Root and tip sections of a typical compressor rotor blade illustrating the 
degree of twist. 



Figure 4.27 Typical vortex and velocity distributions in an axial flow fan as flow rate 
changes (after Eck. 1973). 


It must be said that though this approach works at design flow, consider- 
able changes can occur in flow, particularly for long blades. Figure 4.27 
following material reported by Eck (1973) for fans, illustrates the flow 
changes likely. 


4.6 Actuator disc approach 

The assumption of radial equilibrium before and after each blade row is onlv 
valid if the blade aspect ratio is small; for long blades it ceases to apply. The 
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Figure 4.28 Actuator disc concept. 
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idea of the actuator disc where all the turning of the flow is assumed to 
occur in a very small axial distance appears to approximate to actual flow 
conditions in compressors and turbines. 

Consider Fig. 4.28: radial equilibrium applies far upstream and down- 
stream, but not necessarily in between, and changes in tangential velocity 
occur in the disc. Bragg and Hawthorne (1956) derived equations relating 
enthalpy and velocity changes, and Ruden (1944) extended the treatment by 
assuming that the radial displacement of streamlines is small, and that the 
V u term in equation (4.33) is small. The detailed theory may be followed 
through in the references, but an approximate solution due to Horlock 
(1958) will be discussed. He showed that the axial velocity changes follow a 
relation of the form 

AC, = rj --- (4.37) 


where h is the blade height and Z is the distance upstream or downstream 
from the disc. The axial velocity at the disc is 


C, 


±_K* 

2 


(4.38) 


Thus, following Horlock’s convention, V. d far upstream is given by 


C, = V al + 


V., 


a2 


Jti ) e ("Z) * 


(4.39) 


and downstream by 


C, = ^a2 - 


C,2 - V, 


al \g( — nZ\ h 


2 


(4.40) 
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Interference between rows is developed by Bragg and Hawthorne (1956) 
and Horlock (1958), and their contributions should be consulted for further 
information and the actuator disc technique. 


4.7 Stall and surge effects 


4.7.1 Introduction 

The phenomenon known as ‘stall’ affects axial flow pumps, fans and com- 
pressor stages when flow breaks away from the suction side of the blade 
aerofoil, an event associated with incidence. A typical characteristic is 
shown in Fig. 4.29. Surge occurs as an unstable and usually large fluctuation 
of the flow through a machine, which results in a complete breakdown in the 
flow; however, it is not directly related to stall. 

4.7.2 Stalling of fans and compressor stages 

The operation of single stage fans and compressor stages at flow coefficients 
low enough to cause stall has been studied by many experimenters, one of 
the most detailed studies being that of Smith and Fletcher (1954). Dunham 
(1965) suggests that the stalling phenomenon may be a rotating or pro- 
pagating stall, a band of reversed flow ‘anchored’ by spider supports, a so- 
called centrifuging effect (flow inwards at blade root, and outwards at blade 
tip) and a flow ‘in and out' at all radii. The last two effects probably include 
propagating stall. 

Rotating stall was observed by Whittle on the inducer vanes of a centrifu- 
gal compressor, and described by Emmons, Pearson and Grant (1955). 
When flow through the stage shown in Fig. 4.30 is reduced, the incidence 



Figure 4.29 Typical axial fan characteristic. 
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Figure 4.30 Rotating stall cell in an axial flow rotor row. 


angle increases and one or more blades may stall. Flow leaving the suction 
surfaces causes an effective partial blockage of a blade passage, which 
increases incidence on blade A; this blade will tend to stall, and blade C will 
come out of stall. The stall cell will thus tend to propagate round the row. 
To an observer on the blades the cell propagates in the opposite direction to 
blade motion at about half peripheral speed. Lakhwani and Marsh (1973) 
studied, among other factors, the influence of reducing blade numbers by 
half, and a compressor stage with 22 rotor blades and 33 stator blades. They 
concluded from the single rotor studies that stall cells did not move at 
constant speed while growing. However, once stall cells were established 
their angular velocity was a constant proportion of the rotor speed (/V), but 
it reduced as blade numbers halved from 0.56/V to 0.46/V. They found that 
adding stators downstream had little effect on stall onset. They also found 
that a theoretical prediction by Le Bot (1970) was fairly close for 22 blades, 
but not for 11 blades, indicating a need for more data. The difficulty of stall 
prediction, owing to the complicated cell patterns that occur, is discussed by 



Figure 4.31 Limits imposed by choking and stall on axial compressor performance. 
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Horlock (1958), where he considers the wide variations found in earlier 
experimental studies. 

4.7.3 Surge and stall in compressors 

The surge line on a compressor characteristic is the limit of stability, and is 
determined by testing. Figure 4.31 illustrates a fairly typical characteristic 
and also demonstrates the other limits posed on the performance envelope. 
Also shown is the common condition that for a wide range of lower mass 
flow the first stage at least is stalled, with other stages possibly so, and that 
at higher flows it is quite likely that the last stage will be stalled at higher 
pressure ratios. These phenomena are fully discussed by Horlock (1958) and 
Eckert and Schnell (1961) among classical texts. They argue that surge in 
the low speed region is associated with stalling of the first five stages, and 
in the high speed region with stall in the rear stages. For further infor- 
mation the classical paper by Huppert and Benser (1953) may be consulted. 
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Principles of radial and 
mixed flow machines 


5.1 Introduction 

This chapter is concerned with the basic principles governing the shapes of 
stream surfaces and of the blades for radial and mixed flow pumps, fans, 
compressors and turbines. In these machines flow either starts axial and 
changes direction as in pumps, compressors and bulb turbines, or begins by 
flowing radially inward and is discharged axially as in Francis and Kaplan 
turbines. In each case the rotor and stator systems produce the changes in 
direction while modifying the fluid energy level. 


5.2 One-dimensional approach 

The basic equations relating velocities, blade angles and work interchange 
between machine elements and the working fluid were introduced in section 
1.2, and these may be simply related to machine element design by assuming 
that the velocities are averaged across the passage at the relevant planes in 
the flow' path. These velocities are assumed steady, incidence and deviation 
are absent, and viscous and boundary effects are ignored. This approach is 
one-dimensional, and the ’design' problem is to relate the throughflow or 
meridional velocity to passage shape change. 

Consider Fig. 5.1, which illustrates a typical mixed flow rotor in sectional 
elevation, with the bounding surfaces and two stream surfaces. Figure 5.2 
illustrates the three-dimensional nature of stream surfaces, and has a 
velocity triangle sketched at point A. The velocity V n is simply related to the 
volumetric flow as described in section 1.2. 

If the same type of passage is considered (Fig. 5.3), the one-dimensional 
method of deciding the progress of the stream surfaces, once the passage 
walls are established, may be illustrated. A flow net has been sketched with 
two potential lines and three intermediate stream lines. Each stream line 
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represents a surface of revolution, as does each potential line. Noting that 
V n must be the same across each potential surface, it follows that 

2nRd = constant (5.1) 

for the flow field shown, along each potential surface. Applying this rule to 
each succeeding surface defines the progress of the stream surfaces, which 
are approximately true. 
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Figure 5.3 One-dimensional approach to channel How. 


5.3 Two-dimensional approach 


5.3.1 Passage shapes 

Strictlv, all flow motions will be three dimensional, and the one-dimensional 
approach is thus very approximate as a tool to describe them. In two- 
dimensional theory, flow is assumed to proceed along parallel planes, and 
much hydrodynamic theory may be used. This is based on the concepts of 
irrotational flow, w'here Bernoulli's theory that flow energy is constant is the 
connection between pressure and velocity distribution. 

A relevant flow- pattern is that obtained by superimposing a source and a 
free vortex to produce the spiral vortex of Fig. 5.4. This forms the basic 
model when considering flow' in simple radial flow pumps, as described in 
section 1.4. Since both the radial and tangential velocities vary inversely 
with radius (the free vortex law) it follows that /f is constant, and the 
resulting streamlines are logarithmic spirals. These relate fairly well to flow' 
in non-moving passages like volutes, but do not fit well with rotating blade 
passages, w'here the angle /T at the maximum radius is usually greater than 
/?i at inlet and l/ u2 > as Chapter 1 illustrates. Usually, blade profiles 
follow' a modified spiral, as discussed in section 5.3.2: for more theoretical 
approaches see section 5.4. 

Turning attention now to stream surface, we consider the contrast be- 
tween Figs 5.5 and 5.6. In the first a bend is shown with streamlines based 
on one-dimensional considerations. In this case pressure must increase 
toward the outer wall, and Bernoulli's theory cannot be satisfied since 
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Figure 5.4 Spiral vortex. 



Figure 5.5 One-dimensional stream surface solution to a bend. 


velocity must decrease toward the outer wall. Figure 5.6 shows how stream- 
lines must be adjusted to the requirements of flow net and Bernoulli prin- 
ciples. Acceptable stream surface shapes may be obtained using this 
principle, such as that sketched in Fig. 5.7 for a Kaplan machine. 

5.3.2 Impeller or rotating cascade 

A common technique for obtaining a radial cascade is to use conformal 
transformation of the linear or axial cascades discussed in section 5.4. a 
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Figure 5.6 Adjusted stream surfaces allowing for pressure and velocity distribution. 



Figure 5.7 Stream surface correction for a Kaplan turbine (dotted lines are the 
uncorrected surface shapes (after Pollard (1973), courtesy of the Institution of 
Mechanical Engineers). 


simple one being illustrated in Fig. 5.8. The equation of the vane after 
transforming a cascade of flat vanes is 


UJ = cot /? log,. 


r 2 


(5.2) 



94 Principles of radial and mixed flow machines 


b-d transforms 
to this curve 



Figure 5.8 Conformal transformation of a flat plate cascade to a circular cascade. 


and 


/ 

t 


Z 

2 7tsin jf 


1 

lOgc — 
r i 


(5.3) 


The ideal performance of the transformed cascade may be predicted from 
the relations already described in section 5.4. 

The following must be noted in comparison with the straight vane system. 
For centrifugal cascades, if fi x < /T, the angular momentum decreases with 
radius, and the curvature is less than the logarithmic spiral, the case for 
diffuser systems; and if [i x > /? 2 , angular momentum increases as radius 
decreases and vane curvature is less than the basic spiral (inlet guide vanes). 
In a turbine, if /?, < /T. angular momentum decreases with radius and 
curvature is greater than the basic spiral (this applies to the return passages 
in multistage pumps and compressors). 

The comments thus far relate to static cascades. In rotating cascades the 
irrotational concept breaks down, since flow departs considerably from the 
ideal shapes designed. Typical flow patterns observed are illustrated in Fig. 
5.9, and the velocity distribution across a passage for one flow is shown 
in Fig. 5.10. It will be noted that numerous studies, particularly in com- 
pressors, suggest that this distribution varies with flow rate and is not steady, 
as well as being related to blade passage guidance (blade number Z). 

A number of theoretical approaches have been made to the ‘correction' ot 
the ideal Euler-equation-based rotor outlet velocity triangle to give an actual 
change in V u rather than the ideal (Fig. 5.11). Busemann (1928) examined 
the cascade of logarithmic spiral vanes, and produced a modified Euler 
equation for a conventional backward curved bladed inpeller: 

gH = h {) u\ ± C H u 2 V R2 (col f} 2 ~ cot/?,) 


(5.4) 
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Effect of suction 
and pressure face 
boundary layers 

Figure 5.9 Three-dimensional flow effects in radial passages. 


Figure 5.12 is a plot of C H against sic. The factor C H corresponds exactly 
with the factor in Fig. 4.12 and is called the throughflow factor. The term 
including h„ is a correction for the shut valve energy transfer or displace- 
ment capacity and Busemann presented a series of plots for /t 0 against radius 
ratio for a range of values of fi 2 and demonstrated the effect of blade 
number. Weisner (1967) developed the Busemann approach during a study 
of the various slip theories bv using the approach of Stodola (1945) who 
proposed that the 'slip velocity’ v be given by the relation 


and. since 


1 7 R2 

tan/f 2 


U 2 — V u2 + 


+ V 


(5.6) 
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Shroud 



Figure 5.10 Typical velocity distribution in flow leaving a rotor passage. 



Figure 5.11 Ideal and actual outlet velocity triangles for a radial outflow impeller, 
or 


^u 2 + ^R 2 n Sin fi 2 

u 2 it 2 tan /?2 Z 


he proposed that 

_ V u2 V RZ tan/? 2 _ n sin fi 2 


iii 


Z 


(5.~ 
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Figure 5.12 Busemann's throughflow factor. 


He examined many correlations and compared data from over 60 pumps and 
compressors. From this study he proposed that the following relations 
applied within the limits specified: 


Vtt sin/F 


(5.8) 


up to a limiting value of 
R 2 


In" 1 


^8.16sin/Fj 


c limit 


and thereafter 

<7 


Vsin/F 



(5.9) 


(5.10) 


Figure 5.13 is a plot of the Weisner factor plotted against the blade number 
Z and the radius ratio for an outlet angle of 20 c . The full lines are the 
Busemann factor and the dotted lines are the Weisner factor. The limit of 
the radius ratio given by equation (5.9) is shown also in the figure. The 
dotted lines are those given by equations (5.8) and (5.10). 

Some empirical simple formulae which may be used as an approximation 
approach are as follows. 

• Pfleiderer (1961 ): 

8 ^ Euler 

Euler 


(5.11) 
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Z= 24 Z = 30 Z = x 



Figure 5.13 Plot of Weisner factor against blade number Z and the radius ratio for 
an angle of 20°. 


where gH | uk , r is the corrected value. 

C p = (0.55 -» 0.68 + 0.65 sin (U) 2 

• Stanitz (1952): 

//(or K) = = 1 - O.3150^sin^) (5.12) 

where <j> is defined in Fig. 5.14. For a radial machine where <f>l 2 = 90 c 
equation 4.6 reduces to 

H = 1 - 0.63/Z (5.13) 

• Stodola (1945) suggested an approximation to the slip velocity which has 
already been quoted as equation (5.5). 



Three-dimensional problem 99 




Figure 5.15 Sketch of Whitfield's jet model. 


There are discrepancies between these approaches as Stahler (1965) has 
demonstrated. His article proposes an alternative relation that does not 
correlate with the equations quoted above. 

A more recent approach based on turbocharger design is that by Whitfield 
(1974) based upon radial bladed compressor impellers, draws attention to 
the lack of guidance of the emerging flow, proposes a jet model and suggests 
a modification to the Stanitz approach: 

K, - — = 1 - K,[l - K*j 

U 2 

where K* = 1 — 0.630/2 (6 defined as in Fig. 5.15. A relation of with 

blade number is given in Fig. 5.16. Whitfield's work supported the cor- 
relation. and further refinement is detailed in Whitfield, Doyle and Firth 
(1993). 


5.4 Three-dimensional problem 

Figure 5.9 indicates pronounced three-dimensional flow effects, some of 
which are a product of boundary layer behaviour on the passage walls. If the 
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Figure 5.16 Plot of the velocity distribution factor proposed by Whitfield (1974) 
(courtesy of the Institution of Mechanical Engineers). 



Figure 5.17 Approach to describing a three-dimensional flow in a curved passage. 
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mixed flow model used in the previous section is used again, and flow is 
considered to be a combination of rotational and meridional effects but 
without vanes, the rotation follows the law of constant angular momentum, 
and meridional flow is described by the stream surfaces as before. However, 
the meridional velocity may not be constant but may be affected by the flow 
curvature. Assuming constant energy, and considering the passage in Fig. 
5.17, where both boundary surfaces are concentric about centre C, the 
velocity distribution along a potential line such as AA will follow the free 
vortex law: 

V„X = constant (5.14) 

The constant depends on the line of constant potential. If continuity is 
applied, the total flow rate is given by 

Q = X V„2nrd (5.15) 

or if a defined level at a distance y is V ny , 

Q = T ‘{^x l ) 27ird (516) 

It will be noted that increments of flow between adjacent streamlines must 
be the same: V n 2 nrd = constant. Thus the streamlines cannot be concentric. 
This, although an improvement on the earlier methods described, must also 
be modified to allow for the boundary layer with an injection of an assumed 
velocity profile close to the walls. Practice appears to be that the real fluid 
effects are ignored and the approximate approaches outlined are used. 


5.5 Discussion of theoretical approaches to analysis and design 

The idealized approach to machines has been to assume inviscid steady flow 
and to apply the hydrodynamic equations. Actual flow is three dimensional, 
unsteady and viscous, making rigorous mathematical solution extremely 
difficult. Studies of centrifugal and mixed flow machine theory began by 
assuming steady inviscid flow on asymmetric stream surfaces before and 
after the rotor. Later contributions sought to modify the flow patterns round 
critical surfaces, like the leading and trailing edges, to allow for boundary 
layer separation, and to change stream surface curvatures to incorporate 
corrections for three-dimensional flow effects. 

One large school of solutions has followed the principles discussed by Wu 
(1952) who proposed the use of two surfaces - the SI plane from blade to 
blade across a channel, and the S2 plane from hub to shroud (that is, along 
the length or height of the blade). The solution then proceeds on the basis 
of a two-dimensional stream surface shape through the machine following 
the meridional approach already outlined, with a pseudo three-dimensional 
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correction to the surfaces on the basis of a loss theory. The blade sections 
are then designed along each surface in turn from hub to tip section. An 
early English example of this type of solution is that due to Marsh (1966), 
which was based upon inviscid flow; this was subsequently improved upon 
by Perkins (1970). A later technique introduced the effect of viscosity by 
incorporating a loss model in the form of a polytropic efficiency based on 
experimental relations. This technique was used by Pollard (1973), for 
example, to both analyse a number of machines and to perform design 
studies. He considered a typical pump and a Francis turbine, establishing 
first the stream surfaces, then performing blade to blade solutions on each 
surface in turn, utilizing a conformal transformation due to Wilkinson (1968) 
that produced a planar cascade. This cascade may then be used to study 
velocity distributions with inlet and outlet angles and pressure distributions 
that, when considered satisfactory, are then transformed back to the stream 
surface. This process, repeated for each stream surface in turn, allows the 
blade profiles to be built up from hub to tip. Pollard demonstrated how the 
predicted pressure distributions matched experimental data for a pump, and 
showed quite good agreement apart from the leading edge and trailing 
edge areas. These techniques are being continuously improved, and many 
contributions correlating experimental data to provide corrections to the 
computer-based analytical solutions have been published. Some are referred 
to in context in later chapters. 

An example of the combination of an analytical solution approach with 
experimentally obtained data is that due to Stirling and Wilson (1983). The 
design of a mixed flow pump, illustrated in Fig. 5.18, is described in some 
detail. They discuss their earlier work using a mainframe computer at the 


Stator 



NEL, and demonstrate how a desktop computer system may be used to 
produce the same design. Their method first determines the casing profiles 
and the most suitable positions of the blade leading and trailing edges B), 
B 2 . B 3 and B 4 ; the spacing of blade rows are about a half-blade chord in 
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practice, and the shapes are optimized to give good flow characteristics. The 
method then proceeds to predict the velocity distribution without blades, 
and moves on to select blade energy distribution and to look at the blade 
solution for the rotor, including losses based on experimental and theor- 
etical correlations taken from many sources. The rotor blades are then 
designed, and the stator blades are designed in similar routines, and the 
complete machine drawings may then be produced for manufacture. They 
demonstrate the correlation between prediction and test for pump designs 
having dimensionless specific speeds of 1.8, 2.5 and 3.2. The validity of the 
approach, compared with the experimental and empirical approaches, was 
demonstrated when pumps were designed by several engineers and quite 
close agreement found between the techniques deployed. These are de- 
scribed in the papers by Bunjes and Op de Woerd (1982). Richardson 
(1982), Thorne (1982) Stirling (1982) and Chiappe (1982). 

A survev undertaken by Hughes. Salisbury and Turton (1988) reviewed 
commercial computer-based methods of analysing flow through machines 
pump and compressor passages, and also approaches to design. There are 
now a number of manufacturer-produced software systems that interface 
with computer-based manufacturing systems, and also give user friendly 
method of pump selection. These are reviewed at intervals in such journals 
as World Pumps. 
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Centrifugal machines 


6.1 Introduction 

Centrifugal machines of the pressure increasing type are covered here: these 
include pumps, fans, and compressors or blowers. In all these designs the 
fluid enters axially and is discharged by the rotor into a static collector 
system and thence into a discharge pipe. Pumps and fans have a common 
feature in that the fluid is considered incompressible, and in the compressors 
and blowers there is usually a considerable density change. 

Many machines are single stage designs as sketched in Figs 6. 1-6.4, and 
may either have an axial inlet or a ducted inlet as shown in the double entry 
pump shown in Fig. 6.2. The single stage pump is the workhorse of the 
process and petrochemical industries, producing specific energy rises from 
about 50Jkg _1 to over 1200Jkg~' (head rises from 5 to I25m of liquid) and 
delivering flow rates from 6.3 to 400nr ! h -1 . Single stage compressors in 
turbochargers, for example, give pressure ratios up to at least 4: 1 but more 
conventionally will give about 2.5:1; fans used for ventilation duties give 
pressure rises equivalent to a maximum of around 500mm water gauge. 
Multistage machines like that shown in Fig. 6.5 are used for boiler feed 
and similar duties, and may be called on to deliver an energy rise of 
300000 J kg -1 (or head rises up to about 3000 m of liquid). 

Since the flow paths are the same for the different machines, the main 
considerations involved in the suction area, the impeller and the stator 
delivery system, are discussed in that order, with some attention being paid 
to the design choices open to the engineer. 

The chapter concludes with a discussion of radial and axial thrust loads, 
and a brief discussion of the additional problems posed by the mixed flow 
machines. 


6.2 Inlet or intake systems 

The simplest inlet system found in pumps, fans and compressors is the 
straight pipe coaxial with the impeller centreline, as seen in Fig. 6.6. but 



Inlet or intake systems 105 



Figure 6.1 Simple centrifugal pump. 



Figure 6.2 Double suction centrifugal pump. 




Figure 6.3 Centrifugal fan. 


Radial flow fan 
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Centrifugal compressor 
Figure 6.4 Centrifugal compressor. 


space and suction system layout frequently requires a ducted inlet as shown 
in Fig. 6.7. These can cause undesirable inlet velocity profiles and strong 
three-dimensional flow patterns in the eye of the impeller. This is particu- 
larly true in double suction machines where the shaft passes through the 
flow on both sides of the impeller, and a good treatment of this problem is 
found in the paper by Matthias (1966). Flowever good are the design flow 
patterns, there will be problems at part flow with vortex formation in the 
straight suction pipe which penetrates far into the suction from the suction 
flange in the manner shown in Fig. 6.6. 


6.3 Impeller 


6.3.1 Eye or inducer section 

As the fluid approaches the impeller it slowly takes up its rotation so that 
the absolute velocity begins to increase with a resulting drop in the static 
pressure (Fig. 6.8). This process continues till the fluid reaches the leading 
edges of the rotor vanes, and then has to pass round the nose and into the 
impeller passages. The pressure may tend to drop, depending on the in- 
cidence angle, and may in pumps reach the local liquid vapour pressure with 
resulting cavitation problems. It is conventional to make allowance for 
incidence angles up to 6° in design in order to allow for these effects. 

In compressors, which tend to have radial blade sections at the outlet, the 
inducer or inlet section is given a large amount of twist, as indicated in Fig. 
6.9, to accommodate the relative velocity direction imposed by the inlet 
motion of the gas. The incidence angle rule just noted for pumps is used but 
great attention often has to be given to the shape of the nose. This is 
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Figure 6.5 Multistage centrifugal pump. 


illustrated in Fig. 6.9. which shows the effect of three nose profiles and their 
effect on the local Mach number, shape 1 being preferred to allow smooth 
local velocity changes. 

The axial velocity is also an important parameter, and in pumps is kept as 
low as possible to avoid cavitation problems. The size of this velocity is 
related to the eye or suction diameter, and empirical rules will be found in 
Stepannof (1957b) and Turton (1994) among other texts. The relations 
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Figure 6.6 Simple suction pipe. 



Figure 6.7 Ducted entry suction. 


between axial velocity and suction diameter were discussed in section 3.4, 
and the risks of oversizing were outlined. In compressors the need to control 
the velocity at the maximum diameter of the suction zone to avoid Mach 
number problems is paramount. Figure 6.10 illustrates the probable design 
relations that can be involved and used as design guidelines. 

To illustrate the basic approach, consider the following compressor 
example. A compressor is to deliver air at the rate of 9 kg s ' drawn from a 
plenum chamber where the conditions are 10 5 N m" 2 and 295 K. The plenum 
conditions are stagnation. Thus, if it is assumed as the ideal case that there 
is zero inlet whirl for the impeller and an acceptable axial velocity com- 
ponent of 75 ms -1 , the static conditions in the inducer, and the density and 
other parameters, may be found. Assuming the stagnation conditions 
remain the same along the inlet pipe (approximately true if it is short), 

295 = T t + 75 2 /2C p 
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Figure 6.8 Pressure changes on a stream surface in the suction of a centrifugal 
machine 



Figure 6.9 Inducer of a centrifugal compressor. 
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Figure 6.10 Design plot for a compressor/inducer. 


thus 


T, = 292.2 K 

Following the gas laws. 

Pm'Pi = 

Therefore 

Pt = 0.997 x 10 5 N m -2 
/; = 1 .189 kg m -3 

Assuming a hub to tip diameter ratio of 0.35, 

9 7t , 2 

1.189 x 75 = 4 (D ‘ ,p “ ° r,ub) 

D tip = 0.385 m 
D hub = 0.135m 

The acoustic velocity at the inducer section is given by 
a = \/(yRT) = 342.65 ms -1 

The velocity triangle is drawn for the maximum diameter of the inducer, 
using the assumption of zero whirl and a driving speed of 15 000 rev min -1 , 
i.e. a tip peripheral velocity of 302.38 ms -1 . This gives W t = 311.54m s -1 . a 
Mach number based on W\ of 0.91, and /?, = 13.93' (with incidence prob- 
ably 16°). Similarly, for the hub diameter, W\ = 129.87ms -1 and /f, = 
35.27°. 
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6.3.2 Impeller design 


The inducer (or eve) of the impeller has been discussed, and attention may 
now be devoted to the choice of blade or passage design. 

As in all turbomachines there is a need to balance control with friction 
loss due to guiding surface. It is possible using recent computer techniques 
based on loss data to obtain blade numbers, which can then be fed into the 
suites of programs for passage shapes. Empirical formulae for pumps and 
compressors based on ‘good practice’ will be quoted here. 

For pumps (quoted by Stepannof (1957b), for example) blade number is 

Z = (i 2 l 3 (6.1) 

Another, quoted by Pfleiderer (1961) and others, is 

Di + D * 

Z = 6.5 — = ~sin/f m (6.2) 

D 2 LJ [ 

In these equations. D ] and D 2 are inlet and outlet diameters of the impeller, 
respectively, jl 2 is the impeller outlet angle, and /f m is the arithmetic mean of 
the inlet and outlet angles of the impeller, quoted with respect to the 
tangential direction. 

Equally simple rules exist for compressors, such as 


Z = fUJ 3 

Another, quoted by Eckert and Schnell (1961). is 

7 = 2tt sin /f nl 

(0.35 to 0.45) log c D 2 ID , 


(6.3) 


(6.4) 


Pump design conventions indicate backward curved blades with outlet 
angles f} 2 in the range 15° to 35° or 40°, the lower values being found 
typically in high specific speed designs. Compressors have tended for stress 
reasons to have outlet angles of 90°. but some modern designs have been 
provided with backward curved vanes; McCutcheon (1978) quoted an angle 
of 33° for a 4.4 : 1 pressure ratio machine with a rotational speed of 
70000 rev min’ 1 . 

Varley (1961) studied the effects of blade number on pump performance, 
with the effects of blade angle and surface roughness as additional factors. 
The machine used was a double entry design with an impeller diameter of 
244mm. designed to give a duty of 16m at 6.44 Is -1 when rotating at 
1400 rev min -1 . Figure 6.11 indicates that high numbers of blades resulted in 
the head characteristic falling to zero flow, a condition known as unstable. 
The best shaped curves apply to five and six blade designs, compared with 
the seven suggested by equation (6. 1 ) on the basis of an angle of 27°. Figure 
6.12 demonstrates that higher outlet angles gave unstable curve shapes, that 
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Figure 6.11 Pump characteristics with varying blade number (after Varlev (1961) 
courtesy of the Institution of Mechanical Engineers). 



Figure 6.12 Effect of outlet angle on pump characteristics (after Varley (1961) 
courtesy of the Institution of Mechanical Engineers). 
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Figure 6.13 Angle of blade overlap. 


27° gave a characteristic which rose steadily to shut valve, regarded as a very 
desirable feature, and that this angle gave the best efficiency. Varley studied 
the slip correlations outlined in section 5.3.2. and found that the Stodola 
(1945) and Busemann (1928) formulae compared well with his data for low 
blade numbers and angles, but deviated considerably as both angle and 
blade number rose, particularly that due to Stodola. 

Another design consideration related to blade number is the angle of 
overlap, illustrated in Fig. 6.13 for a pump design. Texts that deal with 
pump design suggest an angle of 45° as a good compromise for pumps of 
conventional designs. Centrifugal fans tend to have more but shorter blades: 
Eck. in his classic text (1973). demonstrated this for the usual circular arc 
blades, and showed the discontinuity in characteristic curves due to local 
stall effects. He gives many working formulae for blade and pitch based on 
good practice. Some data from work by Myles (1969: Fig. 6.14) indicates 
that for fans, as for pumps, the lower angles give rising characteristics. This 
paper also gives useful data about losses and their codification in a manner 
similar to the compressor diffusion factor. 

Optimum proportions for compressors follow the same principles already 
outlined, but these have had to be modified to account for compressibilitv 
effects, and to some extent for stress considerations in high speed machines. 
As has already been stated, the blades tend to be radial at outlet, so that 
equation (6.3) indicates 30 blades. Although this will give good control at 
outlet and reduce slip, there will be too much metal in the inducer area, so 
that the usual arrangement is to make every other blade a half blade as 
shown in Fig. 6.15. McCutcheon (1978) describing the backward curved 
impeller, showed half blades; he reported that empirical rules for blade 
number were checked by finding blade loadings as a part of the design 
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Figure 6.14 Effect of outlet angle on fan characteristics (after Myles (1969) courtesy 
of the Institution of Mechanical Engineers). 



Figure 6.15 Provision of half vanes in a centrifugal compressor. 


process using blade to blade methods, and he found quite close agreement 
with the empirical approach. 

Another important consideration is the shape of the blade passages. The 
side elevations found in good pump practice are shown in Table 6.1, and the 
placing of the inlet edge in elevation is seen to vary with specific speed. In 
blowers the leading edges tend to be parallel to the axis of rotation; the flow' 
is high but pressure rise relatively small. In the majority of compressor 
designs the tendency is to provide inlet edges close to the plane at right 
angles to the axis of rotation, allowing for as long a passage length as 


Table 6.1 Pump side elevations 
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Area of sections before adjustment 



(b) 

Figure 6.16 (a) Method of establishing passage area changes (at point P,. area A\ = 
b\a{lnr\). (b) Cross-sectional area of impeller passage plotted against passage length. 


possible to give good flow control. Blade passages are designed to give a 
controlled diffusion from inlet to outlet, so that the area changes required to 
give this must be determined by the profiles of the shroud and the back- 
plate, and by the blade shapes too. The texts cited give ground rules for 
blade profile: in pumps and compressors the blades will follow an ap- 
proximate Archimedian spiral, either using point to point programmes or 
combinations of radii. When the blade profiles are set and the elevation 
profiles determined the passage area changes must be checked as illustrated 
in Fig. 6.16. The principles outlined in Chapter 5 may be used, but the area 
changes must give smooth changes from inlet to outlet; if this is not achieved 
by the shapes designed, modifications need to be made and the process 
followed through again. 

In compressors the pump approach just outlined may be followed but 
compressible considerations demand care. Figure 6.17, for example, shows a 
compressor design obtained by using meridional and blade-to-blade sol- 
utions that, by successive passes of calculation, achieve a smoother velocity 
variation. Detailed discussion is not possible here, but the contributions 
cited in Chapter 5 may be consulted as well as those given in this chapter. 
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Figure 6.17 Effect of wall shapes on local relative velocities in a compressor. 

The main features governing the design of good performance pumps and 
compressors have been outlined, but before leaving a discussion of impellers 
it is of interest to consider again the paper by Varley i 1961). He studied the 
effects of surface roughness on performance. Figure 6.18 indicates that 
roughness cannot be ignored as unimportant, so that cost reductions which 
lead to deliberately poor surface production give rise to performance penal- 
ties. This is true for pumps, fans and compressors, although the latter, being 
‘high technology' machines, normally have high class finishes on all flow 
surfaces. In fans and pumps, of course, the selling price is also a consider- 
ation that may modify material and surface finish choices. Other factors 
affecting performance are clearances between rotating and static surfaces, 
and their influence on leakage flow. This was briefly remarked upon when 
discussing the approach to performance prediction by Nixon, in section 
2.4.1. For efficient operation clearances should be a minimum, but toler- 
ances in manufacture, and such effects as end float in bearings and thermal 
gradients give rise to differential rates of expansion in many machines, 
which result in large cold clearances in compressor^ and normally w-ide 
clearances in pumps and fans. That care is needed in settling dimensions is 
indicated by a well known manufacturer of turbochargers, who reported a 
drop of 2% in performance as the side clearance between impeller and 
casing (expressed as a proportion of impeller diameter increased from 0.036 
to 0.143. Pampreen (1973) quoted a similar drop in a compressor passing 
0.904 kg of air per second when this clearance ratio changed from near zero 
to 0.06. 
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Figure 6.18 Pump characteristics with various impeller components roughened (after 
Varlev (1961) courtesy of the Institution of Mechanical Engineers). 


Commercial considerations result, particularly for pumps, in the practice 
of providing a range of impellers in one casing, and also in reducing the 
diameter of an impeller (‘turning down') up to values approaching 15% of 
the maximum size. Head and flow may be predicted using the so-called 
scaling laws: 


ft = A '?P-? 
P 2 K2D2 


(6.5) 


Qa = Ni D\ 

q 2 n 2 d 2 

8 = N\ £l 

gH 2 N\Dl 


( 6 . 6 ) 


(6.7) 


These laws allow the likely performance to be predicted from the full size 
characteristics, but do not with accuracy allow prediction if the impeller is 
cut down to less than 75% of the full size. Another well known ‘adjustment’ 
technique is the process of filing the vane tips, which can make small but 
significant adjustments in performance, as Fig. 6.19 illustrates. 

For further illustration the example in section 6.3.1 will be continued. 
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Figure 6.19 Effect of centrifugal pump blade tip shape on performance (after 
Worster and Copley {1961 ) courtesy of the Institution of Mechanical Engineers). 


Since the pressure ratio required is 3.2:1. and the efficiencs suggested is 
80%. 

3.2 = (1 + 0.8A 777;,,)- 1 - 5 

Therefore. 


Assuming zero inlet whirl 


AT = 145.37 
T o: = 440.37 K 


V a2 u 2 = C p AT 

Assuming as a first trial that D z lD , is 2. 

90 + 13.93 


flu 


_ <po 


Then, from equation (6.4). 


27tsin52 c 
0.45 log c 2 


16 


Overall efficiency (%) 


1 20 Centrifugal machines 


Using the Stanitz equation, equation (5.12), 

/ 1 = — — = 1 — 0.315(^sin9o) 
u 2 V 16 / 


which gives 
and therefore 
Thus 


H = 0.876 

0.876//1 = 1.005 x 10 3 x 145.37 

u 2 = 408.4 ms -1 
V u2 = 357.76 ms -1 


and, for 15 000 rev min *, 


D 2 = 0.52 m 


The passage height at this diameter is now needed. The density is required 
as well as the outlet radial velocity, and to establish this some estimate of 
efficiency or loss is required. For example the outlet stagnation conditions at 
the discharge flange are known from the overall pressure ratio already 
calculated. As a first trial calculation it may be assumed that half of the loss 
of 20% implied by the efficiency of 80% occurs in the casing. Therefore a 
loss of perhaps 10% may be assigned to the impeller passages, giving an 
efficiency for the impeller of 90% as a guess. Thus, assuming that the 
stagnation conditions remain the same from impeller outlet to the discharge 
flange, a calculation may be made as follows: 

pmlpm = (1 + 0.9 x 145.37/295) ' 5 = 3.75 

Thus 

P( 12 = 3.75 x 10 5 Nm -2 
r (l2 = 440.37 = T 2 + V\/2C P 

Since in the limit V 2 must not really reach Mach 1. the maximum velocity is 
fixed by assuming that V 2 equals the acoustic velocity; V 2 = 384 ms -1 . and 
hence T 2 = 367 K. Using the gas laws 

p 2 — 1.98 x 10 5 Nm -: 

p 2 - 1.88 kg m -3 

Q 2 = 4.79 m 3 s -1 

If now a diffusion is assumed such that V R2 = 0.95V r Rl , 

U R2 = 71.25 ms -1 
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Using Q 2 shows that an outlet area of 0.067 m 2 is required. Allowing for the 
blockage effect of 16 vanes each 2 mm thick and, if the passage height is h, 

0.067 = n x 0.52 x ft - 16 x 0.002 x h 


which gives 


h = 42 mm 


6.4 Outlet systems 

Fluid leaves the impeller at an absolute velocity considerably higher than 
that in the discharge pipe. The outlet system, therefore, as well as collecting 
fluid from the impeller without affecting its performance, must also reduce 
the flow velocity efficiently and so improve the outlet pressure from the level 
achieved by the impeller working on its own. Figure 6.20 illustrates the basic 
layouts used in fans, pumps and compressors. All the alternatives have an 
annular space outside the impeller before the volute or diffuser ring. The 
vaneless diffuser continues the annular space and will be discussed first; the 
volute will then be described, and finally the vaned diffuser. 




(b) 

Figure 6.20 Alternative diffuser systems used in pumps and compressors: (a) simple 
volute casing with (i) tangential discharge (ii) radial discharge: (b) diffuser system 
with (i) vaned diffuser (ii) vaneless diffuser. 
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6.4.1 Vaneless diffuser 


A two-dimensional model of the flow in a vaneless diffuser is represented in 
Fig. 6.21. The streamlines follow a spiral vortex path that, shown in the 
figure, traces a streamline from the impeller outlet at point 2 to the casing at 
point 3. the space between the points being the diffuser. Application of the 
principle of constant angular momentum and the continuity equation be- 
tween points 2 and 3 gives (ignoring losses) an expression for the pressure 
rise as follows: 


^FUiffuser P?> P2 




( 6 . 8 ) 


This is the ideal pressure rise, and if flow losses between points 2 and 3 are 
accounted for by using the diffuser efficiency, then 

V diffuser = ( 6 . 9 ) 

Ferguson (1963, 1969) discussed flow in vaneless diffusers and the ap- 
proaches to the estimation of losses. He demonstrated (Fig. 6.22) uncertain- 
ties in flow' regime, and strong three-dimensional flow as flow' rate changed. 
He discussed the relevance of the friction loss estimation approaches used in 
normal diffusers, and concluded, as did Johnston and Dean (1966) that all 
the available theoretical approaches have to be modified in the light of the 
flow instabilities found in experimental studies. Of great interest are studies 
published by Sherstvuk and his colleagues (1966, 1969) on the effect of 
profiling the side walls on the diffuser efficiency. Figure 6.23, taken from 
this work, illustrates the effect of quite limited changes in principal dimen- 
sions. Clearly the control of boundary layer behaviour pays dividends. More 
recently, Whitfield and Sutton (1989) have showm how geometry affects 
surge margins and flow range in compressors. 
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Flow increasing 


Figure 6.22 Flow patterns in a vaneless diffuser as pump flow changes (Ferguson. 
1969). 




Figure 6.23 Effect of geometry on vaneless diffuser performance (after Sherstyuk 
and Kosnrin. 1966: Sherstyuk and Sekolov, 1969). 


6.4.2 Volute or spiral casing 

The volute surrounds the impeller, and is a spiral casing whose cross- 
sectional area increases from a minimum at the cutwater to the throat, 
which is the beginning of the diffuser to the outlet pipe, as illustrated in 
Figs 6.24 and 6.25. The cutwater represents the nearest part of the casing 
to the impeller, and is aligned into the general direction of the flow leaving 
the impeller. The cross-sectional area increases to take the volume leaving the 
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Figure 6.24 Volute diffuser. 



Figure 6.25 Volute and a typical cross-section. 


impeller: for example, in Fig. 6.25, if section AA is at 6 = 90° from the 
cutwater, 25% of the pump delivery flow will be passing through the section. 

If ideal flow is considered, the law of constant angular momentum will 
apply, and the flow will leave the impeller on a logarithmic spiral path. In 
many casing designs the cross-sections are circular sections, and the outer 
wall follows an Archimedian spiral. It is possible, using integration across 
each section of small elements indicated in Fig. 6.25. to schedule the area 
changes and get the correct areas for all designated sections like AA. The 
equations and solutions are covered in textbooks like those of Stepannof 
(1957b) and Pfleiderer (1961) though usually the effect of friction is 
neglected as it has been found to affect the areas by very small amounts. A 
much simpler approach is to size each cross-sectional area by assuming that 
the throat velocity is constant round the impeller circumference. Stepannof 
(1957b) also demonstrates this technique, based upon a velocity given by the 
relation 


V = K,V(2gH) 


( 6 . 10 ) 
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Figure 6.26 Effect of the choice of velocity distribution on pressure head variation 
around the impeller: (a) constant angular momentum; (b) constant velocity (Rutschi. 
1961). 


In this equation Ky is an empirical factor plotted against specific speed, and 
gH is the pump nominal design specific energy rise. Rutschi (1961) pre- 
sented a comparison of the angular momentum and constant velocity ap- 
proaches which is shown in Fig. 6.26, and demonstrated little difference in 
pressure head at successive sections round the casing. 

Worster (1963) more recently studied the casing and its effect on pump 
performance and proposed that the most important parameter was the 
volute throat area. He argued, as did Anderson and others referred to in his 
paper, that the best efficiency or design match between the impeller and the 
casing occurs when the impeller and volute characteristics cross, as indicated 
in Fig. 6.27. He based his analysis on the simple casing cross-section shown 
in Fig. 6.27 and maintained that, although the outlet angle / h and impeller 
diameter Z? 2 do determine performance, the correct matching of impeller 
and casing is a decisive factor in the design process. Figure 6.28 illustrates 
his point, and Fig. 6.29 is his design plot that allows the selection of the 
correct volute area once the impeller is designed. This approach has been 
used with some success, and Thorne (1979) has recently underlined its 
utility. 

Although the technique just discussed has been applied to pumps, it is 
possible to argue that it could be applied to fan and blow'er design since the 
density changes are low'. However, as Eck (1973) show's, the principles of 
constant angular momentum or constant velocity are both used with empiri- 
cal factors such as equation (6.10). 

Considering pump design still, Worster (1963) used the loose or large side 
clearance design approach, as against the more conventional tight or close 



1 26 Centrifugal machines 



Figure 6.27 Illustration of the approach to impeller and casing matching proposed by 
Worster (1963). 



Figure 6.28 Illustration of volute area change on pump performance (after Worster 
(1963), courtesy of the Institution of Mechanical Engineers). 
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k s 

0.8 0.4 0.2 



y = Outlet area between vanes of impeller/throat 
area of casing 

Figure 6.29 'Design' plot proposed b\ Worster (1963: modified from that quoted by 
Worster, the symbols are defined in Fig. 6.27). 


side clearance designs. If the loose design is used the width x in Fig. 6.25 
will be the casing width needed to give adequate side clearance for the 
impeller. If the tight design is used the dimension ,v should be fixed, as for 
the vaneless diffuser, to be a little larger than the impeller channel height. 

The diffuser from throat to outlet may be a cone or. if the centreline 
outlet shown in Fig. 6.20 is used, as demanded by the current ISO standards 
for small and medium size process pumps, may be curved. Conventional 
pipe diffuser data are not relevant to this type of diffusion since the flow in 
the casing is complicated, and the presence of distinct jets of fluid issuing 
from the impeller passages in succession will be causing disturbed three- 
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Figure 6.30 Cutwater nose profile. 


dimensional flows into the throat of the casing. The shape of the cutwater 
nose is also important, as indicated by information published by Turton 
(1966). It is also well known that the clearance between the impeller and the 
cutwater nose is important, as too tight a clearance gives rise to noise and 
vibration. A working rule is to arrange that the minimum clearance should 
be of the order of 5-10% of impeller diameter. Depending on casting, the 
cutwater nose may be blunt or sharp (Fig. 6.30). but must be aligned along 
the streamline as shown in the figure. A blunt nose allows for tolerance of 
flow misdirection and thus an accommodation for flow changes from design 
point. 


6.4.3 Vanecl diffuser systems 

The diffusion system is now a vaneless space whose radial depth follows 
the same proportions as already given for the volute cutwater clearance. 
The fluid follows a spiral path up to the diffuser throat, is diffused over the 
length be (Fig. 6.31) and discharged into the spiral casing. Ideally the 
diffuser vane walls should follow an Archimedian spiral but. since diffusers 
are often machined, the surfaces abc and ef are convenient radius ap- 
proximations. Each diffuser passage has an active length L. which a working 
rule suggests can be between 3/ and 4 1 to balance effective diffuser control 
and surface friction loss. Clearly the sum of throat areas will follow the same 
basic rule as the volute throat already discussed, and the number of vanes is 
always different from the number of rotor blades Z. a working rule often 
used being Z + 1. Many working designs that follow the rule for L given 
above appear to have a ratio for D^/D 2 of 1.3-1. 5, a constant width b which 
is a little wider than the impeller tip width, and a maximum diffusion angle 
of 11°. The vanes are as thin as practicable, following the type of profile 
shown in Fig. 6.31. Designs using aerofoil sections have been used, but are 
not conventional. Experiments with ‘pipe’ or drilled diffusers for flow that is 
supersonic at impeller outlet have been used but are not general, as the 
simple generated shapes of Fig. 6.31 are easy to produce and control. 
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Figure 6.31 Vaned diffuser. 
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Figure 6.32 Effect of the choice of diffuser system on compressor performance (after 
Petrie (1964) courtesy of the Institution of Mechanical Engineers). 
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Petrie (1964) reports a study for a helicopter gas turbine where four 
diffusion systems were investigated. Figure 6.32 illustrates the effects on 
performance of three vaned diffusers with different inlet angles, and a 
vaneless diffuser. It can be seen that the larger angle and area both moved 
the surge point to higher mass flows and extended the flow range before 
choking. However, the vaneless diffuser lowered both the maximum 
pressure ratio and efficiency, significantly extended the total flow range, but 
did not significantly improve the range of flow from surge to choking at any 
driving speed. The paper by McCutcheon (1978) also emphasizes the role 
the diffuser plays in determining the range of a compressor, particularly 
when the impeller tip Mach number approaches unity. 

In multistage compressors and pumps the passages that guide the fluid 
from one impeller to the next are partly diffusers and partly simple transfer 
ducts. Typical shapes employed are shown in Fig. 6.33 for a multistage 
pump. 

The example of sections 6.3.1 and 6.3.2 is now continued through the 
diffuser. The height of the diffuser passages (dimension b in Fig. 6.31) will 


Return passage 



Figure 6.33 Comparison of the effect of crossover duct geometry on pump 
performance. 
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be 0.045 (larger than the tip width h already found), and the number of 
passages will be 17 (16 4- 1). Diameter D : in Fig. 6.31 will be 0.52 + 0.05 = 
0 57 m, and the radius of the mean width of the throat will be 0.35 as an 
estimate. 

Assuming free vortex flow in the vaneless space, 

357.76 X 0.52 = K u 0.57 

Thus the tangential component of flow velocity at 0.57 m diameter is 
326.38 ms -1 . The radial component of velocity now must be found. One 
approach is to assume constant density and then correct, or to guess the 
radial velocity. Assume that l/ R = 55 ms -1 . Thus the absolute velocity = 
V(55 2 + 326. 38 : ) = 330.98. Assume that stagnation temperature remains 

constant: 

440.37 = T + 330.98 2 /2C p 
T s ,„„ = 385.87 K 

Pi = /?„;(385. 87 /440.37) 3 5 = 2.36 x 10'' Nm ' 7 
The density thus becomes 2.13kgtrT\ and 

V R = 9/2.13(71 x 0.045 x 0.57) = 52.37 
Thus another iteration is really needed, but the solution will be left at this 
point. 

The throat width i in Fig. 6.31 must now be determined. The radius of its 
centre is fixed, as a guess, at 0.35. The process is followed again, with the 
assumptions that the stagnation temperature remains the same and that V’ R 
at the throat radius is 38ms -1 : this gives the width l = 18mm. The as- 
sumption of a diffuser efficiency then allows the solution. Alternatively a 
suitable change in area of the passage may be designed, and the design 
checked to give a pressure rise. The solution will stop at this point. 

6.3 Thrust loads due to hydrodynamic effects 

6.5.1 Radial thrust forces 

In vaned diffuser designs the pressure distribution round the impeller 
periphery remains approximately uniform for a range of flows, and to a 
reasonable degree this is so for vaneless diffusers. For the volute tvpe of 
design it is found that the pressure distribution is approximately constant at 
the design flow, but as the flow departs from design the distribution changes 
as shown in Fig. 6.34. 

Stepannof (1957b) suggests for pumps a simple equation for the radial 
thrust force: 


F r = KpgHD 2 B 2 


(6.11) 
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Figure 6.34 Variation of pressure round an impeller with flow rate in a volute pump. 

Here gH is the specific energy generated by the pump, and the other 
dimensions are shown in Fig. 6.35 for a conventional design. K is an 
empirical constant defined as 

K = 0.36J1 - (<2/0 dcsign ) 2 ] (6.12) 

Other texts give more sophisticated equations that allow for recirculating 
flow effects in the spaces between impeller and casing. All indicate that the 
resulting radial loads to be absorbed by the bearings supporting the shaft 
tend to be large, and in all cases may not be neglected. The formulae are 
used to size shafts and bearings, but it is found in many pumping instal- 
lations that even though the shaft has been carefully sized, with a margin for 
error, shafts do break. These events may be due to poor materials or to 
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Figure 6.35 Definition of dimensions used in radial force calculation. 


misuse or poor fitting. In other cases the failures are due to effects not 
allowed for in the formulae, which are based on steady-state flow con- 
ditions. This circumstance was studied and many researchers agree that, 
since flow conditions are not steady, there is a substantial fluctuating force 
component that must be added to the steady value. In papers by Turton and 
Goss (1982, 1983) it is shown that the fluctuating forces may be of the same 
size as the steady load, and this factor as well as the fluctuating nature of the 
load effects gives rise to vibration and bearing and shaft failure in some 
cases. The reason for the fluctuating components is felt to be the varying 
velocity distribution leaving the impeller and interacting with the casing. 

One device for reducing the radial loads is the double volute; two flow 
paths are provided in the volute so that each receives half the flow from the 
impeller, and the radial thrust loads cancel out. The penalty is increased 
friction loss, which may result in probably a 2% reduction in efficiency. A 
few other solutions exist, but the only tried method is the double volute. 

6.5.2 Axial thrust loads 

Axial hydraulic thrust is the summation of unbalanced forces acting on the 
impeller in the axial direction. Reliable thrust bearings are available so that 
this does not present problems except in large machines, but it is necessary 
to calculate the forces. These forces arise due to the distribution of pressure 
in the spaces between the impeller and the casing. In Fig. 6.36 the as- 
sumption is illustrated for the discharge pressure acting oxer the baekplate 
and front shroud in the single entry pump, and down to the wear rings for a 
double entry impeller; the suction pressure is assumed to be distributed 
across the area up to the wear ring diameter in both designs. It can be 
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Figure 6.37 Unbalanced force due to pressure distribution on a centrifugal impellei 


argued that the pressure forces are balanced in the double suction design 
and this is approximately the case in practice: the single suction impelle 
clearly has an unbalanced force. In practice there is always an unbalance! 
force acting on a double entry impeller due to such factors as unequal flov 
distributions in the two entry passages, and interference due to externa 
features such as bends in the suction line, so that thrust bearings are alway 
needed. Figure 6.37 illustrates the probable pressure distribution on tla 
impeller surfaces. The actual pressure variation will depend upon surface 
roughness of the pump surfaces on the side clearances, and on leakage flow 
through any wear rings fitted, so that front and back net forces may var 
from the design conditions assumed. The classic solutions assume that th> 
fluid in the clearance space rotates at about half the impeller speed as a sola 
mass, and the texts by Stepannof (1957b) and Pfleiderer (1961) may bi 
consulted for the equations. The net result of the pressure forces and th> 
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Figure 6.38 Rear wear ring, balance chamber and balance holes method of axial 
force balance. 


fluid change of momentum is a force towards the suction flange; the force 
magnitude depends on size, outlet pressure, rotational speed, and whether 
the impeller is provided with a shroud, as in many pumps, or is open, as in 
many compressors. The most common wav of reducing the axial thrust load 
is shown in Fig. 6.38. A wear ring is formed on the impeller backplate. thus 
creating a chamber vented to the suction through ‘balance holes’. Suction 
pressure is in this way applied to the backplate, up to the wear ring diam- 
eter. and hence reduces the delivery pressure related force. This increases 
the risk of cavitation inception but is tolerated because it reduces the load. 
Other methods connect the balance chamber to the suction by using a 
balance pipe, or use the so-called pump-out vanes which centrifuge fluid 
back to the casing. These also work quite well but there is a critical clear- 
ance: vane height ratio, as too tight an axial clearance results in vibration. In 
multistage machines these measures cannot be used, and in some designs 
impellers are used back-to-back to balance thrusts, but the usual technique 
is to use balance pistons or balance discs as illustrated in Fig. 6.39 for 
compressors. The sizes of the components, such as wear ring diameter, 
balance piston diameter and balance hole sizes, depend on the calculated 
loads and the probable leakage flows considered acceptable. A common 
design approach is to consider that 10% of the thrust load is taken by the 
bearings and the remainder absorbed by the balance system being used. The 
equations published in the standard texts refer only to steady loads, and the 
fluctuating load problems already introduced in section 6.5.1 again apply. 
As Turton and Goss (1982. 1983) among others have indicated, the fluctuat- 
ing component in some cases tends to be twice the steady value, and the 
shafts and other components should be sized accordingly. A comprehensive 
series of papers on Hydraulic forces in Centrifugal Single and Multistage 
Pumps was given during an Institution of Mechanical Engineers Seminar in 
1986. and this group of papers should be consulted for more detail on thrust 
causes and cures. 


(b) 


(c) 


Figure 6.39 Alternative systems of axial force balance: (a) back to back lavout: (b) 
balance piston: (c) balance disc. 


6.6 Exercises 

6.1 The second stage of an intercooled centrifugal compressor in a helium 
line is illustrated in Fig. 6.40. The design rotational speed is 24000 
revolutions per minute and the static pressure and temperature are 4 x 
10 5 N m -2 and 300 K at inlet to the stage. If the inlet angle at point A is 
30° to the tangential direction determine the mass flow rate and maxi- 
mum mach number for zero incidence at point A if the absolute velocity 
makes an angle of 3° to the axial and its tangential component is in the 
same direction as the rotation of the blades. If the impeller total to total 
efficiency is 90%. determine the discharge stagnation pressure. 

6.2 The double suction water pump shown in Fig. 6.41 is driven at 1750 
revolutions per minute and delivers 0.8 m 3 s -1 . There are seven blades 
per flow path of backward curved design with an outlet angle of 25°. 
Flow can be assumed to be divided equally between the flow paths. 
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Figure 6.41 Exercise 6.2. 


Assuming zero inlet whirl, estimate the Euler specific energy rise and 
the blade angles at points A and B on the suction edge of the blades. 

If the flow conditions change to give an inlet swirl of 5ms -1 in the 
direction of rotation at point A, estimate the reduced specific energy 
rise and the fluid angles at points A and B. 
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6 3 A single entry centrifugal compressor provided with a vaneless diffuser 
delivers air at a mass flow rate of 10 kgs -1 when rotating at 30000 
revolutions per minute. Air is drawn in from a large plenum chamber 
where the stagnation temperature and pressure are, respectively, 290 K 
and 1.25 x 10 5 Nm 2 . 

The impeller diameter is 270mm, the tip width 25 mm, the diffuser 
passage height is 27 mm and its outer diameter is 550 mm, there are 
eight impeller blades 3 mm thick radial at outlet, and the compressor 
efficiency can be assumed to be 79%. 

If the Mach number based on impeller outlet absolute velocity is 
limited to 0.95. construct the outlet velocity triangle. 

Calculate the static pressure and temperature at outlet from the dif- 
fuser using a diffusion factor Cp of 0.5, and discuss the limitations of 
this factor (Cp = Apl0.5pV 2 . Where V is the gas absolute velocity at 
outlet from the impeller). 

6.4 A single stage single entry centrifugal pump is to deliver water at the 
rate of 100 m’ h _1 when driven at 2900 revolutions per minute. The 
pump is illustrated in Fig. 6.42. There are five backward curved impeller 
blades each 3 mm thick, and the nominal outlet blade angle is 20° to the 
tangential direction. 



If a conventional volute casing is fitted, determine the specific energy 
rise generated at design flow and the probable input power needed. 
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Figure 6.43 Exercise 6.5. 


If a vaneless diffuser is fitted with the same volute design placed 
around it. calculate the probable extra specific energy rise generated. 

6.5 Figure 6.43 shows a centrifugal pump impeller. It is to deliver water at 
the rate of 15()m 3 h _1 . If the value of k. is to be 0.75. suggest the 
specific energy rise probable. It may be assumed that the impeller has 
five blades, that the outlet angle is 25\ and that the pump has a 
conventional volute casing. 

For the design shown estimate the axial thrust, if the suction pressure 
is 0.85 x 10 5 Nm" 2 . 

6.6 An air compressor is to be of centrifugal design, is to be driven at 
1575 rad s -i and will draw directly from a plenum chamber where the 
stagnation conditions are 1.25 x lO^Nm -2 and 300 K. 

The impeller has an outer suction diameter of 300mm. a hub of 
100 mm. a tip diameter of 600 mm and a tip width of 20 mm. 

Estimate the maximum mass flow rate if the inlet relative Mach 
number is limited to 0.85. Also determine the probable outer stagnation 
pressure if the outlet radial velocity at limited flow is 100ms~‘ and the 
impeller has 12 radial blades. Assume the efficiency is 82%. 

Suggest how usual slip factors can be improved upon and also how the 
flow effects can be modified. 
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Axial machines for 
incompressible flow 


7.1 Introduction 

Axial fans and pumps may be treated as incompressible flow dev ices, pro- 
ducing large flow rates and small pressure rises. Axial pumps in service 
deliver up to 80000m 3 h~‘ against resistances equivalent to a level change of 
20m, and have been applied to irrigation, draining, dock dewatering and 
circulation duties. Fans have been applied to ventilation duties and cooling 
tower circulation, and commonly produce pressure rises of the order of 
500 Nm" 2 ; special designs give higher pressure rises. Sizes clearly vary with 
duty, from kitchen ventilation fans to machines up to 4 m in diameter. 

Axial turbines may have their guide vane and runner blades coaxial, as in 
the tubular or bulb turbine, or have guide vanes disposed in the radial plane 
and the runner of axial design. Hydroelectric plant using this type of machine 
will be working with level drops from a few metres to over 50 m. producing 
powers from 100 kW to over 100 MW. the maximum runner tip diameter 
being of the order of 8m. 

The principles underlying pumps, fans and turbines of the axial type have 
been introduced in Chapter 4. In the sections following, two examples are 
used to illustrate the way the principles are applied to fans and pumps, with 
the practical implications indicated. 

A brief description of the axial water turbine follows, covering the es- 
sential points concerning the shapes of the hydraulic components. This 
section concludes with simple examples to illustrate typical velocities and 
velocity diagram shapes. 

The chapter concludes with a discussion of the forces acting on static and 
rotating blades, and their complex nature is illustrated. 

7.2 Axial flow pumps and fans 

The machines described have stream surfaces that are a.visymmetric, and in 
their simplest form consists of a rotor only; in some cases guide vanes are 
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added downstream, and in rare applications are used upstream to provide 
control. In many instances the prime mover is an electric motor, so that the 
rotational speeds tend to be synchronous and constant; variable-speed drives 
are rarely used except for the largest machines. 

One approach to design is outlined by using a fairly typical fan application 
with appropriate comments. The references cited may be consulted for other 
approaches and for more advanced techniques. 

A fan draws air in from the atmosphere (pressure HP Nm” 2 and 290 K) 
and is to supply at a rate of 6m J y' against a pressure of 1 .005 x lO" N m” 2 . 
The applications engineer suggests that the probable duct diameter will be 
750mm. 

The hub diameter needs to be defined first. Assuming a hub tip ratio of 
0.5 (a typical value). D h = 375mm. and the axial velocity through the rotor 
ignoring blade blockage effects is 18. 108 ms -1 . The peripheral velocities are 
now to be found, bearing in mind that for noise level limitation tip velocities 
should not exceed 70-75 ms -1 . The driver speed may be assumed to be 
1450revmin”', so the peripheral velocity at the tip is 56.94ms” 1 and the 
hub peripheral velocity is 28.97ms” 1 . 

The density of air at the conditions specified is 1.2kgm"'\ so the specific 
energy rise is 416. 67 J kg” 1 . The Euler specific energy rise, assuming a fairly 
typical efficiency of 80% will thus be 520.8 J kg” 1 . (If the machine were a 
pump the overall efficiency would be found from Fig. 2.6 after the specific 
speed has been calculated. The text by Balje (1981) presents the same sort 
of picture for fans as well, and this may be referred to for additional 
information.) 

Using the Euler equation (1.4). 

520.8 = u&V u 

Considering the hub section first. W u = 16.14ms '. The resulting velocity 
triangles are shown in Fig. 7.1. and the air angles relative to the axial 
direction are /), = 58° and (i 2 = 31.44°. and (S m = 44.72°. If the incidence and 
deviation are ignored these will be the blade angles too, so that 6 = 26.56. 
This will be assumed for the purpose of this example. 





Figure 7.1 Velocity triangles for axial flow fan example. 
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Figure 7.2 Plot of a, and blade data camber angle (Hay, Metcalfe and Reizes (1978) 
courtesy of the Institution of Mechanical Engineers). 


A choice of profile must now be made, and in this treatment the circular 
arc cambered plate will be used, the data for which is given in Fig. 4.5. 
Cambered plates are often used in low performance fans as they are cheap 
but have sufficient strength. Where stress levels are higher, as in high 
performance fans and in pumps, the blade profiles discussed in section 4.3 
must be used. Whatever profile is used. C L and C\ JC D are read from the 
data sheet for a suitable angle of attack that gives a good margin before 
stall, or are related to the maximum lift to drag ratio. These values can then 
be fed into equation (4.7) to find solidity, and then choices would be made 
of blade number, chord and stagger angle. 

An alternative approach that gives an optimizing view is the paper by 
Hay, Metcalfe and Reizes (1978), from which Fig. 7.2 is derived. They 
approached the problem of blade design by producing carpet plots for 
several profiles of a, against (9, where 

tfi = /fi - y (7.1) 

As Fig. 7.2 illustrates, the carpet plot has two sets of curves for C L and for 
C L /C D . From the figure, an optimum value for a, is 4.5°, with C, = 0.95 and 
C L /C D approximately 30. 

The design approach then goes on to correct for the fact that the blades 
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Figure 7.3 Plot of interference coefficient against pitch: chord ratio (Hay, Metcalfe 
and Reizes (1978) courtesy of the Institution of Mechanical Engineers). 


will be in cascade, the appropriate correction factor being shown in Fig. 7.3. 
The angle of incidence cq is replaced by a c , and ft by /ft, with 

etc = /ft - 7 (7.2) 

They relate ft to ft,, and a c by using the equation 

ft = (ft m - «c) - 0/2 + S (7.3) 

Here <> is defined by Howell (d = m0(s/c)"), and for a circular arc camber 
line 


Thus 


m = 0.23 4- ft/500 


ft = (/ft 


etc) - 0 



(7.4) 


Then, assuming n = 0.2, using this treatment and equation (7.4) gives c/s = 
1.48 and s/c = 0.676, and 7 = 44.72 - 4 = 40.72°. 

From Fig. 7.3 the proximity factor k = 0.8, the corrected C L = 0.76, s/c = 
0.451: from equation (7.4), ft = 33.6° and <5 = 5.724°. 
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The method then goes on iterate again, but for this example the numerical 
treatment will stop. If iteration continued a correction to u c of the order of 
the difference in /? 2 would be made, a corrected 6 obtained, C L and the 
other qualities found, resulting in a new value for /? 2 , and so on until the 
differences became very small. 

Reverting to the example, the probable blade angle will be /i 2 = 27.88° 
and - 58°. With c/s = 1.85, using 12 blades gives c = 182; from this a 
simple calculation gives the radius of the circular arc as being 400mm. 

Similar calculations may be performed for the tip and other sections, and 
the blade profiles related to one another as outlined in section 7.4. The 
numerical values will not be given here, however, but the concept of radial 
equilibrium may be used to relate the basic velocity diagrams, and the 
iteration procedure used as needed. Alternative approaches are detailed by 
Eck (1973) and Wallis (1961). With the advent of suites of computer pro- 
grams it has been possible to develop approaches based upon the principles 
outlined and to inject corrections for blockage and for secondary flow 
effects, as well as to build-in procedures for stacking the blade profiles on 
selected lines for manufacture. Such approaches are found in contributions 
by Myles and Watson (1964), Myles, Bain and Buxton (1965) and by 
Hesselgreaves and McEwen (1976), among others. 

The preceding discussion has dealt with rotor blades, but the same pro- 
cedures will follow for stator blades, for the absolute fluid vector angles are 
found in the velocity diagrams and the profiles will emerge in the same way 
as described above. 

The mechanical problems in fans of fixing and location are discussed 
briefly by Eck (1973), and the pump applications are well dealt with in the 
standard text by Stepannof (1957b) and bv other authorities. A recent 
conference on fans updates design principles and gives some interesting 
examples of applications (Woods-Ballard, 1982 and Lack, 1982) 

As a further illustrative example consider an axial flow fan with a tip 
diameter of 1 m and a hub diameter of 0.4 m, which rotates at 1450 revolu- 
tions per minute. Five cambered plates have been used as blades, at the root 
section the chord is 100 mm. and the camber is 2%. For a flow rate of 
10m 3 s' 1 estimate the probable pressure rise, the stagger angle at the hub 
section. Assume that the free vortex principle has been used in design of the 
blades, that the air density is 1.2kgm _ \ and the piates that form the blades 
are bent to a circular arc form. 

From the data given 

flow area = 7r/4 (0.1 2 - 0.4 2 ) = 0.66 m 2 


V A = 10/0.66 = 15.15 ms 1 

"h = 1450 x n x 0.4/60 = 30.37 ms' 1 


Therefore 
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Figure 7.4 Axial flow fan example, geometry. 


The spacing at the hub section is given by 

n x 0.4/5 = 0.251 


so that 


c/s - 0.398 


For a 2% cambered plate C\ = 0.5 at a = 3.5 for best performance (from 
Fig. 4.5), and using geometry for the hub section the camber radius is 
629mm and the camber angle is 9.12° (Fig. 7.4). If zero inlet whirl is used 
the blade inlet angle referred to the tangential direction is /f, = 26.5°. From 
Fig. 7.4 a is 3.5° to the chord line, so the relative velocitv makes an angle of 
incidence of - 1.06°. 


1 4 (> 
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Figure 7.5 Axial flow fan example, velocin triangles. 


The stagger angle will be 90 — 30 = 60°, as shown in the figure. If zero 
deviation is assumed. (i 2 = 34.56° and 

() m = 26.5 + 34.5o 2 = 30.53° 

From the velocity triangles (Fig. 7.5) 

V / u2 = 4 . m s - 1 

and gH E = 4.68 x 30.37 = 142. 17 J kg . 


7.3 Axial water turbines 

Water turbines of this type may be either of the Kaplan propeller type, or 
the fully axial or bulb type (Fig. 7.6). As can be seen the turbine rotor is 
axial in both designs, and the guide vanes (or wicket gates) may be disposed 
in the radial plane or be coaxial with the rotor. 

For both flow paths the problem is the same - the establishment of the 
flow surface shapes and the stream surfaces on which the blade profiles are 
to be disposed, and the design of the blade shapes themselves. The basic 
hydraulic problems are the same as those met in the axial pump, with the 
essential difference that energy is being extracted from the fluid. The meri- 
dional approach to stream surfaces has been covered in Chapter 5, and the 
methods outlined may be used with empirically determined proportions, 
based upon experience and presented against specific speed as shown in 
texts such as those by Nechleba (1957) and Balje (1981). 
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Figure 7.6 Bulb and Kaplan turbines and their efficiency envelope. 


The free vortex approach to radial equilibrium introduced in section 4.5 is 
often used in relating velocity triangles and rotor blade sections, and once 
the stream surfaces have been established the approach outlined already in 
this chapter may be used to determine the necessary geometry and blade 
profiles. The stress considerations are discussed in section 7.4. Turbine 
blades tend to be of long chord and low camber at the tip and relatively 
thick and cambered at the root section, as sketched in Fig. 7.7. The root 
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Figure 7.7 Hub and tip profiles for a water turbine. 



Figure 7.8 Typical guide vane profiles for a water turbine. 


section is dictated by stress considerations and by the need to give a good 
blend of section into the hub attachment; a typical layout is illustrated in 
Fig. 7.7 for a Kaplan design, in which the rotor blades must be rotated over 
a range of angles to give the control illustrated in Fig. 7.6. 

Whether the machine has fixed rotor blades or moving Kaplan runner 
blades, the main control of output is achieved by the provision of moving 
guide vanes or wicket blades. These have alternative shapes as indicated in 
Fig. 7.8 where pit vanes are the more normal shapes used when water 
enters from the spiral casing are shown. The profile must be chosen to give 
good passage shapes for all openings as well as to control the outlet direction 
into the runner system. The nose shape must be such that it will tolerate the 
wide range of incidence it meets as the vane is moved over the control 
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Figure 7.9 Construction for relating guide vane setting to velocity triangle (as 
outlined by Nechleba, 1957). 


range, and in most cases is a 'blunt' large radius which blends well into the 
main profile, as seen in the sketches in Fig. 7.8. When the machine is truly 
axial, with the guide vanes coaxial with the runner blades, there is very little 
difficulty in relating the blade settings. However, when the layout is that 
usually used in the large output Kaplan designs there is a change in direction 
between the radial plane in which the guide vanes are placed and the axial 
direction which is referenced for the runner blade sections. Figure 7.9 
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Figure 7.10 Some draft tubes used in water turbines. 


illustrates this three-dimensional problem, and demonstrates a simple cor- 
rection technique that allows for the stream surface effect on the flow 
direction, and thus allows the setting angle measured in the radial plane for 
the guide vanes a to be related to the Euler triangle angle a'. 

The draft tube or outlet diffuser, fitted after the runner to recover some of 
the kinetic energy leaving the rotor blades, may follow any of the forms 
sketched in Fig. 7.10, variation (c) being often used with vertical axis 
machines. If the cone option (a) is considered to establish the principle, 
applying the equation of energy to the two planes 1 and 2 gives 

2 w2 

— + — + gZi = — + -Z- + °Z 2 + (losses between 1 and 2) (7.5) 

If plane 2 is datum, gZ 2 is zero, and the draft tube efficiency is defined as 
// DT = pressure regain/( V — In - 
equation (7.5) becomes 


P2 ~ Pi 


= gZ 1 + t/DT {V\ 


Vi)/2 


p 


(7.6) 
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Figure 7.11 Typical draft tube and the pressure and velocity distributions. 


Kaplan performed tests on simple draft tube details like cones, bends and 
combinations of simple elements, and later developments have resulted in 
the shape illustrated in Fig. 7.11: typical pressure and velocity variations are 
also shown. In many cases the outlet from the runner is very low, below 
the tail water level in the case of the reversible pump turbine units, 
where cavitation limitations when pumping fix the runner level. The text bv 
Nechleba (1957) outlines some of the history of water turbines and draft 
tubes, and recent work reported at the International Association for 
Hydraulic Research (IAHR) for hydraulic machinery symposia suggests a 
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Figure 7.12 Spiral and draft tube shapes in a barrage installation. 


contribution to the overall losses of about 8% from the draft tube, falling as 
the head drop over a Kaplan machine increases. 

The spiral casing acts as the device that directs the How into the guide 
vanes and thence the runner, so its function is the exact reverse of the volute 
in a centrifugal pump. In many machines the spiral shape is used, but in the 
run of river barrage machines the shapes may be as in Fig. 7.12, the profiles 
being formed in the concrete mass of the dam. 

To illustrate typical values of velocity and the shape of the velocity 
triangles a short solution to a student problem follows. 

Each Kaplan turbine in a large river barrage scheme rotates at 65.2 rev 
min' 1 , has a tip diameter of 8m and a hub : tip ratio of 0.4. During prov- 
ing trials, with a river level difference of Urn, the gauged flow rate was 
500m 3 s' 1 and the electrical output was 45 MW. If it is given that the 
alternator efficiency was 96% and the mechanical efficiency 97%, determine 
the plant efficiency, the hydraulic efficiency and the velocity triangles for the 
rotor root and tip sections, assuming that the blading followed the free 
vortex principle. If the draft tube has an area ratio of 1.5:1, has its exit 
centre line 6 m below the tailwater and its entry section 2 m above tailwater, 
estimate the turbine outlet pressure, and comment upon this figure. It may 
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Figure 7.13 Tip velocity triangle. 


be assumed that the vapour pressure of water is 2500 N m -2 . and the draft 
tube efficiency is 80%. 

The hydraulic power is given by 

g x 11 x 500 x 10 2 = 54 MW 

45 

I/O = err = 0.834 


Thus 


Also 


'/ H 


0.834 

().% x 0.97 


0.896 


gH K u i u .|. = g x 11 x 0.896 = 96. 69 J Kg 1 

= it 1 1% , (since V u2 assumed zero) 


The velocity at the runner tip is u, = 27.31 ms -1 , and at the runner hub l' 
tt h = 10.92ms -1 . Therefore at the tip section 1%] = 3.54ms -1 . Using the 
annulus dimensions and neglecting blade blockage. U, = 11.84ms -1 . The 
tip velocitv triangle resulting is Fig. 7.13 and = 26.48°. /F = 23.43°. 

Assuming a free vortex blade flow system, at the hub U,,, = 8.85ms '. 
Figure 7.14 results for the hub velocity triangle, and (1 { = 80.1°. fi 2 = 47.3°. 
This is based on zero outlet whirl, so V 2 = V a = 11.84ms -1 . 

In Fig. 7.15 the area in plane A is given by (7r/4) (8) 2 = 50.2 m 2 , and that 
at plane B is 75.4 m 2 . Applying the energy equation between A and B. 
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The theoretical energy regain in the draft tube is 


9.95 2 - 6.63 2 
2 


27.52 J kg -1 


The actual regain is 

0.8 x 27.52 = 22. 018 J kg~' 

Thus the draft-tube loss is 5.502Jkg~'. Substitution in the energy equation 
yields gH A = 58. 36 J kg -1 or p A = 0.584 x 10'' NitT 2 . 

Since at normal temperatures the vapour pressure for water is about 
2500 N m~ 2 , cavitation at the turbine runner outlet is unlikely. 


7.4 Forces on blades and their implications for design 

7.4.1 Static blades 

Clearly, guide vanes are subject to varying loads as the angles change during 
the control cycle. When closed, for example as in Fig. 7.16, if Ap is the 
pressure difference across the row of blades, the force exerted F is given by 

F = A ptli (7.7) 

Ap depends upon the inlet pressure and the discharge pressure levels in a 
pump, and upon the difference between the spiral casing pressure and that 
in the turbine space for water turbines, which is a function of the tailwater 
level and whether the space is full of water or empty. The torque needed to 
keep the vanes closed is 

T = FX (7.8) 

This, with an allowance for dynamic or shock effects which is based on 
experience, is used to decide the blade fulcrum pin sizes and profile thick- 
nesses when using acceptable stress levels. 


F 



Figure 7.16 Force on a guide vane. 
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Figure 7.17 Streamlines for open guide vanes. 



Figure 7.18 Variation of guide vane moment with flow rate. 


If the vanes are open, as sketched in Fig. 7.17 for turbine guide vanes, it 
is possible by plotting velocity changes round the profiles to estimate the 
pressure levels and thus to determine the net fluid torque applied to the 
vane about any chosen pivot point. Figure 7.18 illustrates the way this is 
likely to vary, the pivot point being chosen to balance maximum moments 
M\ and M 2 . These values are used to determine pivot and blade thickness 
sizes, and also to design the linkages used to transmit the moving effect from 
actuator to vanes. 


7.4.2 Rotating blades 

In Fig. 7.19 water turbine blades are used again to illustrate the pressure 
distribution round a blade profile and also how the hydraulic load varies 
along the blade length. As well as these loads there are the mechanical 
loads, both centrifugal and tangential, due to the mass and rotation of the 
blade, so that considering the blade as a cantilever will lead to the stresses at 
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Figure 7.19 Diagrammatic representation of fluid forces on a rotor blade. 


all sections to be established. Clearly the root section will be most highly 
stressed, and the profile must be adequate mechanically as well as aerody- 
namically. The stress situation is extremely complex as is the root fixing 
loading and stressing, since in variable geometry blades the loads change 
considerably, resulting in the need for a complete study of the variations in 
stress during the design stage. The article by Barp. Schweizer and Flury 
(1973) illustrates the empirical nature of the problem in water turbine 
design. 


7.5 Concluding remarks 

As in the discussion of centrifugal machines, it is necessary to consider the 
thrust loads that occur due to the interaction of the blades with the fluid. 
Clearly the main effect will be axial thrust. This may be calculated using the 
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equation given as the preamble in section 4.2, or by assuming that the 
pressure change acts over the annulus swept by the blades. In single-stage 
machines the thrust bearing is sized to take this calculated value plus an 
operating margin, or in the case of very large machines a balance disc or 
piston is fitted, as described in section 6.5. It will be noted that if the axis of 
the machine is vertical, the whole weight of the rotating assembly must be 
taken by the thrust bearing in addition to the hydraulic load. In this simple 
approach, complications are met during start-up. where the water column in 
pumps, for example, can give a very large component that may differ in 
direction from the usual load. Stepannof (1957b) and other hydraulic texts 
deal with this problem. 

Since this is not intended to be a treatise on the complete design of 
machine elements, the texts and papers cited may be consulted for further 
detail. For a very useful exposition of the application of computer techniques 
to water turbine design, Chacour and Graybill (1977) may be consulted, as 
can articles in such journals as the Escher Wyss News and Water Power and 
Construction. 


7.6 Exercises 

7.1 An axial flow water turbine is to generate 750 kW when installed in a 
small ‘run of river’ station where the level difference is usually 2.5 m. 
Propose a suitable outline design if the angular velocity of the runner is 
to be 7. 15 rad s -1 , and the tip diameter is limited to 5 m. It may be 
assumed that tj m = 95.5%, ^ a | tC mator = 97% and a/ h = 92%. 

7.2 An axial flow fan is to deliver air to a ventilation system at the rate ol 

4 m 3 s' 1 with a pressure rate of 175 Pa. The ventilation system is designee 
with a mean velocity of 15 ms -1 and the fan intake is at lO'Pa ant: 
290 K. Perform an outline design of a suitable fan. assuming a ‘hydraulic 
efficiency of 75% and a rotational speed of 1450 rpm. Discuss am 
assumptions you may make. Can the factory use a constant section tc 
produce the fan? 

7.3 A small axial water turbine is fitted to a ‘run of river' station. The leve 
difference is 2.5 m. the rotational speed 7.15 rads -1 , the tip diameter b 

5 m, the hub diameter 2.5 m. When rotating at 7.15 rad s -1 . the estimatec 
flow rate was 50m s -! . If 12 Clark Y-section blades with a chord at th<. 
hub section of 500 mm are fitted, estimate the runner blade setting angk 
at the hub (noting that this angle is 90° - the stagger angle), and the 
power developed at the shaft, assuming that hydraulic efficiency is 92°/c 
and ideal flow conditions apply. 
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7.4 Each Kaplan turbine in a large river barrage scheme rotates at 65.2 
revolutions per minute, has a tip diameter of 8 m and a hub tip ratio of 
0.35. When on proving trials with a river level difference of 11m the 
gauged river flow through each machine was 500m 3 s _1 and the electrical 
output was 45 MW. Assuming that the mechanical efficiency was 96% 
and the alternator efficiency was 97%, determine the hydraulic efficiency 
and the velocity triangles at the tip and root sections assuming free 
vortex principles and zero outlet whirl. If the draft tube has an area 
ratio of 1.5:1, has its exit centreline 6m below the tailwater level, and 
has an efficiency of 80%, estimate the outlet pressure from the turbine 
and comment on this level. The water vapour pressure is 2500 N m“ 2 . 



8 Axial turbines 

and compressors for 
compressible flow 


8.1 Introduction 

The principal difference between axial turbines passing gases and their 
incompressible counterparts is the compressibility of the fluid being moved. 
The basic fluid mechanics is the same, and the profile behaviour information 
has already been covered in Chapter 4. The same can be said of compressors 
and blowers. 

Gas and steam turbines are in wide use as the main driver in electricity 
generation, in gas pipeline booster stations, and as prime movers. Gas 
turbines are also in use in such industries as the production of industrial 
process gases, and in the utilization of surplus energy in blast furnace plant 
as an aid to fuel economy. Sizes of machine vary from one or two megawatts 
to 660 MW in the base load stations supplying the grid system in the UK. 
Gas turbines are at the lower end of that range, typical of the size being the 
Tornado machine producing 6 MW currently being produced by Ruston Gas 
Turbines. 

Compressors of the axial type are used in aircraft engines, blast furnaces, 
petrochemical plant, nitric acid production, natural gas iiquifaction systems, 
and in the process industries. Machines have been supplied to give pressure 
rises in excess of 25 bar, and for flow ranges from about 20m's _1 to over 
350m 3 s _1 . The largest power input is of the order of 88 MW to a single unit, 
though a string of machines may be supplied to allow for smaller units and 
for the provision of intermediate pressures. 

Figure 8.1 illustrates a simple axial compressor, indicating the way the 
annulus changes, and Fig. 8.2 shows a gas turbine and the annulus area's 
change through the flow path. The simple steam turbine shown in Fig. 8.3 
illustrates the greater area change needed by the large specific volume 
increase through the flow path. The principles underlying the layout of these 
machines are introduced in this chapter, but for detailed examination of the 
problems the references cited must be consulted. 
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Figure 8.1 Simple compressor. 


One stage 



Figure 8.2 Simple gas turbine. 

8.2 Approach to axial compressor principles 

In multistage machines like axial compressors the total pressure rise is 
shared between the stages, so that an important consideration is the stage 
pressure rise. Other matters of fundamental interest are the blade profiles 
used and their behaviour, the numbers in successive rows, the axial spacing 
between sets of blades and the interaction of the flow leaving one cascade on 
the pattern of flow at inlet to the next. Radial equilibrium, and the effects of 
boundary layer development on the inner and outer annulus on blade 
work distribution, are other factors to be considered in this treatment. The 
principles have already been introduced in Chapter 4. and these will be seen 
applied in the w-orked example which is included in the discussion. The 
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Figure 8.3 Simple steam turbine (in wet steam units the increase in blade length is 
more marked). 


problem of designing profiles, or of applying known information to a given 
situation, has not yet been covered, so that the classical English approaches 
due to Howell (1945) and Carter (1948) will be introduced, and Howell’s 
correlation applied to a sample compressor problem. 

Howell approached the problem by basing his correlation on nominal 
conditions, which he defined as those pertaining to a cascade deflection 
which is 80% of its maximum 'stalling’ deflection. His view was that the 
nominal deflection c* for a number of cascades studied was a function of the 
nominal gas angle a 2 . the space: chord ratio and the Reynolds number. He 
therefore proposed a correlation of these factors with the stalling deflection 
(defined as £ s = 1.25 f;*) and claimed that the camber was apparently not a 
dependent variable. The use of his approach will now be explained by partly 
solving a compressor design problem with explanatory comments where 
necessary. It will be seen that in following this line of solution, the choice of 
the camber angle is arbitrary. 

A compressor is required to deliver air at the rate of 50 kgs -1 and provide 
a pressure ratio of 5:1, the inlet stagnation conditions being 288 K and 
10 5 N m 2 . The target efficiency is 86%. 

pressure ratio = 5 = (1 + 0.86A77r m ) 3 ’’ 

Thus AT = 195.5 and T 02 — 483.5 K. 

The number of stages must now be chosen. In aircraft compressors, to 
limit length and weight, there is a tendency to go for as low a number as 
possible consistent with surge effects due to early stage stall, but in indus- 
trial machines there is a tendency to settle for lower ratios. Choosing here 
10 stages gives 19.6° rise per stage. As Horlock (1958) discusses, the early 
stages may be unloaded to reduce surge risk at low speeds and thus increase 
the temperature rise in later stages, as illustrated in Fig. 8.4; however, for 
the purpose of the example this will not be implemented. 
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Figure 8.4 Temperature rise in late stages. 



Figure 8.5 Velocity triangle for axial compressor. 


Good practice suggests an axial velocity of 150ms" 1 and a mean periph- 
eral speed of 200 ms" 1 (the latter is related to noise generation, which is not 
within the scope of this treatment). 

The application of these values to the first stage, and the assumption of 
50% reaction at mean blade height, allows the velocity triangles to be 
constructed. 


CpAT = QUS\' U 

Howell's work done factor (see Fig. 4.20) is Q = 0.86. Therefore. V,, nu .. in 
114.5ms" 1 . The velocity triangle of Fig. 8.5 results, and 


a, = /A = 15. 
a 2 = P\ = 46.35 : 

£ = 30.45 

V x = 156.05ms" 1 

W, = 217.31 ms" 1 
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Figure 8.6 Tip velocity triangle for axial compressor. 


Since stagnation conditions are quoted, the static conditions at the inlet 
must be calculated to allow density and flow area to be found. 

p | = 0.86 x 10 5 Nm -2 
T ] = 275.88 K 

thus p\ = 1.087 kg m -3 , and the flow area is given by the equation 

A = 50/1.087 x 150 

Therefore the annulus dimensions are tip diameter 0.674 m. mean diameter 
0.472m and hub diameter 0.27m. Thus, since the mean peripheral velocity 
is 200ms~\ the rotational speed is 8093 rev min '. 

If these sizes and speed are not acceptable, owing to space or driver 
choice, another trial must be made and new values determined. 

The mean height triangles were found; now the tip section will be ex- 
amined. Figure 8.6 results from applying the free vortex principle, from 
which the maximum fluid velocity at the tip appears to be 322.59 m s 1 for 
the first stage blades. The acoustic velocity is 

a = V(1 .4 X 287 X 275.88) = 332.94ms 1 

The corresponding Mach number is 0.97. This is rather high, suggesting that 
measures to adjust the velocities need to be taken to reduce this value to 
possibly 0.85 or lower. 

We now return to the mean height section. Howell s correlation will be 
applied to the blade section design. The deflection e*' = 30.45, and fi 2 = 
15.9°, so from Fig. 8.7 the space: chord ratio is approximately 0.7 (fi 2 is the 
outlet angle from the blades in the triangles, and is a 2 in Fig. 8.7). Since the 
blade height is 202 mm, then for an aspect ratio of about 3-4 a chord ot 
50mm may be acceptable, yielding a value for spacing of 35.4mm; thus the 
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Figure 8.7 Deflection against outlet angle with varying sic ratio (Howell (1945) 
courtesy of the Institution of Mechanical Engineers). 


number of blades is approximately 42. Selection of a prime number to avoid 
vibration interaction problems suggests 43 blades, resulting in a chord of 
49.3 mm and an aspect ratio of about 4.1. If a C4 profile is used (as in 
Howell's work) the gas and blade angles may be related with the camber 
line, with assumptions of incidence: using Fig. 8.8 the optimum incidence 
appears to be 2.5°, so the blade inlet angle becomes 46.35 — 2.5 = 43.85°. 
To determine the deviation, equation (4.22) may be used 

d = mdV(slc) 

where 

m = 0.23(2«/c) 2 + 0.1(«f/50) 

as proposed by Howell (a is the distance of the point of maximum camber 
from the leading edge, and «f is the nominal outlet angle). 

Assuming a circular arc camber line. 

m = 0.23(1)' + 0.1(15.9/50) 

to give S = 0.22 0 and (i' 2 = fi 2 ~ 0.22 0. Thus 

9 = ft\ - ji 2 = 43.85 - 15.9 + 0.22 6 

and 6 = 35.8°. The stagger angle is 

y = 43.85 - 35.8/2 - 26° 
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Figure 8.8 Performance of a cascade (Howell (1945) courtesy of the Institution o 
Mechanical Engineers). 


Following Howell (1945) and using his data (Fig. 8.8). then assuming 

C L =0.85 
C D = 0.025 

and using equation (4.20), the probable maximum row efficiency is 94.1. 

The camber radius and other geometric functions may be calculated usin; 
simple trigonometry, and similar calculations performed for other sections 
The sections may be related radially, as discussed in section 7.4, to give th> 
required stress level calculations. 

Carter's (1948) correlation proceeded to relate the optimum lift : dra; 
ratio to the deflection, gas incidence angle and outlet angle. A desigi 
following his approach would be as in Table 8.1. It should be noted tha 
Carter considered variations in outlet angle to be unimportant over thi 
range 0 < fi 2 < 40°. 

Both the correlations and others like them are effectively two dimen 
sional, based on static idealized cascade rig test data, and the only cor 
rections for three-dimensional flow are to use the work done factor am 
radial equilibrium approaches to relate the blade sections. One school o 
thought is to use the actuator disc technique pioneered by Hawthorne am 
developed by Horlock (1958). This lends itself to computer-based design, a 
is well illustrated by Railly (1961) and Railly and Howard (1962). Th> 
alternative approach is to use data obtained from actual tests on machines 
McKenzie (1980) used a low speed four-stage compressor with a hub : ti] 
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Table 8.1 Steps in the calculation of optimum performance, following Carter (1948) 


• Estimate a value of optimum incidence t opt ; thus 

optimum inlet angle = gas angle + / op , 

• Find the optimum deviation angle using the equation for deviation and Fig. 3.21 

• Calculate the outlet blade angle from gas angle and deviation. 

• Find the optimum incidence / opl for an isolated aerofoil. Carter (1948. 1961). 

• Calculate C L for each blade. 

• Find the change in / opl from foil in cascade to isolated foil: 

/=;_/ — 

1 ‘ppt 'opl(original) 2y/ c 

where C„ is the induced velocity perpendicular to the mean velocity due to unit 
circulation and pitching, found using Betz rules as f(s/c)(y). 

• Find i l>pI = i„ pl (<,riginai) + >• and compare with the value estimated. Iterate if 
necessary. 


ratio of 0.8, provided with constant section blades with C 5 profiles on 
circular arc camber lines. He showed that for constant stagger angle the 
maximum efficiency occurred at a constant flow coefficient irrespective of 
the camber and space: chord ratio, from which he derived a correlation of 
stagger angle to mean gas angle: 

tan . = tan fi m - 0.213 (8.1) 

He showed from exhaustive traverse tests that the deviation followed the 
equation 

d = (2 + e/3)(s/c) l!3 (8.2) 

He compared actual data available for low hub: tip ratios giving alternative 
information, and also proposed two design charts which could be used. He 
also made a comparison with the fan approach of Hay. Metcalfe and Reizes 
(1978) discussed in Chapter 7. and found fairly close agreement. He also 
considered the difficulty of accurately predicting stall, which the other 
approaches share, and demonstrated considerable scatter for the perfor- 
mance plots covered; he remarked that low speed performance prediction 
is reasonable but is of doubtful validity w'hen applied to high speed com- 
pressors. The approach is as valid as the Howell and Carter correlations, 
however, as it is based upon actual machine data, and may be applied to 
compressors in the same speed, flow and pressure range. 

The first stage has been examined, and the design of the rest of the blade 
rows would proceed in the same way, with allowance being made for density 
change - the blade heights must reduce progressively. Some designs use an 
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annulus based on constant hub diameter, used in some aircraft designs, 
where a ‘waisting’ effect accommodates auxiliaries in the area of the outlet 
casing. Many others keep the outer diameter constant to simplify stator 
blade fixings; blade sections are constant over a number of stages, the length 
being changed to suit annulus height requirements. This reduces costs with 
very little reduction in efficiency. 

When the blade profile dimensions have been determined, the axial spacing 
of the blades, the root fixings, and the stress levels are found. The axial 
spacing has to allow for good fluid coupling between the blade rows and at 
the same time provide for relative expansion between the rotor assembly 
and the stator; a typical industrial axial spacing appears to be up to 25% of 
the blade chord, but one company has stated that, on low pressure rise 
stages, varying the spacing from 14% of chord to as much as 50% gave 
changes in efficiency that were considered to be commercially acceptable. 
Tip clearance too must be chosen to minimize flow loss but still allow for 
radial growth of the rotating assembly into the stator. Blade stressing was 
covered in Chapter 7 in outline, and root fixings using pins or fir tree root 
types are matters of company policy. Carchedi and Wood (1982) illustrate 
how a commercial design may be carried out with the working restrictions 
on blade profile choice, and how general mechanical design decisions are 
reached. 





Figure 8.9 Simple impulse stage, typical velocity triangles and the ideal overall 
efficiency envelope. 
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8.3 Axial turbine principles 


8.3.1 General principles 

The basic principles were introduced in sections 1.5.2 and 1.7.3, as were the 
usual parameters used. Before proceeding to discuss the correlations used 
for loss and performance prediction, it is necessary to cover the simple 
relations governing the enthalpy change in a stage with the velocity triangles. 

Consider first the 'impulse' stage, illustrated in Fig. 8.9. The nozzle 
discharges gas at an angle of a, and this impinges on the rotor in such a way 
that the gas angle is the blade inlet angle (deviation and incidence are 
ignored). All the enthalpy drop is assumed to take place in the nozzle. The 
diagram work is 


itAV u = 2it(V, cos a, - it) 

(8.3) 

The energy available is V]/2 (assuming no loss in the 
to the enthalpy drop as shown in section 1.7. and F, 
velocity). 

The diagram efficiency is 

nozzle; this is related 
is then the isentropic 

2/d l 7 | cos a, - u) 
>,n ~ Vj/2 


4 " / it \ 

"■> = r | cosa| - vj 


Differentiating with respect to it Vfi and equating 
maximum ;/ D occurs when the blade /speed ratio is 

to zero suggests that 

cos a | 

K/V 1 = , 

(8.4) 

the maximum efficiency is 


no = cos 2 a. 

(8.5) 

and at optimum speed ratio 


work output = 2ir 

(8.6) 


As Fig. 8.9 also shows, the diagram efficiency varies with the speed ratio, 
reaching a maximum and then falling away again as the speed ratio con- 
tinues to increase. 

If stages are placed in series, in the velocity compounded design where 
only the first stage nozzle is an energy convertor, the succeeding stator 
rows are designed to change velocity direction only (called velocity 
compounding). Figure 8.10 demonstrates how the efficiency curves change, 
and also the maximum efficiency, as the number of impulse stages increases. 



1 70 Axial turbines and compressors for compressible flow 


<1 stage 



Figure 8.10 Effect of the number of stages on the efficiency envelope. 


Now consider the 50% reaction case. By similar reasoning, and following 
the velocity triangle shapes shown in Fig. 1.11. the optimum speed ratio is 
given by the relation: 


— = cos a | 


The maximum diagram efficiency is given by: 

2 cos 2 a | 

'Id = ~ J 

1 + cos a | 


and 


(8.7) 


( 8 . 8 ) 


work output at optimum speed ratio = tr (8.9) 

Figure 8. 1 1 . following Kearton (1958) compares the efficiency envelopes for 
multistage reaction and multistage impulse machines based upon a nozzle 
angle of 20° and a nozzle efficiency of 94%. 

This discussion, together with those in Chapter 1. allows the construction 
of velocity triangles. The next stage in analysis or design consists in relating 
the gas directions to blade angles, and giving an approach to the estimation 
of losses and efficiencies, and allowances for three-dimensional flow and 
other effects. 

The main correlations are based on data obtained using perfect gases such 
as air; all gas turbines, and many steam turbines, pass such fluids, because 
superheated steam follows the same laws as air. The only exception to this 
rule is the turbine passing wet steam, as in the low pressure cylinders of a 
large main load machine and in machines such as those used in nuclear 
generators; the particular problems associated with these machines will be 
discussed. 
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Figure 8.11 Comparison of multistage impulse and reaction machine efficiency 
envelopes. 


Horloek (1966) in a useful background study, brought together two main 
streams of information in comparing the correlations from steam turbine 
and gas turbine technology. He compared the early steam turbine data 
reported by Kearton (1958). Stodola (1945) and Gun (1939). for example, 
with the gas turbine studies published by Zweifel (1945). Ainley and 
Matheison (1957) and Soderberg (1949) among others. Horloek demon- 
strated the wide differences between the approaches reviewed, and a private 
study comparing them with some company data, indicated an almost 10% 
band in efficiency prediction for low reaction and 'impulse' blading designs. 

Zweifel (1945) used data obtained for both accelerating and decelerating 
cascades with typical large deflections, and argued that for typical zero 
reaction blading a tangential blade loading coefficient y/ T should be used 
rather than the more usual pressure coefficient or lift coefficient. He also 
adv anced a correlation based on profiles with a value of <// T = 0.8, and used 
it to predict the number of blades that w'ould give effective performance for 
axial turbocharger turbines. 

Ainley and Matheison (1957) suggested that the total loss of a cascade of 
blades consisted of profile loss T F ' secondary loss > s and tip clearance loss 
Y ( (Fig. 8.12). They correlated T P for zero-incidence flow conditions with 
blades having a thickness: chord ratio of 20%, and applied corrections for 
incidence and thickness. They proposed a combined formula for Y s and Y c : 


Y s + V'c 


;. + 



C, \ 2 COS 2 « 2 
,s/cj cos ’a m 


( 8 . 10 ) 


B and <) are defined in Fig. 8.13, and / is also shown. The data were 
presented for low reaction blading, and in correspondence Ainley suggested 
the corrected line as being a better correlation for higher reaction turbines. 
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Figure 8.12 Loss coefficient proposed by Ainlev and Matheison (after Ainley and 
Matheison, 1957). 


Soderberg, using Zweifel’s information and other data, related losses on 
the basis of space: chord ratio, Reynolds number, aspect ratio, thickness: 
chord ratio and profile geometry, assuming zero incidence and correcting for 
flow conditions with some incidence and for other deviations from the basis 
on which the data was obtained. Horlock (1966) may be consulted for a 
thorough comparison of this and other approaches, and demonstrates the 
wide differences between them and the essential complexity of the problem. 

An alternative approach is to use information obtained from machines to 
complete the correlation. One example of this approach is that of Craig and 
Cox (1970) for steam turbines. They identified two loss groupings. Group 
one comprised losses due to profiles, secondary flows, clearance effects, 
annulus and cavities for both stator and rotor rows; group two losses were 
due to gland leakage, windage, wetness effects, and partial admission. The 
group one effects are effectively corrected for, as in the approach of Ainley 
and Matheison (1957), with company experience fed into the databank. The 
discussion of the group two losses highlights the problems facing the steam 
turbine manufacturer, for example in the exhaust casing systems. These 
must carry the wet steam leaving the low pressure turbines into the con- 
densers with the maximum effectiveness in the minimum axial distance. 
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Figure 8.13 7. factor proposed by Ainley and Matheison (after Ainley and 
Matheison, 1957). 


owing to the need to keep the total length of the turbine to an acceptable 
value. Flow is strongly three dimensional, and steam turbine manufacturers 
have used extensive air model testing to improve their designs. 

The comparisons of Horlock referred to above, and the later study bv 
Dunham (1970). related to the profile data then available, which was ob- 
tained at Reynolds numbers based on chord of about 10 A and at Mach 
numbers well below sonic velocities. As power output has increased, annulus 
dimensions and turbine cylinders have also become larger, with low pressure 
turbine blade tip velocities attaining supersonic values. This, and the cost 
pressure causing designers to reduce machine sizes and optimize the number 
of stages, has led to higher velocities, and the work of designers has gone 
hand-in-hand with analytical and experimental studies. 
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Figure 8.14 Effect of partial admission. 


Wilson (1987) approached the preliminary design of axial flow gas turbines 
with the Craig and Cox correlation in mind. He identified twin objectives: to 
correlate efficiency wfith reaction, loading coefficient (defined as uAVJtc) 
and flow coefficient; and to provide an informed choice of vortex flow 
distribution across the turbine annulus and allow an optimum annulus 
profile to be chosen for turbines working at various expansion ratios. He 
had as a further objective the improvement of the Craig and Cox loss 
estimation method, already discussed, along with other published ap- 
proaches. The paper draws extensively from the earlier work already dis- 
cussed in this chapter and proposes some simple guidelines for reaction 
choice, hub tip ratio, stage number, and loading coefficient. He also suggests 
an equation for swirl variation to give an almost constant nozzle outlet angle 
across the turbine annulus. The paper ends with loss correlations in terms of 
pressure loss ratio statements based on Craig and Cox and upon Brown 
(1972). 

The economic and engineering aspects of steam and turbine design are 
well illustrated by the contributions of Burn (1978). Parsons (1972) and 
Smith (1975), the latter being an extremely well illustrated treatment of wet 
steam problems. Hesketh and Muscroft (1990) detail the engineering of the 
steam turbine generators for the Sizewell B nuclear station giving an insight 
into the problems to be solved and the reasoning behind the solutions 
selected in constructing the plant. Recent studies of wet steam problems are 
discussed by Laali (1991) and Tanuma and Sakamoto (1991). 

For further reading refer to Fielding (1981) which is a look at losses from 
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a different standpoint from those quoted above; Hill and Lewis (1974) is an 
experimental study of a very high divergence turbine; a study of secondary 
flow losses is provided by Gregory-Smith (1982) and. in later papers by the 
same author and his colleagues. 

8.3.2 Partial admission problem 

Partial admission is the term used to describe the flow situation where only a 
proportion of the annulus is supplied with gas, so that the rotor blades are 
moving from a live throughflow sector into a ‘dead' zone where the gas is 
carried bv the rotation from the live flow and its energy is dissipated in 
windage and other losses. 

Partial admission occurs in many turbocharger turbines, owing to the need 
to keep discharging pipes lront individual cylinders of the engines being 
blown from interfering with one another. In this case a different sector of 
the annulus is live at different times in the engine cycle, so that the flow is 
ver\ three dimensional and complex, giving rise to large losses. In small 
steam turbines a complete sector may be blanked off permanently to adjust 
steam flows to give a set output, so carry-over fluid flow losses from the live 
to dead sectors account for a large part of the losses. Evidence of this is 
given in Fig. S. 14. based on work published by Stenning (1953) and by Suter 
and Traupel (1959). Clearly seen are the effects of greatly reduced active 
flow annulus area. 


8.4 Other problems 

Since the temperatures are high there are large differential expansions, so 
that clearances, which are cold during erection, have to allow' for the closing 
up due to expansion. This is particularly important with tip clearances, 
where adjustment is needed to prevent tips touching, particularly for com- 
pressor blades made from titanium where local fires may be caused. In gas 
turbines the problem is one of simple calculation owing to the light con- 
struction of the casings and blades, though arrangements have to be made to 
allow blades to take up their own attitude under load and expansion move- 
ment to reduce induced stresses. In steam turbines the high pressures require 
very heavy pressure vessel construction: the rotor heats up more quickly 
than the casing, so that the axial clearances during warm-up are more 
critical than those when cold or at running temperature. 

Aircraft turbine temperatures are now exceeding 1600 °C so that elaborate 
cooling methods using hollow blades are now standard practice. Industrial 
machines like the Ruston Tornado are now utilizing the technique, as higher 
efficiencies are sought and extended life is a high priority. 

Axial thrust loading is clearly a problem, and in industrial compressors 
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with many stages is usually compensated for by introducing a balance piston 
or disc as described in section 6.5.2. In aircraft gas turbines and industrial 
gas turbines the thrust produced by the compressor is partly balanced by 
that from the turbine so that a thrust bearing is only needed for small thrust 
loads and for axial location. Where a turbine is used as a drive other 
measures are needed. For example, in multicylinder steam turbines the high 
pressure and middle pressure turbines are arranged so that steam inlet ends 
and the axial thrusts from the two units are opposed; the low pressure 
cylinders have their steam inlet belt in the centre with expansion in both 
axial directions, thus ensuring thrust balance. Gas or air turbines driving 
compressors use balance as in the gas turbine, but need thrust bearings in 
most other cases. As in the machines described in Chapter 7, radial thrusts 
are low and well within the capacity of the bearings conventionally used. 


8.5 Computer-aided solutions 

There are a number of CFD and CAD programmes now in use. Some are 
simply the provision of computer power using existing empirical data and 
the approaches noted above, but the reader is referred to a number of very 
recent papers and articles that indicate the power and capability of the soft- 
ware being developed. Typical are the papers by Scrivener et al. (1991). 
Hart et al. (1991) and Kobayashi et al. (1991) who were all contributors to a 
recent Institution of Mechanical Engineers conference. 


8.6 Illustrative examples 

8.6.1 Axial compressor example 

An axial flow compressor is to deliver helium and has six stages equalh 
loaded, with stage temperature rise 25 K. If the overall efficiency is to be 
87% , determine the overall pressure ratio if T (ll is 288 K and p m is 1(F N m _; . 
and the stage pressure ratio for the last and first stages. Calculate the blade 
height at entry to the last stage and the rotational speed if F, is 165 ms -1 
and makes an angle of 20° to the axial direction, the mass flow is 13 kgs -1 
and D mc;in is 680mm. The Howell work done factoi is 0.83 and the mean 
section has been designed to be 50% reaction. What is the maximum Mach 
number in the last stage? 

From the data given 

T m = 288 + 6 x 25 = 438 K 
C p = 5.193 x 10 ; J kg“ 1 KT 1 
M = 4.003 



Illustrative examples 177 


thus 


and 


R = 8 ' 314 = 2.077 x 10 3 J kg -1 K“‘ 

4.003 


1 V -\ 


C v = 3.116J kg - K 
Q> , 

n = — = 1.67 


n - 1 


= 2.493 


Therefore 


POZ 


P 01 7 overall 


= 1 + 0.87 X 


150\ 


2.493 


288/ 


= 2.54: 1 


p t)2 = 2.54 x ION 5 m“ 


At the first stage 


Pul | 

p { \ i / first stage 


and. from Fia. 1.17. 


25 \-- J93 

1 + °- 88 x m ) 


//,, = 0.88 


1.2 


At the last stage 7",,. at inlet to this stage, is 438 - 25 = 413 °C and 

^ = 1.138 

Am 

The stagnation pressure at inlet to the stage is 

2.54 x lO-Vl.138 = 2.232 x 10 ? Nm“ 2 
To construct the velocity triangles at the sixth stage: 

25 x 5.193 x 10 3 = 0.83 x aAl u 

Turning to Fig. 8.15 

F u , = 60.055 

V A = 165 sin 20 = 56.43 ms-' 

F u: = (u - 60.055) 

AF U = (u - 60.055) - 60.055 
thus. V u = u — 120.11, and 

25 x 5.193 x 10 3 = 0.83 x u(u - 120.11) 




Figure 8.15 Axial compressor example - velocity triangles. 


Solving this equation gives u = 460.055 and 
. 460.055 x 2 x 30 

0.68X, 


V u2 = 400 ms' 1 
V 2 = 432.7 ms" 1 = W, 


The static temperature and pressure at inlet to the last stage are now- 
needed: 


413 


r 5 + 


175. 6 2 

2 x 5.193 X 10 3 ' 


p s lpo s = (410/413) 2 493 = 0.98, p s 


r 5 = 410 K 


2.192 x 10 5 Nm“ 2 


thus 


_ 2.192 x 10 s 

p5 ~ 2.077 x 10 3 x 410 


0.257 kg m" 3 


and 


<2 = 


13 

0.257 


= 50.5 ms 


Since the annulus height is small compared with the mean diameter 


50.5 

56.43 


= n x 0.68/t, 


h = 0.419m 
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Thus the tip diameter is 1.099m and the hub diameter is 0.261 m. 

The maximum velocity through the stage is Wj and, since the acoustic 
velocity a is 

a = 1.67 x 2.077 x 10 3 x 410 = 1192.52 
the maximum Mach number is 


432.7 

1192.52 


0.363 


8.6.2 Axial flow gas turbine problem 

An axial flow gas turbine has 16 stages, equally loaded and is supplied with 
gas at the rate of 12.5 kgs -1 . Determine the shaft power, the inlet tip 
diameter and the gas angles at the root and tip sections of the first-stage 
rotor. It may be assumed that the rotational speed is 5000 rpm, the inlet 
static conditions are 1150K and 7 x 10 5 Nm -2 , the outlet static conditions 
are 1.05 x 10 5 Nm -2 and 750 K, the turbine rotor inner diameter is lm for 
all stages, the reheat factor is 1.02, is 90%, the axial velocity is 100ms -1 
in all stages, and all root sections are designed to the 50% reaction principle. 

For the combustion gas products n = 1.333 and C p = 1.145 x 
10 3 Jkg -1 K -1 . For isentropic expansion across the turbine the gas laws give 

T ow = 703.36 K 


from the 16th stage, thus 

A/! 1S en,rop,c = (1150 - 703.36) x 1.145 x 10 3 
= 511.4 kJ kg -1 

A/i acIua i = (1150 - 750) x 1.145 x 10 3 
= 458kJkg -1 

The static to static efficiency is 

458 

= 5T77 = 0896 


Since 


R = 


Coverall 

7stage 


7 stage 


0.896 

1.02 


0.878 
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Since 


reheat factor = number of stages X ^- 1| * :nlr 2Bi c P cr . st ^8 e 

overall A/j lsentropic 

i 02 = 16 Ahisen,r °P ic per stage 
511.4 

A/tisentrop.c per stage = 32.6 x 10'Jkg -1 
and 


^actual P er stage = 32.6 x 0.878 = 28.62 x 10 3 Jkg 1 
At inlet to the stage 


7.5 x 10 s 
287 x 1150 


2.27 kgm" 3 


the annulus area is 

m/pV = 0.055 m 2 

and, since D h = lm, D tip = 1.034. Also, u u = 261.8ms" 1 and u t = 
270.7 ms" 1 . 

At the hub section reaction = 50% , which gives 

28.62 x 10 3 = 261.8 x AV a 
AV 7 ,, = 105.73 ms" 1 

The velocity triangles are shown in Fig. 8.16, from which 
V ul = 183.77. V u2 = 78.05 

/?,= 37.97°, (i 2 = 61.45 

At the tip section, since V u x R = constant, 

V ul = 177.73, V u2 = 75.48 

(},= 37.05°, fi 2 = 60.64 

The shaft power is given by 

shaft power = 1.145 x 10’ x A/i actua | x ij m 

= 1.145 x 10 3 x 28.62 x 10 3 x 0.9 
= 29.49 MW 


8.7 Exercises 

8.1 A reaction steam turbine rotates at 1500 rpm and consumes steam at the 
rate of 20000 kgh" 1 . At a certain stage in the machine the steam is at 
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Figure 8.16 Axial flow gas turbine problem - velocity triangles. 


a pressure of 0.4 bar with a dryness fraction of 0.93, and the stage 
develops 0.25 MW. The axial velocity is 0.7 of the mean blade section 
peripheral velocity and is constant through the stage. The reaction at 
the mean radius is 50%, and the blade profiles are the same for both 
rotor and stator blades with the outlet angles being at 20° to the per- 
ipheral direction. Determine the inlet annulus dimensions for this stage, 
if the total losses are estimated to be 30% of the developed power. 

8.2 A four-stage axial steam turbine is designed on the assumption of equal 
loading per stage and a constant axial velocity through the machine, 
with a mean diameter for all stages of 1.5 m. A reaction of 50% was 
used for each stage at the mean section with the nozzle angle at this 
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diameter of 20° and a rotor blade inlet angle of 55° to the peripheral 
direction. The design reheat factor was 1.04, the rotational speed 
3000 rpm. When tested under steam inlet conditions of 220 bar and 
973 K. the exhaust conditions were 0.5 bar and 373 K. For these condi- 
tions calculate the steam mass flow rate, the diagram work, and the 
overall efficiency, if the blade height at the last stage is 35mm. 

8.3 A 10-stage axial flow gas turbine is to run at 3000 rpm when passing gas 
at a mass flow rate of 70 kgs -1 . The gas is that normally produced by 
hydrocarbon combustion. It is assumed that the blading is designed to 
the 50% reaction concept at the mean diameter of 1.5 m. the diagram 
output is to be 16.5 MW, the axial velocity of 100ms -1 is constant 
through the machine, and all stages are equally loaded. If the inlet 
stagnation conditions are 7.5 x HP N m -2 and 1000 K for the first stage, 
determine the annulus dimensions and the mean blade height velocity 
diagrams. Assuming free vortex conditions across the leading and trail- 
ing edges determine the tip and hub velocity triangles and discuss 
whether a constant section could be used, specifying which section 
would be the most suitable. 

8.4 The first stage of an axial flow compressor passing air has a hub: tip 
ratio of 0.6, and F, = V A = 140 ms -1 . The rotational speed is 
6000 rpm, the temperature rise is 20°C. and the Howell work done 
factor is 0.925. The inlet stagnation conditions are 1 .01 x HP N m -2 and 
288 K and the stage efficiency is 89%. Find for the flow condition that 
the maximum value of Mach number (relative) is 0.95. the tip diameter, 
the blade angles at the tip section, the mass flow rate, the stage pressure 
ratio, and the blade angles at the root section. 

8.5 A 12-stage axial flow compressor is to compress helium gas at the rate of 
350 kgs -1 . The overall pressure ratio is to be 3.75, the stages are equallv 
loaded, the design efficiency is to be 85%. and the inlet stagnation 
conditions are 5 bar and 350 K. Determine for the fifth stage the hub 
diameter, the tip and hub velocity triangles, and the blade angles and 
stagger angle at the tip section if the tip chord is 150 mm and the 
incidence angle is 2°. It may be assumed that the rotational speed is 
3000 rpm, the tip diameter is 3.5 m and constant through the machine, 
the blading follows the free vortex principle and is 50% reaction in 
design at the tip section, the axial velocity is 120ms -1 and constant 
through the machine, and circular arc camber lines have been used with 
a radius of 500 mm at the tip section. 
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Radial flow turbines 


9.1 Introduction 

The flow path in radial turbines may be centrifugal as in the Ljungstrom 
design or, more commonly, centripetal as in most of the machines found in 
both water power installations and in gas power generation plant. The two 
types of flow path are illustrated in Figs 9. 1-9.3. 

The water turbine in Fig. 9.1 is a typical Francis design, and may be 
installed vertically as shown or horizontally depending on the station. The 
head or level range applied to these machines at present is from 10 m to over 
500m. the power developed ranging from 10 kW to over 450 MW. with 
runners varying in diameter to well over 6m. Some of these machines are 
used in pumped storage schemes as will be described in Chapter 10, but do 
not differ in essentials from those used purely for power generation. 

The inward flow gas turbine in Fig. 9.2 may be used in turbochargers, in 
small gas turbines, or as expanders in the cryogenic industry. They are 
rugged in design and. particularly when fitted in turbochargers, may rotate 
at speeds over 100 000 rev min -1 . Their flexibility in terms of flow range and 
power range control is much less than the axial turbine designs, so they tend 
only to be used for small power installations. 

The outward flow machine is less common, but has been used to some 
effect with steam as the large increase in specific volume is more easily 
allowed for. The Ljungstrom turbine (Fig. 9.3) is one example of this 
design, and has been supplied to generate powers of up to 100 MW. but is 
offered for small inhouse generating plant delivering 10-35 MW typically, 
and for pass-out turbine applications. 

In the following sections these machines will be discussed in the order of 
this introduction, with a worked example to assist understanding. 

9.2 Water turbines 

The early designs proposed by Francis were purely radial, with both guide 
vanes and runner blades disposed in the radial plane, and with the rotor 
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Exhaust 



Figure 9.3 Outward flow radial turbine. 


having no shroud as in the Kaplan and propellor turbines. As heads and 
powers increased, the flow' path became completely three dimensional, with 
the runner passages being mixed flow as Fig. 9.1 indicates. The machine is 
still known as the Francis turbine. 

Francis machines are used for a wide range of level differences from 50 m 
to well over 500m, which is into the Pelton impulse turbine range. Figure 
9.4 shows how the elevational profile changes with specific speed, the low'er 
value being associated with the highest level drop through the turbine. Since 
the flow path is three dimensional, the stream surfaces may be established in 
any of the ways outlined in Chapter 5. Some method like that discussed in 
section 7.2 for the Kaplan is used to relate guide vane setting angles to the 
runner inlet angles, and also the inlet and outlet angles of the rotor blades, 
since there is a considerable curvature of the stream surfaces and of the 
leading and trailing edges shown in Fig. 9.4. The shapes of both guide vanes 
and runner blades follow the same principles as those in the axial turbines, 
and the routine step by step method that may be used is well illustrated by 
N'echleba (1957). The more recent computer-based approach is discussed 
by. among others, Pollard (1973) and Chacour and Graybill (1977). The 
house journals of the leading hydroelectric machinery manufacturers may 
also be consulted for articles on the application of computer methods to 
machines. An example of this is the article by Dubas and Schuch (1987). 
This details the advanced techniques used for analysing the pressure distri- 
bution and resulting stress levels in a Francis turbine runner, and predict the 
probable vibration levels. 

Since the only control method is by the use of the guide vanes, the 
efficiency envelope is ‘peaky’, as Fig. 9.5 shows for one guide vane setting. 
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(a) 



(b) 


Figure 9.4 Illustration of the typical Francis turbine runner shape: (a) for low 
specific speeds; (b) for high specific speeds. 



Figure 9.5 Efficiency envelope. 


The envelope is not so wide as that of a Kaplan machine, but the maximum 
efficiency compares very well, as Fig. 9.6 demonstrates. 

Operating problems are discussed in the house journals and in such 
journals as Water Power and Dam Construction. Typical of the former are 
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Specific speed 

Figure 9.6 Comparison of the Francis and IFR efficiency envelopes. 


papers by Grein and Staehle (1978). Anghern, Holler and Barp (1977) and 
Grein and Bachmann (1976) who detail the development problems as- 
sociated with these machines. 

9.2.1 Francis turbine problem 

Figure 9.7 illustrates a Francis turbine. When the angular velocity was 
29.3 rads' 1 and the flow rate was 12m ’s' 1 the total energy at inlet to the 
guide vanes was 500 J kg' 1 , the pressure loss through the guide vanes was 
estimated to be 0.3 x 10 3 Nm' 2 . the pressure loss through the runner 
passages was estimated to be 0.15 x 10 5 N m~ 2 . and the pressure at inlet to 
the runner was found to be 3.4 x 10’’ Nm -2 . The pressure at inlet to the 
draft tube was measured to be 0.86 x 10 s Nm' . Estimate the hydraulic 
power and the necessary setting angle for the guide vanes referred to the 
tangential direction, if the hydraulic efficiency may be assumed to be 94%. 

Draft tube inlet area = 1.25 2 /4 = 1.227 m 2 . Therefore. 

V 3 = 12/1.227 = 9.78 ms' 1 

Energy levels through the flow path: 

total energy at inlet to guide vanes = 500 J kg' 1 
static energy at inlet to runner = 340 J kg' 1 

static energy at inlet to draft tube = 86 J kg' 1 
Applying the energy equation to the guide vane passages 

500 = 340 + V\!2 + 0.3 x 10 5 /10 3 + 1.5g 
V 2 = 15.185ms' 1 
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Figure 9.7 Francis turbine problem. 


Similarly, across the runner. 


340 + 


15.185 2 


9.78 : 


4- l.5g = 86 + — — + 1 g + 


0.15 x 10 5 


in 3 


+ work done 


thus 


work done = 31 1 J kg 1 

The hydraulic power = work done x mass flow = 3.74 MW. Since rj hyd = 
0.94: 


and 


shaft power = 3.51 MW 


gti e = = 330.85 = u 2 V u2 - u 3 V u3 

Taking zero outlet whirl conditions, since u 2 = 23.44 ms -1 

V u2 = 14. 11 m s -1 

The inlet velocity triangle becomes as shown in Fig. 9.8, and 
a 2 = 21.69°, /? 2 = 31.02° 
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Figure 



Using the setting angle correction shown in Fig. 7.7, 

a ' = 20.49° 


9.3 Radial inflow gas turbine 

This machine, used extensively in turbochargers and in cryogenic appli- 
ations. is illustrated in Fig. 9.2. The flow path is centripetal, and unlike the 
Francis design the majority of turbines have fixed geometry nozzles and 
runners, the output being varied by altering gas flow rate in cryogenic 
systems and by controlling fuel rate in gas turbines. The machines are, like 
the Francis turbine, typically single stage, having a similar efficiency envel- 
ope. and are much used in low flow applications where their overall ef- 
ficiency is better than the equivalent axial machine. 

Typical idealized velocity triangles and h-s charts were discussed in 
section 1". as were the common efficiency definitions. The conventional 90° 
inlet angle for the rotor blades gives rise, when zero outlet whirl is assumed, 
to large outlet angles in the exducer section (Fig. 9.9). The velocities also 
tend to be large, as the worked example that ends this section shows, and 
this leads to one of the limiting considerations discussed in section 9.3.2. 

Approaches to loss correlation and design criteria follow those outlined in 
Chapter 8 for water turbines, and studies by Benson (1970) and by Rohlik 
(1968) supplement the working formulae proposed by Balje (1981). 

9.3.1 Nozzle systems 

Two nozzle systems are available, guide vanes and a vaneless space. In the 
former case the guide vane angle seems to be in the range 10-30°, but there 
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Figure 9.9 Sketch of an IFR turbine showing the twist in the exducer. 



Figure 9.10 Effect of increasing guide vane throat area on the mass flow coefficient. 
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is need for correction to allow for (low in the vaneless space between the 
vanes and the rotor, as Knoerschild (1961) has discussed. The most import- 
ant parameter is the throat area (the minimum area), and this is the ap- 
proach used by turbocharger manufacturers. For example. Fig. 9.10 shows a 
plot of pressure ratio to mass flow parameter for four nozzle areas: it is 
typical, and indicates the way that the mass flow range of the machine is 
affected bv nozzle design and geometry. The way in which turbine and 
compressor matching is carried out using changes in geometry in the nozzles 
is discussed in a number of papers presented at an Institution of Mechanical 
Engineers conference on Turbocharging held in 1978, and the various 
methods are well illustrated. 

Where there are no vanes the flow follows the Archimedian spiral path 
approximately, the direction being modified by real fluid effects, but the 
height of the nozzle passage is used to control the area imposed on the gas 
in the same way that the nozzle ring is manipulated. The smaller turbo- 
chargers use the vaneless designs as they are cheaper to produce, and give a 
wider range of flows even though the efficiency is somewhat reduced. 

9.1.2 Rotor geometry 

Apart from the flow direction, the flow problems met are the same as those 
found in centrifugal compressors. The blade angle at inlet is normally 90°. 
but secondary flow is a problem as in the compressor (Fig. 9.11) so that 
there is a balance needed between control and friction loss due to too manv 
vanes. Jamieson (1955) Wallace (1958) and Hiett and Johnston (1964). 
among others, have considered the choice of vane numbers. Jamieson sug- 
gested the minimum number to avoid flow reversal in the passages, and Fig. 
9.12 is based on his correlation, suggesting the minimum number of vanes 
against nozzle angle. Wallace used a more rigorous approach, and com- 
ments that more vanes are needed than Jamieson's minimum. Both ap- 
proaches suggest a large number of vanes at the inlet, but flow blockage 
considerations mean that only half the blades continue through the exducer, 
the others stopping short as illustrated in radial compressors. Hiett and 
Johnston (1964) reported tests on a rotor with a nozzle angle of 77° and a 
12-blade rotor that gave a total to static efficiency of 84% measured at 
optimum flow conditions. They found that doubling the number of blades 
gave an improvement in efficiency of 1%, which suggests that the data of 
Fig. 9.12 form a reasonable compromise. 

Wood (1963). following the lead of Balje (1981). discussed the limits on 
performance imposed by the flow conditions that resulted from the zero 
outlet whirl design concept. Fie argued that choking would occur when the 
maximum outlet relative velocity approached the local speed of sound. 
Wood produced the plot (Fig. 9.13) of stagnation pressure ratio against a 
ratio (VA/Vo) 2 , showing two lines, one for a Mach number of unity and one 
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Need for more blades 


at this radius 



Figure 9.11 Need for flow control at inlet to an IFR rotor. 



Figure 9.12 Variation of blade number with nozzle angle (after Jamieson, 1955). 


for a more typical value of 0.7 which is better for use in machine design. The 
velocity V () is the so-called spouting velocity, defined as 

V 0 = V[2C p (7oi - 7 o 3ss )J (9.1) 

(the two temperatures are defined in Chapter 1). The plot is for air as the 
working fluid, and an efficiency of 90% was assumed, but the principle is 
well illustrated. It may be noted that this attempt at correlation is based on 
the design condition of zero outlet whirl, and that there is very little in- 
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Overall 



Figure 9.13 Correlation proposed by Wood (1963; courtesy of the American Society 
of Mechanical Engineers). 


formation available on the probable effects for other flow conditions that 
may occur. 

It is not possible in this discussion to do anything but introduce the basic 
theory of machines which are the essential component of modern high speed 
turbochargers. A recent contribution to an analytical approach to optimizing 
such turbines was written by Chen and Baines (1992). Using data and 
techniques published earlier by their co-workers and others, they present an 
optimization procedure for a preliminary design study which depends on the 
proposition that minimum exit loss is associated with zero outlet whirl from 
the rotor, and design proceeds from a specified loading coefficient. 

9.3.3 Worked example 

An inward flow' radial gas turbine is provided w-ith a nozzle ring designed to 
give a nozzle angle of 25° to the tangent. The impeller diameter is 160mm, 
the exducer maximum diameter is 115mm and the hub: tip ratio at outlet is 
0.2. The rotor vanes are radial at inlet, and the design was based on the 
50% reaction concept. When running at 50 000 rev min -1 the inlet stagnation 
conditions were found to be 3.6 x 10 5 Nm -2 and 1100K. Assuming a nozzle 
efficiency of 93% a total to static efficiency of 82.5%. suggest the probable 
mass flow rate, the exducer blade angles at root and tip, and the maximum 
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Figure 9.14 Worked example - flow path. 



Figure 9.15 Worked example - velocity triangle. 


Mach number at outlet. It may be assumed that the axial velocity in the 
outlet is 150ms -1 (C p = 1.145 x 10 3 Jkg -1 ). 

The flow path is shown in Fig. 9.14, and it is assumed that there is zero 
outlet whirl: u 2 = 418.9 ms -1 ; ljp = 301.1 ms -1 : ity hub = 60.22 ms -1 . 

Following the velocity triangle in Fig. 9.15 and assuming that there is no 
deviation of flow from the radial direction, the velocity V 2 from the nozzle is 
found from the equation 

V 2 = V(/7 n 2CpA7' aclua |) 

From the triangle, V 2 = 462.2ms -1 and W 2 = 195.34ms -1 . Thus ATac^i = 
100.31 K, and since 50% reaction applies the total temperature drop over 
the turbine is 200.62 K. Therefore, since T in = 1100K, 

t] T s = 0.8 = 200.62/1100 - T 3s 
T 3s = 849.23 K 
T m = 899.38 K 
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Figure 9.16 (a) Tip triangle; (b) Hub triangle. 


With isentropic expansion, following Fig. 1.26, (k/k —1=4) 

Pm = (M 4 = / HOP V 
Pi \tJ \ 849. 23/ 

and 

/849 '>3\ 4 

Pi = 3.6 x 10 5 - 1.28 x 10 s Nm" : 


Since 13 = 150 ms 


T m = 73 + 150 2 /2C P 
T y = 889.56 K 


py 


1.28 x 10 5 
287 x 889.56 


0.501 kg nt’ 


Thus 

m = 0.051 x 150 x ?r/4(0. 1 15 2 - 0.02875 2 ) 

= 0.732 kgs" 1 

For the exducer section, the velocity triangles are as in Fig. 9.16. Consider 
the maximum diameter first: fix = 26.48° and W 3 = 336.39ms~'. Since the 
acoustic velocity is 

a = 1 .33 x 287 x 889.56 = 582.71 ms H 


the maximum Mach number is 0.577. 

For the hub diameter. /j 3 = 68° and W 3 = 161.64ms -1 . 

The change in fU along the exducer trailing edge is typical of the usual 
twist in the vanes needed at the rotor outlet. 
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9.4 Ljungstrom or radial outflow turbine 

The working fluid in this machine is usually steam, which enters the centre 
and then flows outward through rows of blades that are arranged so that 
alternate rows rotate in opposite directions (Fig. 9.3). The two discs each 
drive an alternator. The arrangement offers a compact design of short 
length, but because of overhang problems has not been developed to 
provide the large powers that the axial machines currently deliver. The 
principles and practice of the design are covered by Kearton (1958) and a 
simplified discussion only will be given here. 

If it is assumed that all blades have the same profile, camber and spacing, 
apart from the first and last rows, that the velocity ratio is the same for all 
rows, and that the clearance radially between the rows is neglected, a simple 
outline of theory may be attempted. Consider the three rows shown in Fig. 
9.17, and the velocity triangles that are also drawn there. The inlet triangle 
is based upon the assumption that there is zero inlet whirl, and the others 
have been drawn using the assumptions just stated. The triangles have been 
drawn with the blades in the manner shown to emphasize the interrelations. 


Ring 2 
W 2 in 



U 0 


Figure 9.17 Build-up of velocity triangles for three rings of moving blades in a 
Ljungstrom-type machine. 
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Appyling the Euler equation for the first row, the energy change, since 
V u() is zero, is 

§H ] = U \V ul 

For the second row the energy change is 

gHi = U 2 V u2 - (-t/,V ul ) = U 2 V u2 + t/,V UI 
Similarly, for the third row. 

gH 3 = U 3 V u3 + U 2 V u2 
and for each succeeding row up to row n, 

gH„ = U„V un + U„- l V un _ ] 

The sum of all the energy transfers is the total energy change from inlet to 
outlet. If the energy transfer for a stage for maximum utilization is con- 
sidered, it can be shown that this condition is met for a radial stage when 
U/W = cos a/2 and gH (maximum) = 2 U 2 , which is the same as for an 
impulse machine of the axial type. This demonstrates that, whereas two 
rows of blades are needed in an impulse axial turbine, only one is needed in 
the Ljungstrom machine. 

Early designs were only provided with a radial path, with the height of the 
blade rows increasing with radius. However, with high superheat available, 
modern designs are provided with one or two stages of axially disposed 
blades at the maximum radius of the discs. A number of machines of up to 
60 MW were installed by the UK Central Electricity Generating Board, but 
present designs offered by a number of companies appear to be of a size 
tailored to in-house generating plant for large factory sites; typical powers 
having been outlined in the introduction to this chapter. Kearton describes 
the ingenious techniques used to attach blade rows and allow for lateral and 
radial relative expansion and for the very long blades needed at outlet with 
radial stress and movement. 


9.5 Exercises 

9.1 A mixed flow water turbine of the Francis type installed with its axis 
vertical, is to rotate at 18 rads -1 . Estimate the power developed and the 
overall hydraulic efficiency of the machine given the following in- 
formation relating to dimensions and flow losses. 

• Guide vanes: discharge angle 20° to the tangential direction elevation 
of mean passage height above tailwater level 4 m. 

• Runner: inlet diameter 3300mm, mean outlet diameter 2150mm. 
blade inlet and outlet angles referred to the tangential direction are 
85° and 25°, respectively, inlet flow area 3m 2 outlet flow area 3.1 m 2 . 



1 98 Radial flow turbines 


• Draft tube: inlet 3.25 m above tailwater level, outlet area 11.5 m 2 . 

• Tests at the duty flow rate indicate that losses related to velocity 
were: guide vane loss, 0.05 x outlet absolute velocity energy; runner 
loss, 0.2 x outlet relative velocity energy; draft tube loss, 0.5 x inlet 
absolute velocity energy. 

9.2 An inward radial flow gas turbine has been designed on the 50% 
reaction principle to give 100 kW at 50000rpm when the inlet gas 
stagnation pressure was 3 x 10 5 Nm~ 2 and the stagnation temperature 
at inlet was 850 K. When on test at the specified power the stagnation 
outlet temperature was found to be 670 K and the nozzle efficiency and 
total to static efficiencies were calculated to be 97% and 86%, respect- 
ively, and there was zero inlet whirl. The rotor outer diameter was 
165 mm, the rotor blades were radial at inlet, the mean diameter of the 
nozzle throats was 175 mm, the passages 30 mm high and the blockage 
due to the blades was 15%. 

Determine for the conditions stated (a) the gas mass flow rate, and 
(b) if the axial velocity at outlet is measured as 300 ms 1 find the outlet 
tip Mach number based on the relative velocity. 

9.3 An inward flow gas turbine designed for the 50% reaction case is 
provided with a 150 mm rotor and is to rotate at 50000rpm. The nozzle 
guide vane outlet angle is 35° referred to the tangential direction and the 
rotor blades are radial at inlet. 

When operating at the design point flow rate of 0.56 kg s~' the inlet 
stagnation conditions were 3 x lO^Nm -2 and 850 K, and the nozzle 
and total to static efficiencies were calculated to be 97% and 83%, 
respectively. 

For this operating condition, assuming zero outlet whirl, estimate the 
diagram power, the outlet pressure and temperature and the necessary 
nozzle guide vane peripheral throat area. 

9.4 An inward flow gas turbine is provided with a vaneless casing and a 
rotor that is 150 mm in diameter. The outlet section has an outer 
diameter of 105mm and a hub: tip ratio of 0.25. The design approach is 
to use 50% reaction, and the blades on the rotor are radial at inlet. 
When delivering 120 kW at 50000rpm the casing inlet stagnation con- 
ditions were 3.5 x 10 5 Nm -2 and 1000 K. If the inlet casing ‘nozzle’ 
efficiency is assumed to be 93%, the total to static efficiency is 83% and 
the shaft mechanical efficiency is 80%, suggest the exducer angles at the 
tip and hub sections and the maximum Mach number at outlet. Assume 
that that the inlet casing delivers the gas to the rotor at an average angle 
of 35° referred to the tangential direction. If a zero reaction design 
approach is used with the geometry nozzle velocity and power output 
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quoted above, suggest the effect on the exducer section angles and gas 
flow through the rotor. 

9.5 Figure 9.18 shows part of the longitudinal section of a reversible 
pump-turbine in a hydroelectric station. When pumping, the unit 
delivers 74m's -1 under an entry head of 318 m when rotating at 
333 rpm with an efficiency of 88.5%. 

When generating, the flow is to be 76.5 m’s -1 with an efficiency of 
93.5%. The draft tube efficiency when turbining is 80% and when 
pumping the friction factor / in the formula gH = 4 fLV : 2D may be 
assumed to be 0.005, V being the mean between the inlet and outlet 
draft velocities, and the diameter D the mean of the inlet and outlet 
diameters. 

Determine the minimum pressure at inlet to the draft tube in Nm -J 
and suggest the value of the Thoma cavitation parameter for this 
condition. 

If a tenth-scale model is used to study performance suggest, when 
operated in the pumping mode, the flow rate, head rise, and power 
required, if the model speed is limited to 3000 rpm. 

9.6 Figure 9.19 illustrates a Francis turbine installation. When com- 
missioned the machine was supplied with a total head above atmos- 
pheric level of 100m in the spiral casing, the runner speed was 
29.3 rads -1 and the generated power was 50 MW. The computed ‘wire 
to water' efficiency was 90%. Determine the flow rate, and the pressure 
at inlet to the runner and the draft tube. 
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It may be assumed that the draft tube efficiency is 83% and there is 
zero inlet whirl, the guide vane setting relative to the tangential direc- 
tion is 25°, and the loss through the guide vanes is estimated to be 0.06 
x absolute kinetic energy at inlet to the runner. 



Special machine 
applications 



10.1 Introduction 

This concluding chapter discusses the effects of viscosity change, of fluid 
properties, and of gas content on the operation of pumps, the ways in which 
output control of machines can be achieved, the steady state interaction 
between pumps and systems, and the hydraulic problems to be solved in 
pumped-storage systems. 


10.2 Problems involved in special pumping applications 

10.2.1 Gas suspension problems 

Pumps find difficulty in coping with gas in suspensions and tend to vapour 
lock and cease to pump when gas fills up the suction area. Merry (1976) 
illustrates the effect on pump performance of air in measured quantities, 
and the dramatic effect on the curves is clearly demonstrated. A general rule 
often followed is that if air is present in a proportion above about 10% the 
pump performance is unacceptable. In the process industries, problems arise 
in plant when the pump is shut down for a period and heat soaks from the 
pump case into the fluid. The local temperature rises to a level that causes 
vaporization and a local pocket of gas forms; this causes the effect known as 
vapour locking, and the pump will not pump. Cooling passages in the pump 
casing are usually provided to avoid the problem with those chemicals 
susceptible to this problem. In other cases provision is made to ensure that 
fluid is not retained in the casing. 

Where, as part of the system operation there is a risk of air ingestion, as 
in site water-level maintenance pumping equipment, for instance, special 
provision is made to ensure that air is evacuated from the pump, continu- 
ously if necessary. Arrangements which do this are known as priming sys- 
tems, and are of various types. The simplest consists of a tank in the suction 
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line, and the pump draws from its bottom, so if flow is stopped enough 
liquid is retained in the tank to allow the pump to attain its normal suction 
pressure and thus stimulate flow from the suction source. More sophisticated 
designs have a special pump casing, with a large chamber in the suction zone 
which retains liquid. This liquid is drawn into the impeller on start-up and is 
delivered to the volute from where it recirculates, after giving up entrained 
gas, to the suction to entrain gas from the empty suction line, and then 
passes through the impeller. This process goes on until the suction line is full 
of liquid (the pump is primed) and the pump gives a normal delivery flow 
rate. There are a number of patented recirculation systems, described by 
Rachmann (1967), and reference to this contribution and to pump hand- 
books may be made for more detail. These types are called self-priming as 
they rely on the dynamic action of the impeller. 

A special pump design which is inherently self-priming is the regenerative 
pump (Fig. 10.1). The rotor has radial channels machined in opposing faces, 
and the casing has channels provided as can be seen. Flow spirals in a 
manner similar to the fluid flywheel. The pump as applied to bilge pumping 
aboard ship is an example of self-priming action. Papers by Burton (1966), 
Yamasaki and Tomita (1971) and Tomita, Yamasaki and Sasahara (1973) 
may be consulted for further reading. 

For many years contractors’ pumps have been provided with external 
priming systems. Typically the diesel engine driving the pump provides air 
extraction through its exhaust manifold. A cut-off valve operates when the 



Figure 10.1 Typical regenerative pump. 
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pump is primed, preventing water ingress to the engine. Jet pump extraction 
systems are used and, in vehicle-mounted fire pumps, liquid ring pumps 
(Fig. 10.2) are used to provide prime. 


10.2.2 Liquid-solid suspension pumping 

Gravel, limestone, coal and power station ash are now pumped as a suspen- 
sion in water over considerable distances. In addition, sewage and chemical 
suspensions need to be pumped. Some of the solids are aggressive and some 
of the liquids are corrosive, so that pump parts need to be able to resist 
abrasion, corrosion and erosion. The pump designer has therefore to 
provide for wear, to give a suitable service life and an additional problem is 
to design to avoid blockage. 

Figure 10.3 illustrates typical solids handling impellers designed to avoid 
blockage: if a solid mass can get into the impeller it will pass through. 
Typical centrifugal machines are similar to clean-water pumps but. with the 
special impellers mentioned, the provision of easy access (split casings) and 
easily replaced casing components in areas of high wear. A special type of 
machine (the 'Egger' Turo design) with a recessed impeller is shown in Fig. 
10.4. 

Pumps may be constructed of ‘Ni-hard’ or of chromium-molybdenum 
cast iron, which can be hardened after machining for longer life. Linings of 
rubber, nylon, polyurethane and more exotic materials have been used, but 
the consensus appears to be that they are not often cost effective, so that 
metallic materials tend to be favoured. Ceramic sleeves are used in some 
designs to protect the shaft from abrasion under the seals, and give good 
service: the material is not suitable in moving parts owing to its brittle 
nature. 

Reference may be made to contributions by Burgess and Reizes (1976). 
Odrowas-Pieniazek (1979), Ahmad, Goulas and Baker (1981). Willis and 
Truscott (1976) and Bain and Bonnington (1970) for background infor- 
mation. For current developments the continuing series of hydro-transport 
conferences may be consulted. 

10.2.3 Effect of viscosity change 

The discussions in previous chapters have been devoted to fluids that have 
viscosities similar to that of water. Figure 10.5 illustrates the fall-off in 
efficiency and pressure rise as viscosities increase. The rapid fall-off in 
efficiency begins when kinematic viscosities exceed 100 centistokes. This 
effect arises due to increased friction, so that as well as the loss in perfor- 
mance there is a heating effect that can have undesirable effects on the seal 
svstem. 



Suction 
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Figure 10.5 Effect of viscosity change on constant speed centrifugal pump character- 
istic behaviour. 


10.3 Pumped storage systems 

The size of the generating units supplying electricity grid systems has 
steadily increased in the pursuit of efficiency, so that the base load is now 
being supplied in the UK by a few large stations, some of which are nuclear. 
The efficiency of these units is high, but the fuel burn of a thermal station is 
higher at half load than at full load, so that it is not economical to run at 
anything but full load. This poses problems for the grid control and planning 
engineers, for as Fig. 10.6 shows, the grid demand varies throughout the day 
from very low levels during the night period to two peak periods during the 
normal day. There is an additional problem posed by sudden changes in 
demand, such as those resulting from sudden drops in ambient temperature 
or from surges in consumption at the end of popular TV programmes. 
Keeping plant in spinning reserve to respond is expensive, and response is 
slow due to the need to avoid thermal shock. One solution to the problem is 
the installation of pumped storage hydroelectric plant, where power from 
the grid at times of low demand is used to pump water into high-level 
storage reservoirs, which is then allowed to flow through turbines at time of 
peak loading on the grid. Two advantages are offered - fast response to 
demand, and the use of off-peak electricity to provide storage capacity. 
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Figure 10.6 Typical load variation on the National Grid over 24 hours. 


Early installations used the four-machine layout, that is. a turbine, an 
alternator, a pump and a motor: both hydraulic machines were arranged to 
run at the optimum efficiency. Many plants have been used with very high 
level differences, so that multistage pumps provide storage, and Francis or 
Pelton turbines absorb the water power available. More recently three- 
machine layouts have been used, both vertical and horizontal installations 
being employed, with a single-electrical machine doubling as alternator and 
motor, and using separate turbine and pump maehines. Typical of these is 
the station at Blaenau Ffestinniog in North Wales, provided with a Francis 
turbine designed to give about 79 MW and a multistage pump to lift water to 
a storage reservoir about 300 m above the lower lake. This machine string is 
vertical, with the pump below the turbine for cavitation submergence 
reasons: a coupling is installed that allow's the turbine to generate without 
having the drag of the pump, which is dew'atered and filled with air during 
the power generation cycle. In other sets the pump is always run in air 
during the generation phase, but in some not disconnected. Other sets are 
horizontal. For example, the machines at Vianden, which form phase one of 
the station, are placed with the electrical machine between the turbine and 
the pump, and a small Pelton turbine is provided to reduce electrical loading 
during the start-up period. The operating sequence is shown in Fig. 10.7. 

Alternative machine layouts are illustrated in Fig. 10.8. As can be seen, 
the two-machine arrangement permits lower excavation costs and less 
expensive machinery because only one hydraulic unit is used. Extra elec- 
trical equipment is needed: a phase changer is required to change over the 
direction of the turbine-pump when moving from the pumpine mode to the 
generating mode. There is a compromise required between efficiency when 
pumping or turbining. and a number of studies have been published on the 
problem. The pump-turbine, when generating, has a lower efficiency than a 
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Figure 10.8 A comparison of the machine layouts for pumped storage stations: (a) 
two-machine system - pump plus turbine; (b) pump which operates in reverse as a 
turbine: (c) turbine which operates as a pump - the usual single hydraulic unit system 
used. 


Francis machine under the same level drop and peak efficiency occurs at a 
lower value of pressure coefficient but the efficiency ranges being com- 
parable. The condition for this is that the machine was designed as a turbine 
but capable of being run as a pump. A machine designed as a pump but run 
as a turbine has about the same turbine efficiency as the reversible machine 
just discussed but a very restricted range of effective operation, and tends 
only to be used for high level change applications where several pump stages 
are required. In practice there is a choice of the ratio of pumping and 
turbining flow rates to optimize the operation in both modes, and this is 
linked with grid loading patterns to accommodate economically the periods 
of peak and low demand. A good discussion of this and related problems 
will be found in Meier (1966) Meier et al. (1971) and Muhlemann (1972). 
The Biannual proceedings of the International Association for Hydraulic 
Research (IAHR) section for hydraulic machines, cavitation and equipment 
may be consulted for up-to-date information. 


10.4 Comments on output control of rotating machines 

Turbines, whether water, gas or steam driven, are used to a large extent to 
generate power at a constant speed because of the need for frequency 
control. Marine propulsion units, and those used for some forms of traction 
drives, work over a range of running speeds. The control of steam and gas 
turbines is achieved either by fuel scheduling or by suitable combinations of 
valve settings. Water turbines use guide vane control, and in Kaplan 
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machines finer control is provided by the use of variable settings of the 
runner blades, as has already been outlined in Chapters 8 and 9. 

Output control of compressors and blowers can be provided by speed 
control or variable stator or rotor geometry. Speed control is the simplest, 
particularly if the machine is large and driven by a steam turbine, but the 
drive is expensive if a thyristor-controlled AC motor is used or if a Ward 
Leonard DC set is fitted The alternative systems, however, make the 
compressor more expensive, though stator control can give good control 
without a large penalty in efficiency loss over a range of flows down to about 
60% of the design value. Variable rotor-blade geometry can give good 
control over a wider range of flows, down to perhaps 30% of design, but at 
much greater expense and complication since the control system needs to 
provide a connection between the static actuator and the rotating assembly. 

Pumps and fans may be simply controlled by using valves in the outlet 
system, but at the expense of power since this method creates control by 
increasing pressure loss. Stepannof (1957a) surveyed the problem and 
demonstrated that such machines may effectively be controlled, with 
significant power savings, by using variable geometry inlet guide vanes. 
Figure 10.9 illustrates how varying the vane angles affects the inlet velocity- 
triangles, the important vector being the inlet whirl component of the 
absolute velocity. The Euler equation (1.4) states that 

gH = u 2 V u2 ~ t<iF ul 

The zero inlet whirl ‘design’ condition suppresses the second term on the 
right-hand side of the equation, and if the whirl component is in the op- 
posite direction to rotation (prewhirl) the result in the equation is to in- 
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Figure 10.10 Sketch of pumped system curves superimposed on the constant speed 
curve for a centrifugal pump. 


crease the energy rise, since the niV ul term changes its sign and adds to the 
other part of the equation. There is also a small change in axial velocity, as 
can be seen from the figure. If, on the other hand, the whirl component is in 
the same direction as the peripheral velocity the effect is to reduce the 
energy rise with a small reduction in flow rate. Stepannof demonstrated for 
a fan that a vane movement of 53° from the neutral position, producing 
whirl in the direction of rotation, reduced flow to 70% of the design value 
and the power to about 70% of the rated level. 

Flow control in fan-supplied systems leading to energy savings was 
considered by Woods-Ballard (1982). He showed that a fixed-output fan 
may cost £1200 and an inverter motor drive for speed control would cost 
£4500 involving a payback time of 2.6 years in terms of the saved energy 
costs at 1982 prices. The cost of an inlet guide vane system could be of the 
order of £1750 with a payback time of about 1.3 years. The application was 
a variable air-volume type of air conditioning system, and indicates the 
order of relative costs involved when choosing the best method of providing 
variable flow rate and energy input. 

Similar considerations apply for pumps, though most are driven at 
constant speed, so that changes in system resistance illustrated in Fig. 10.10 
are achieved by valving in the discharge line to bring the crossover point to 
the level needed. Larger pumps, for example boiler feed pumps, tend to be 
driven by steam turbines, and output control is provided with variable- 
geometry rotor blading where the economics allow for the extra first costs. 
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1.1 See Fig. A.l. 

oj = 71 = 157.1 rads -1 , u x = 5.89ms -1 , u 2 = 23.96ms -1 . 

From inlet A, F R1 = V\ = 5.89 ms -1 , therefore 

Q - 5.89 x 0.049 x n x 0.075 = 0.068 m 3 s -1 

gH E = 23.96 V u2 , V U 2 from outlet A = 21.45ms -1 (since F R2 = 
1 .45 ms -1 ) and gH E = 513.91 J kg - 1 . Therefore 

A p = 513.91 x 10 5 = 5.139 x 10 5 NrrT 2 

Hydraulic power = pQ x gH E - 34.95 kW. 



Figure A.l Solution 1.1. 
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Figure A. 2 Solution 1.2. 



Figure A. 3 Solution 1.4. 


(b) 


1.2 See Fig. A. 2. 

Mean diameter = 1.125m, therefore u m = 25.31ms -1 
From the outlet triangles. V u2 = 3.865 ms -1 . therefore 

gH E = 25.31 x 3.865 = 97.82Jkg -1 

1.3 From Exercise 1.2 V A = 10 ms -1 , so that 

Q = 10 x ^(1.5 2 - 0.75 2 ) = 13.25 mV 
Ap = 97.82 x 1.2 = 117. 38N m -2 


1.4 See Fig. A. 3. 


oj = 1450 x ^ = 151. 84 rads -1 
U u p = 151.84 x | = 151.84ms -1 


Ap = 5 x 10 3 Nm -2 
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Figure A. 4 Solution 1.5. 


therefore 

gH E = 5 x 10 3 /1,2 = 4i67Jkg _I = 151.84 l/, j; 
V u2 = 27.44 ms' 1 
(j, = 19.91°, (i 2 = 23.85° and 

Q - 55 x ^(2 : - i.i. 8 : ) = 145.15 nr s _l 


therefore 


power = 4167 x 145.14 x 1.2 = 0.726 MW 

1 .5 See Fig. A. 4. 

co = 45rads _l , t/ tip = 33.75ms' 1 , T T hub = 13.95ms' 1 . U mca „ = 
23.85 ms -1 and gH = 8 x g. Therefore 

gH E = 0.92 x S x g = 72.2.1kg' 1 

72.2 = 23.85 V uI , therefore V ul = 5.03 ms' 1 . /f, = 25.66°, (i 2 = 22.75°, 
flow rate = 14.65 m 3 s' 1 . 

Electrical power = 87 x 14.65 x 10 3 x 0.95 
= 1 .092 MW 


2.1 80% of design speed = 185 x 0.8 = 148 rads ‘. 

4.5 x 10 5 A P . s 

iPPSTP " Tip X &• ,here,ore Ap = 2te x ia N m " 

0.28 O 

I85D 5 * 148D 3 ’ therefore Q = 0.224 try's - 
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At 1 85 rads' 


P = 4 ^ 10 f x 0.28 x 10 3 = 0 . 107MW 


10 3 


0.85 


At 148 rads 


p = 2.8 _ 8x10 ! x 0.224 xl0^ aQ76Mw 


10 3 


0.85 


2.2 


4.5 x 10 5 _ 4,5 x 10 s 
185 2 x D\ ~ 148 2 x DV 

0.28 _ Q 
185 D\ ~ 148 Dl’ 

„ 0.28 x 10 3 4.5 x 10 5 

P = — r-rz — x — — — = 0.108 MW 


D 2 

— = 1 . 23: 1 

D, 


Q — 0.28 itPs 


0.85 


10 3 


2.3 


20 = 


H 


2900V0.03 

(H) 


, 3/4 


/2900 \ 4/3 

^-^-V0.03j = 73.54 m per stage 

820 

number of stages = j = 11.15 (i.e. 12 stages) 


2.4 


37.5 x 10 6 


0.93 


18 


90 2 x Dp 2 NmfDm 1 


Q 


and 


18g x Q x 10 3 

228.35 m 3 s“ 1 

37.5 x 10 6 45 x 10 3 

10 3 20 3 D/P “ 10 * NmfDm 5 


P, = 37.5 MW, Hi = 18m, rj x = 0.93. N i = 90rpm, Hi = 6 m, Pi = 
45 kW and >h = 0.93. 

Hi _ H 2 Pi _ P 2 Qi _ Q 2 

N\D T N\D\' pN]D\ pN\D 5 2 A \D\ N 2 D\ 

and 


ks 


ojVP 

= constant 
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are all valid if losses are neglected (implied by rj x = rj 2 ). Then, 


18 

90 2 D\ 

90 x V(37.5 x 10 6 ) 

(18g) 5/4 


6 

A l\D\ 

/V 2 V45 x 10 3 
(6g) 5/4 


N 2 = 658 rpm 


Therefore 


Di 
D 2 

45 x 10 3 


18 6 
90 2 D\ ~ 658 2 D\ 

658 
90 


= 12 . 66:1 

6 , 


)- 

= 10 3 X Q 2 x 6g 


0.93 

Q 2 = 0.822 m 3 s' 1 


2.5 D,/D 2 = 3:1, Ni 


= 1450 rpm, <2i = 

_ 9.89 x 10 s 
lh _ 287 x 283 

- io 5 

Pl ~ 287 x 298 


5m 3 s 1 , A pi = 450 and ij t = 0.78. 
= 1.218 kg m -3 


= 1.169 kg m 3 


From tables, i'j = 14.5 mm 2 s 1 and v 2 = 16mm : s '. Since Re must be 
the same (for dynamic similarity), then 


1450 x 3 2 N 2 x 1 2 
14.5 ~ 16 


N 2 = 14400 rpm 


Since 


Q 

ND 3 


= const. 


N 2 D 2 


= const 


Q 2 ~ 1.839m 3 s 1 , &p 2 


To find ij 2 , use equation (2.10): 


1 ~ 7i 

1 - 7 2 


0.3 + 0.7 



4.38N m~ 2 


0.2 
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therefore ij 2 = 78.5% and 


1.839 x 1.169 4.38 ( . A p 

p _ x = m — 

2 0.785 1.169 V P 

= 10.26 W 


2.6 Pi = 100 MW output, therefore hydraulic power full size 

108.696 MW. The full size flow rate is given by 


108.696 x 10 6 = 20g x Q x 10\ Q = 554 m 3 s' 


Using 


gH 

nW = const ’ 


554 


0.9 


Q 

ND 2 

20 g 


= const 


5 g 


93.7 x 7 3 N 2 DV 93. 7 2 x 7 2 /VsD? 

Solving gives D 2 ~ 0.399 m and N 2 ~ 821.92 rpm. 
Using the Hutton formula for i] 2 

, 0.2 


| — = 0.3 + 0.71 , 

1 - i/2 \Re u ' 




Since R c = ojD 2 /v and r same for both machines, )] 2 = 87.57%. 
Therefore 

P 2 = 0.9 x 10 3 x 5g x 0.8757 = 38.66 kW 
If air is used as the flow medium in modelling 
' water = 1 miTTS -1 at 20 °C 
v air = Hninrs' 1 at 20 °C and 1 bar 
If Re is to be the same, then 


«->air 

17 


1 


w., jr = 13973 rpm 


Clearly this will pose stress problems, so one solution could be to follow 
other similarity laws, including characteristic number and to correct for 
efficiency variation. 


2.7 For dynamic similarity 


P«2 = J rhVjRTni) x HP 
Pm D 2 Po\ 


x R e x M n x . 


R = RJM = 8.3143/4.003 = 2.077 kJkg^'K 


i v -i 


For helium 
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From tables, /^ciium = 20.8 x 10 6 kgm 's 1 , p hclium = 0.1626 kg m 3 , 
t’ ajr - 1.568 x 10 -5 m 2 s _l and p an = 1.16 kg m -3 . 


P 02 


P 02 

Pm 

helium 

Pm 


mVRT 0 


m V/?r ()1 

D 2 P oi 

helium 

D 2 P oi 


12000V2.077 x 300 _ m A V0.287 x 300 
D 2 x 10 x 10 5 D 2 1 x 10 5 

Therefore m air = 3228.2 kg h' 1 


ND 


j ND 

V RT m 

helium 

VRT 01 


3000£> ND 

V2.077 x 300 ~~ \ 0.287 x 300 

from which N is 1115. 2rpm. But R E should also be constant since 

M 

0.1626 x 3000 x 2n x D : _ N x D 2 2n 
20.8 x 10~ 6 x 60 " 1.568 x 10“ 5 X 60 


Clearly, both criteria for N cannot be satisfied! One approach would be 
to test at 1115.2rpm and correct power (efficiency) by using R E cor- 
rections. Clearly there are problems when modelling using air. 


2.8 Again, as for Exercise 2.7 


P 02 _ r 

P<n ~ 1 


mV(RT 0l ) ^ ND ^ ^ 

~nIZ x x/JWt \ x Re x M n x . . 

D Pm \ (RT (n ) 


will be the same in both cases 


ND 


6000 D 


V(2.291 x 660) 


Therefore N = 25 662.5 rpm and 


\ (0.287 x 288)’ 


mV( 2.291 x 660) 

15 V(0.287 x 288) 

1.4 x 10 6 D 2 

gas “ ! 10 x 10 4 D 2 


m gas = -19.1 kgs 
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3.1 Vapour pressure head from Fig. 3.8 = 0.345 m: therefore, vapour pres- 
sure energy = 3. 385 J kg -1 . 

Since atmospheric pressure = 10 5 Nm -2 total energy at the pump is 

1 (p 

-5.5g + — r - 8.5 = 37.55 J kg 1 
6 10 3 5 

Now, NPSE a = 37.55 - 3.385 = 34. 165 J kg -1 (or NPSH A = 3.48m of 
water). If atmospheric pressure falls to 0.85 x lO^Nm^ 2 , then 

0.85 x 10 s 

NPSE a = -5.5g + ^ 8.5 - 3.385 

= 19. 16 J kg -1 (NPSH a = 1.95 m of water) 


3.2 Vapour pressure energy = 3. 385 J kg 
Case (a). Total energy at pump is 

10 s 

4g -l ^ — 50 = 89. 24 J kg“ 1 

* jq3 

NPSEa = 89.24 - 3.385 = 85.86Jkg-' 
(NPSHa = 8.75 m of water) 

Case (hi. 


NPSEa 


10 5 

4 s+ U ? XI - 3 


50 - 3.385 = ] 15.863 ks" 


(NPSH a = 11-81 m of water) 


NPSEa = 100 = 4g + - 100 - 3.385 

P = 1.64 x 10 5 N m -2 


3.4 If the tank is at vapour pressure NPSE A = 5.5 g - 50 = 3. 955 J kg 
and if free surface is at 3.5 bar gauge, then 

Kinrr c . 0.95 x 10 5 3.5 x 10 5 3.59 x 10 5 

NPSE x = 5.5? H 5 + T -50 r 

s 0.56 x 10 0.56 x 10 3 0.56 x 10 3 

= 157.51 J kg -1 (NPSH a = 16.06 m of butane) 


3.59 x 10 5 

NPSE - = ^ + 

loss = 29.24 J kg -1 


loss - 


3.59 x IQ 5 
0.56 x 10 3 


10 J kg" 


3.5 
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Figure A. 5 Solution 6.1. 


3.6 For the pump 


14507? // 80 \ 
30 V \3600/ 
(g x 50) 3/4 


0.217 


From Fig. 3.9, a = 0.07, therefore 

NPSH r = 0.07 x 50 = 3.5 m, NPSE R = 34.34 J kg” 1 
NPS&a = ~Xg + 72. 65 J kg” 1 and, in the limit. 

NPSEr = NPSE a - 10 = 34.34. NPSE A = 44.34 


Therefore 


44.34 = -Xg + 72.65 
X = 2.9 m suction lift 


6.1 C p = 5.193 kJ kg -1 K -1 and R = 8.3143 4.003 2.077 kJ kg” 1 K” 1 . 
Therefore 


C v = S.llOkJkg-'K” 1 , A: =1.67 


0.09 

"a = — x 
A 2 


24000?t 

30 


113.1ms 1 


Figure A. 5 shows the inlet velocity diagram and V ul = 3.32 ms *, V A = 
63.34 ms” 1 , V { = 63.43 ms” 1 and H 7 ! = 126.86 ms” 1 . Therefore 

p A = 4 x 10 5 /2.077 x 10 3 x 300 = 0.642 kg m” 3 


m = 0.642 x 63.34 x ^(0.09 2 - 0.06 2 ) = 0.144 kgs” 1 
a = V(kRT) = V(1.67 x 2.077 x 10 3 x 300) = 1020ms” 1 
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Therefore M nmax = 0.124 at A, outside diameter = 0.25 and 
24 000 t: 


u 2 = 


30 


x 0.125 = 314.16ms" 


Using Stanitz for slip (any of those quoted could be used) 
V U 2 , 0-6371 

= 1 — — , V u i = 262.34 ms 

u 2 12 

A p 


= 8.242 x 10 Jkg = u 2 V u2 (=h (p - h { „) 

P 

Therefore 

0.9 = 5.193 x 10 3 x T (n (T l)2 /T (n - l)/8.242 x 10 4 
T (n = 300 + 63.43 2 /2 x 5.193 x 10 3 = 300.39 K 
and T m !T m = 1 .048/1 . 


Since 


then 


~ = ( 1 . 048)' 1 '’ 7 1124 

' 01 


p... = 4 x 10 + 


63. 43 2 x 0.642 

9 


P (P = 4.013 x 10 5 Nm“ 2 


6.2 iv = 183.26 rads thus 

u 2 = 183.76 x 0.4 = 73.3 ms -1 

«ia = 32.07ms"" 1 . u iB — 18.33ms" 1 

V R2 = 0.8/te x 0.8 x 0.075 = 4.24 ms" 1 

Theoretical V u2 = 64.2ms" 1 . An estimate of ‘slip’ is needed for which a 
simple formula is 

77 ^ = 1 - 7 sin 25. V u , = 52.02 ms " 1 
64.2 7 

Assuming zero inlet whirl V ut = 0 

gH E = 52.02 x 73.3 = 4580Jkg" J 

considering the inlet annulus, flow through one inlet = 0.4m 3 s" 1 , 
therefore 

U a = 0.4/Q(0.35 2 - 0.2 2 ) = 6.17ms" 1 



222 Appendix 


V, 


(a) 


(b) 

Figure A.6 Solution 6.2. 

Thus at A, with VAj, = 0, the velocity triangles are (Fig. A. 6(a)) 

/i, A = 10.9°, /J 1B = 18.6° 

If 5 ms -1 swirl is imposed in the direction of rotation, the velocity 
triangles are now (Fig. A. 6(b)) 

/A a = 12.84°, /Ab = 24.84° 

(The incidence at A is tolerable, but that at B could lead to flow 
separation.) 

If point A is considered 

gH E = 4580 - 5 x 32.07 = 4424 J kg -1 
Thus energy rise reduction = 156 J kg -1 . 

6.3 (o = 30000 x nl 30 = 3141.6rads _1 , therefore 

u 1]p = 0.27 x 3141.6/2 = 424.1ms' 1 

— =1-0.63^. V u -> = 319.2ms' 1 
u 2 8 

C p AT = 319.2 x 424.1 - 1.354 x 10 5 Jkg“‘ 

AT = 134. 7K 

7"q2 = T m + 134.7= 424.7 K 






Using Stanitz 
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and 


Pn2 

Pin 


1 + 0.79 x 


134.7\ 3 ' 5 
290 1 


3.01 

1 


Using the simple assumption that p m = p m = 3.76 x 10"' Nm' 2 and. 
since eight blades 3 mm thick are used, the impeller peripheral flow area 
is reduced to 0.021 m 2 


U R2 = 4.5/0.021/; 


(A. 1 ) 


Since the only design data available is that Mn is not to exceed 0.95, it is 
necessary to assume a value of U R2 = 100 ms -1 . Also, since V 2 = 
334.5 ms -1 


334.5 

0.95 


a = V(1 14 x 287 x T 2 ), 


T 2 = 308.56 K 


Applying the gas laws, 

— = (^r 2 ) , p, = 1.23 x H) 5 N m -: 

PZ \ P 2 ' 

therefore 

p 2 = 1.23 x KF/287 x 308.56 = 1.39kgm -3 

This gives l 7 R; from equation (A.l) as 154ms -1 . Assuming = 
175 ms -1 . this vields V 2 = 364 ms -1 , T 2 = 365.4 K, p 2 = 2.22 x 
lO^Nm -2 and p 2 = 2.12kgm - \ which compares with p 2 from equation 
(A.l) of 1.224 kg m - -\ The iteration ceases, and 1/ R2 is 175ms -1 . 
Considering the diffuser 


,, 319.2 x 270 , N 

k u 3 = (free vortex law') 

V u? = 156.5 ms -1 


C p = 0.5 


Vi x 2.12 x 3 64 2 


and 

p, = 2.22 x 10 5 + 7.354 x 10 J 
Pi = 2.9554 x 10 ? Nm -2 

The Cp approach does not reflect the true changes since both p and p 
vary with temperature. 

6.4 At design, zero inlet whirl is assumed allowing for blade blockage, the 
flow area at outlet normal to flow = 7.455 x 10 -3 m 2 . Therefore 


V R2 = 100 x 10 3 /3600 x 7.4155 = 3.75 ms -1 
u-> = 37.96 m s -1 
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Figure A. 7 Solution 6.4. 


From the velocity triangle (Fig. A. 7) 

V u2 = 27.66 ms -1 

gH E = 27.66 x 37.96 = 1049.9 J kg' 1 

A slip allowance is needed. In this example Karassik’s formula is used, 
which gives 

F u 2^u 2 ideal = 0.65, V u2 = 18.04 

gH E (corrected) = 648. 9 J kg" 1 

Since the pump is of conventional design, r/ H will be used to get gH actu . A . 
Using gH = g //corrected the characteristic number k = 0.37. From Fig. 
2.6 // = 0.73 and, assuming a mechanical efficiency of 95% ^hydraulic 
= 76%. 

Recalculating k using 520 J kg" 1 yields k = 0.43. giving a slightly 
low'er r ] H of 75%, rj = 71.25% and gH = 514 J kg" 1 . Therefore 

power * S x 1C|5 ) x 57175 * 20 04 kw 

Turning to the diffuser and using Fig. 6.23, /?diff USC r = 0.75. Therefore, 
actual ^ - V\) 

Since V 2 = 7.03, V u2 = 5.95, U u3 = 3.72 and V 3 = 4.33 (using the free 
vortex principle). This gives 

— = ^t ; -( 7.03 2 - 4.33 2 ) = 1 1 .5 J kg~ 1 = extra generated pressure. 

P 2 

6.5 Q = 150/3600 = 0.0417 m 3 s _1 . 

The rotational speed is not specified, but since the pump is not of a 
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special type, a typical rotational speed would be 1450 rpni on a 50 kHz 
supply. Thus, 


0.75 


14507c 

30 


VO. 04 17 


(gH) 3 ' 4 


gH = 142.94 Jkg -1 
Ap = 1.429 x 10 5 N m -2 

Since the suction pressure is 0.85 x 10 2 Nm -2 the outlet pressure = 
2.2794 x 10 5 Nm“ 2 . 

As the discussion in Section 6.6 indicates there are a number of 
solutions available, and here for illustration it will be assumed that the 
outlet pressure applies to the back plate and shroud down to the wear 
ring. Suction pressure applies to the area inside the wear ring diameter. 
Area subjected to outlet pressure = 7i/4(0.14 2 — 0.065 2 ) = 0.0121 m 2 . 
Area subjected to suction pressure = 7i/4(0.14) 2 , therefore 

axial thrust = 2.2794 x 10 5 x 0.0121 — 0.85 x 10 5 x 0.0154 


= 1450N due to pressure 


Force due to atmospheric pressure of shaft end = 10 2 x 7i/4(0.065) 2 = 
331. 83 N. Total end thrust = 1781. 83 N. 


6.6 In the plenum 7<n = 1.25 x 10 5 Nm 2 . 7 0 1 = 300 K. At maximum 
radius in the inlet 


u n = 1575 x 0.15 = 236.25 ms 1 
For zero inlet whirl 

Wj = V\ + 236. 25 2 

VF, W, 

« = Vfrflro. = — = -ttff = °- 85 

3] vyRT\ 

Therefore 

W,= 0.85 V(1 .4 x 287 x 7,) 

and 

300 - T'i + 2 x 1.005 x 10 3 
Solution of equations (A.2)-(A.4) yields: 

V A = Vi = 165.3ms -1 , 


(A. 2) 


(A. 3) 


(A. 4) 


7, = 286.4 K 
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and from the gas laws p x = 1.063 x 10 5 N m 2 . This gives 
p l = 1.29 kg m -3 

rh = 1.29 x 165.3 x tt/ 4(0.3 2 - 0.1 2 ) = 13.43kgs -1 
Turning now to the outlet diameter of the impeller 
u 2 = 0.3 x 1575 = 472.5 ms -1 


Using Stanitz, 


u2 


394.54 


0.62^ \ 
12 / 


Therefore 


472.5 x 394.54 = 1.005 x 10 3 AT 


AT = 185.5 °C 


and 



(l + 0.82 x 


185.5\ 3 5 
300 / 


= 5.25 x 10 5 N m -2 


4.2 1 


7.1 Power output = 0.75 MW, therefore 

0.75 x 10 6 

hydraulic power = 0 955 x „ 92 x ‘ 0 , 97 
= 0.88 MW 


Since head = 2.5 m 


0.88 x 10 6 
2.5 g x 10 3 


36.2 m’ s 1 


Assuming hub:tip ratio = 0.4, then D H = 2m and 


V A 


36.2 

7t/4(5 2 - 2 2 ) 


2.195ms _1 


« lip = 7.15 x 5/2 = 17.875 ms 1 
u H = 7.15 x 2/2 = 7.15 ms -1 
gH E = 2.5g x 0.92 = 22.56 J kg -1 
Therefore, at the tip section, 

22.56 = 17.875 V ul (V u2 = 0), V ul = 1.262ms -1 
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Figure A. 8 Solution 7.1. 

From Fig. A. 8(a); (i 2 = 7°. //, = 7.53°, fi m = 7.265° (Ji m = (//, + /A)/2). 
a, = 60.1°. 

Using Clark ‘Y’ data (Fig. 4.8). C L (optimum) = 0.4" at a = 5.5°. 
which gives 

2.195 c 

22.56 = x 0.97 x - cosec 7.2-0 

2 x 17.870 s 

c/s = 0.149, s/c = 6.67 

At the hub section (Fig. A. 8(b)) V u] = 3.195. /i, = 2S.~9\ /A =17.07° 
and /?„, = 22.93°. Therefore 

2.195 c 

22.56 = — — x 0.97 x - cosec 22.9? 

2 x 7.15 5 

c/s = 1.15, s/c = 0.866 

Assuming h/c = 3 and c hub = 1.5/3 = 0.5 m then 

s = 0.866 x 0.5 = 0.43 m 


and 


Z = ' ^ - = 14.6 (i.e. 15 blades) 
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(a) 


53.15 


V u2 = 7.31 



Figure A. 9 Solution 7.2. 


to give .5 = 0.42 and c huh = 0.482m. 
At the tip section 


s = 


I X 5 


26 


c = 0. 163 m 


Stagger angles are shown in Figs A. 8(c) and (d) 


7.2 Since Q = 4m 3 s _l and V A = 15ms _1 , the flow area = 4/15 = 0.267m 2 . 
If the hub:tip ratio = 0.5 then D T = 0.673. Rounding: D r = 700mm, 
D h = 350 mm and V A = 13.86 ms -1 . 

Since A p = 175 N m~ 2 , gH = 175/12 and 

175 

gH E = — = 194. 4 J kg 1 

6 1.2 x 0.75 B 

Since w = 1450 X 7t/30, u H = 26.6ms _1 and u T = 53.15ms- 1 (zero 
inlet whirl is assumed). 
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At the tip section V u2 = 3.66 and (Fig. A. 9(a)) ji x = 14.62°, B 2 = 
15.65° and [S m = 15.14°. Thus 


194.4 
53. 15 2 


■ 13 - 86 . 2 xC L 
2 x 53. 15 2 L 


x - cosec 15. 14 

5 


Cy- =0.1378 

lip 


Similarly, for the hub section, V u2 = 7.31 and (Fig. A. 9(b)) /?, = 27.52°, 
p 2 = 35.7° and (J m = 31.61°. This gives 

C L - = 0.553 
s 
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Since the deflections are very small, flat plate data can be used at the 
tip, and at the hub a cambered plate. Assuming a camber of 4%, 

Q 

maximum — - is at a = 4° and C. = 0.8 
C D 

Therefore 


0.8 x - = 0.553, - = 1.447. 

s c 

Since h = 175 mm if blc = 3/1, c = 60 mm. This gives j hut , = 0.087 and 
Z = 12.64 = 13 blades. 

At the tip, C L = 0.4, sic = 2.903 and C,, p = 0.76m. For the tip 
section, Fig. A. 9(c) shows y = 71.35° while, for the hub section. Fig. 
A. 9(d) shows 7 = 62.48° with a camber radius of 162.1 mm. 

Comment : A constant section may be possible, with a 2% camber at 
the hub. 




w H = 7.15 x 2.5/2 = 8.9375ms ' 1 
V/ A = 50 j (5 2 - 2.5 : ) = 3.4ms " 1 


Hub spacing 


n x 2.5 

12 


0.654 


Therefore 


£ 2.5 

0.92 x 8.9375 2 


jc/s| hub = 0.765 

- — — — — — X 0.96 x 0.765 cosec// m 
2 x 8.9375 m 


(C L = 0.96 at C L /C D max @ 5° a) 
In Fig. A. 10(a) /? m = 24.75, and fi 2 = 20.83°. Now 
P m = = 24.75 


So that /?i = 28.67°, V u) = 2.72 ms 1 and 

gHz = 2.72 x 8.9375 = 24.3 J kg ' 1 
Hydraulic power = 50 x 10 3 x 24.3 = 1.22 MW 
Shaft power = 1.22 x 0.92 = 1.12 MW 

90 - y = 23.67° = setting angle in Fig. A. 10(b) 
y = 66.33° 
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Figure A. 10 Solution 7.3. 


7.4 Electrical output = 45 MW 

Hydraulic power = g x 11 x 500 x 10 3 = 53.955 MW 


q„ = 0.834 


45 

53.955. 


>h\ = 


0.834 

0.96 x 0.97 


0.896 


Therefore 

gltp = 11 g X 0.896 = 96.69 J kg -1 
= u\V ui , since V u2 = 0 


u,, n = 65.2 x x 4 = 27.31ms 1 
u hub = 9.56 m s _l 

V'uuip 3.54ms~ 1 and V A = 11.336ms _i 

From the tip velocity diagram (Fig. A. 11(a)). = 25.5° and ji 2 = 

22.54°. Assuming the free vortex ruie applies. 
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Figure A. 11 Solution 7.4. 


From the hub velocity diagram (Fig. A. 11(b)), = 92.8° and fi 2 = 

49.86°. For the draft tube 

inlet area = n 8 2 /4 = 50.27 m 2 
outlet area = 1.5 x 50.27 = 75.4 m 2 


Therefore the inlet velocity is 9.95 ms -1 and the outlet velocity is 
6.63ms' 1 . Theoretical draft tube regain = V2(9.95 2 - 6.63 2 ) = 
27. 52 J kg -1 ; actual regain = 0.8 x 27.52 = 22. 018 J kg -1 ; draft tube loss 
= 5. 502 J kg' 1 . 

Applying the energy equation to the draft tube 


Pa 

P 


9.95 2 

2 


10 5 

+ % " To 5 + 


6.63 2 

2 


+ 0 4- 5.502 


and, taking the draft-tube exit centreline as datum, 


^ = 58.361kg' 1 , p A = 0.584 x 10 5 Nm' 2 (abs) 
P 
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Figure A. 12 Solution 8.1. 


Comment: Since the vapour pressure is only 2500 Nm 2 there is no 
general risk of cavitation. 


S.l From steam tables iv = 3.992m’ kg 

tit = 20000/3600 = 5.56 kg s“‘ 

which gives 

20000 


2 = 


x 3.992 x 0.93 = 20.64 nFs 


3,-1 


3600 

Power output = 0.25 x l0 f1 Js~' and losses = 30%, therefore 


hydraulic input 

0.357 x 10 A 
5.56 


0.25 x 10° 
0.7 


= 0.357 x 10° Js' 


= 0.642 x 10 5 = uAV„ 


V A = 0.7// = 0.7 x x ^ = 54.98D n 


and it = 78.54 D m . 

From the velocity triangles at mean blade height (Fig. A. 12) 
F ul = V A cot 20 = 151.06D m 

F u2 = V A cot 20 - u = 151.06 D m - 78.54 D m = 72.54£> m 


u 


Figure A. 13 Solution 8.2. 




The velocity diagrams can now be drawn for the mean section, and 
Figure A13 results 

.'. V A = 15.1ms -1 
AV U = 396.8 m s -1 
From tables, v g = 3.42 m 3 kg -1 
at the outlet 
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115.1 x 71 x 1.5 x 0.035 


= 5.55 kgs 


Diagram work = 235.6 x 396.8 

= 93.48 kJ kg -1 per stage 

Since 4 stages 
Total diagram work 


Total diagram work = 4 x 93.48 

= 373.94 kJ kg’ 1 

from tables inlet enthalpy = 3793 kJ kg -1 
outlet enthalpy = 2683 kJ kg -1 
.'. enthalpy change = 1110 kJ kg 1 
373.94 

'/thermal — “TTSTT — 33.7% 


8.3 Total diagram work = (muAV n ) x 10 

= 16.5 MW 


A1 , 16.5 x 10 h = 23571 J kg 1 for one stage. 

uAV.. = 

u 10 x 70 


300071 1.5 

" m = ~ x — = 235.62 ms 


Figure A14(a) can be drawn. 

.-. AV U = 100.04ms -1 
V 7 ,,: = 67.81 in s -1 
l 7 ul = 167.81 m s - 1 
fix = 55.86° 

(U = 30.79° 

7 : = 25.09° 

V 7 ! = W 2 = 195.34 ms -1 
V 2 = W t = 120.83 ms -1 


195. 34 2 

7111 = 7 + F x U45 = 100 °- 

7, = 983.3 K 

a = V( 1.33 x 287 x 983.3) = 612.66ms" 


maximum AT. = 0.32 


Pi) ~ P 
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Figure A. 14 Solution 8.3. 


Since M n is low, p can be assumed incompressible and k = 1.33. Thus 

P °. ~~ 2 P = 1.0256 p = P „„„ „ = 3.543 x 10“ 6 p 

\pv 2 y 287 x 983.3 ^ 


7.5 x 10 5 - p = 1-0256 x 195.34 2 x 3.543 x 10“ 6 p 

p = 7.0138 x lO-'Nm 2 
p = 2.485 kg m“ 3 


e = 2 ^ = 28 - 17m3s-1 


Therefore 
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This gives 

28 17 

blade height = — — = 0.06 m 

6 100 x n x 1.5 

Annulus dimensions are 1.56m and 1.44 m. Since free vortex blading 
applies. V u x R = const. n, jp = 245.04ms ' 1 , 17 U | (ip = 161.35ms -1 , 
V u 2 tip = 65.2ms - ’, n hub = 226.2ms -) , F ul hub = 174.8ms -1 and 
V u2 hub = 70.64 ms -1 . 

For the tip section (Fig. A. 14(b)) /?, = 50.08°, /T = 29.08° and n = 
21 °. 

For the hub section (Fig. A. 14(c)) /?, = 62.8°, /? 2 = 32.7° and c = 
31.1°. 

A constant section based on mean section will give acceptable in- 
cidence angles. 


8.4 For air. C p = 1.005 x 1 0 3 J kg 1 so that for the first stage 
1.005 x 10 3 x 20 = 0.925«AF U 
uAV u = 21.73 kJ kg -1 
Since F a = V t = 140ms -1 


288 = T, + 140 2 /2 x 1.005 x 10 3 , T { = 278.25 K 


_ 1 1 iv. | 1 11 

/>, \ i 


/;, = 0.895 x KPN M -: 


P\ 


0.895 x 10 s 


= 1.121 kgm' 


287 x 278.25 
a = V(1.4 x 287 x 278.25) = 334.4 ms' 


to give 

VF, = 0.95 x 334.4 = 317.65ms -1 
Since V A = then V u] = 0 and 

u = V(317.65 2 - 140 2 ) = 285.13ms -1 
AF = F u2 = 76.21 ms -1 

From Fig. A. 15(a) (j lT = 26.15° and /? 2 y = 33.83°. Since u v = 
285.13 ms -1 


285.13 


600071 Dr 
X — - 

30 2 ’ 


D t = 0.908 m 


D H = 0.6 x 0.908 = 0.545 m 

p — 1.121 kgm -3 
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Figure A. 15 Solution 8.4. 


and therefore 


m = 1.121 x 140 x -(0.908 2 - 0.545 2 ) = 65 kgs -1 


1 + 0.89 x 


= 1.233:1 


If free vortex conditions, F u i H and V u2H are found: 

V ul = o, 

V u2 h x 0.545 = 76.21 x 0.908 = 126.97 ms -1 
u H = 171.14ms -1 

From Fig. A.15(b) ft lH = 39.285° and [i 1H = 72.49° 

8.5 Since the fluid is helium C p = 5.193 x 10 3 , M = 4.003 and 
R = ^ = 2.077 x 10 3 J kg -1 K -1 


Thus C v = 3.116 and 


n = yr = 167 

C-v 


n - 1 


= 2.493 


/ A7V" 1 *' 

3.75 - (l + 0.85 x -) 

Therefore AT = 287.8 °C and A T slagc = 23.23 °C. 

From data sheets = 0.87 (Fig. 1.17) and fl = 0.855 (Fig. 4.20). 


stage 1 — 1 = 11 + 0.87 x 
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For the first stage P, )2 = 5.71 bar and T (a = 


= (l + 0.87 x 

MG' stage 2 V 


373.23 K. 

23.23 \ 2 43 
373.23/ 


P 03 = 6.559 bar, T m = 396.46 K 


and so on up to stage 5: 

P ()5 = 8.348 bar, 7 05 = 442.92 K 

These are inlet conditions (stagnation) to the 5th stage. 
At the 5th stage tip diameter 



300071 

30 


549.78 ms" 1 


Assuming 50% reaction at the tip, 

5.193 x 10 3 x 23.23 = 0.855 x AF U x 549.78 
AF U = 256.63 ms -1 

also (Fig. A. 16(a)) l/, = 189.45 ms -1 , e = 22.74°, /?, = o 2 = 73.43° (to 
axial) and fi 2 = a \ = 50.69° (to axial). Applying the free vortex laws 

F ul x R = const 

F u2 x R - const 


therefore the annulus dimensions are needed. 

At inlet to stage 5. T {)5 = 442.92 K and P (l3 = 8.348 bar which gives 

189. 45 2 


442.92 = Fs + 


2 x 5.193 x 10 3 


Therefore 


T 5 = 439.46 K 

Pos __ / 442.92 \ 2 49: 
p s ~ V439.46/ 


Ps = 8. 157 bar 


P 


8.187 x 10 5 
2.077 x 10 3 x 439.46 


0.897 kg m -3 


Q 


350 

0.897 


390.23 m 3 s -1 


390.23 

120 



d h = 2.524 m 


Thus F u , hub — 203.26ms ', V u2 hub = 559.13ms 1 and 
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^hub 


300071 

30 


2.524 

2 


396.47 ms 1 


From Fig. A. 16(b) ci\ = 59.44°, /?, = 58.16°. a 7 = 77.89° and = 
143.58°. 

Returning to the tip section, Fig. A. 16(c) shows the camber section, 
where 0/2 = 8.63°, 6 = 17.26°, y = 48.69 - 8.63 - 40.06° to the axial 
direction, = 48.69° to axial and P 2 = 7 ~ 8.63 = 31.43° to axial. 


« 2 = 


3.3 


x 18 = 29.7 ms" 


9.1 
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2.15 

Uj = — — x 18 (mean dia) = 19.35 ms 1 

The velocity triangles (Fig. A. 17) may be drawn at the mean diameter, 
to give 

V u2 + W 2 cos 85 = 29.7, 29.7 - IV 2 cos85 = V u2 

V / u2 tan20 = V A2 = VF 2 sin85 

therefore 

(29.7 — W 2 cos 85) tan 20 = VF 2 sin85 
F 2 = 30.64 ms- 1 . 

V u2 = 28.79 ms-', V A2 = 10.48 ms-' 

A 2 — 3 m : . Q = 3 X 10.48 = 3 1.44 mV 
A 3 = 3.1m 2 , V A3 = 10.14 ms - ', F 4 = 10.73 ms - ' 

Thus 

V u3 = W 3 cos 25 - 19.35 

W 3 = 10.14 cosec 25 = 23.998ms-' 

V u3 = 4.65. V 3 = 11.16ms-‘ 

Therefore 

Euler gH = 28.79 x 29.7 + 19.35 x 4.65 
= 855.06 + 89.98 = 945.04 Jkg -1 

and 

runner power = 945.04 x 31.44 x 10 3 = 29.71 MW 
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The total ‘head’ drop is found by applying the energy equations. Using 
(1) to indicate inlet to guide vanes, (2) as inlet to runner, (3) as inlet to 
draft tube and (4) as draft tube outlet, then from (1) to (2) over guide 
vanes 

Vj Vj l/s 

g^i + y + gZ 1 = gh 2 + y + gZ 2 + 0.05 y (A. 5) 

From (2) to (3) over runner 

v? v\ wl 

gh 2 + y + gZ 2 = ghy 4- y 4- gZ 2 + 0.2 y 4- work done (A. 6) 

From (3) to (4) over draft tube 

V\ V 2 4 V\ 

ghy + y + gZy = gh A + y + gZ 4 + 0.5 y (A. 7) 

From equation (A. 7) 

gh 3 = g6 A - 31.14 4- 115.13 - 31.88 
= gh a + 52.11 

From equation (A. 6) 

g'U + + g X 4 = (g/, A X 52.11) 4- 

7 3 998 : 

+ 945.04 + g x 25 4- 0.2 x ~ 

w 1 


gh 2 = gh A ■+ 640.26 


From equation (A. 5) 


ghj = gh A + 640.26 4- + 0.05 x 

= gh A 4- 1133.135 

Total energy drop = 1133. 135 J kg -1 , therefore // hydraul , c = 0.834. 
9.2 From Fig. A. 18(a) 

C P 


1 TS 


0.86 - 


"01 ~ 

hoi hy ss 

850 - 670 

a r« 


To, - r t 


C n VTa, - 7\. 


AT is = 209.3 K 


AT js for nozzle = 104.65 since 50% reaction. 
„ = Cp ( r<)1 Z 7 2 ) = n 07 

/N C p (T m - r 2s ) 09 ’ 


Tq 1 Tj — 0.97 x 104.65 
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(b) (c) 

Figure A. 18 Solution 9.2. 


T 2 = 748.5 K 

Thus nozzle velocity = V(1.145 x 10' x 101.51 x 2) = 482.14ms -1 . 
The inlet triangle becomes Fig. A. 18(b). Since 

u T = 50 000^/30 X = 431.97ms -1 


a, = 26.37° 


W 2 = V R2 = 214.15 ms -1 

Peripheral nozzle area = n x 0.175 x 0.03 x 0.85 = 0.014m 2 . Now, p 2 
is required in order to find p 2 . 



P2 


850 

,748.5. 


= 1.67 
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P2 


Thus 


P 2 


3 x 10 s 
1.67 

1.8 x 10 5 


= 1.8 x 10 5 N m~ 2 


= 0.838 kg m 


-3 


Assuming V A3 

670 = 7 3 + 
foi / 

P 3 V 


287 x 748.5 
m = 0.014 x 214.15 x 0.838 = 2.52kgs~' 
V 3 = 300 

300 2 


2 x 1.145 x 10 3 ’ 
850 N4 


r 3 = 630.7 K 


P 3 
2.52 


640.7/ 

= 300 x A 3 


p 3 = 0.968 x 10 N m 
/; 3 = 0.502 kg m -3 


0.502 

A 3 = 0.0167 m 2 


D >‘ V; 


x 0.0167 = 0.146 m 


Thus 


50000tt 0.146 „„„ „ , 

u 3 tip = — — — x — - — = 382.2 ms 


30 2 

From Fig. A. 18(c) 

W 3 = 485.9 ms' 1 , a 3 = V(1.33 x 287 x 630.7) = 490.66ms' 1 
Thus M n = 0.99 

9.3 Since zero outlet whirl w out is not needed and 

u T = x 50 000 x ^ = 392.7ms' 1 

The inlet gas triangle (ideal) is shown in Fig. A. 19. 

479.4 = V(l. 145 x 10 3 x AT x 2) 

AT = 10 °C = T 0I - T 2 
T q , - T 2 


7n 


To, - T 2s 


= 0.97, 


7% = 839.69 K 


Since 50% reaction, isentropic turbine drop = 20.62° and 

C P (850 - T 3 ) 


*7 ts — 0.86 — 


C p (20.62) 


r 3 = 832.27 K 
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^Rl 

<=W,) 


Figure A. 19 Solution 9.3. 

and T 3s = 850 - 20.62 = 829.38 K 


Pm 

P.1 


850 


= 1.103 


V829.38 
p 3 = 2.72 x 10 5 N m” 2 
For the nozzle guide vanes 

Pm = /_850_' 4 
p 2 V839.69 


= 1.05 


and 


p 2 = 2.857 x 10 5 Nm” 2 


2.857 x 10 , 

p 2 = ~ = 1 • 186 kg m ~ 


287 x 839.69 

Now, m = 0.56 kgs -1 , therefore Q = 0.472 m's” 1 and, for the throat 

0.472 


A = 


479 


= 9.86 x 10” 4 m 2 


The ‘diagram' power is 

m x uV ul = 0.56 x 392.7 2 = 83.36 kW 

9.4 At plane 2, A - = 50000, u 2 = 392.7ms” 1 and u 3T = 274.9ms” 1 . From 
Fig. A. 20(b) 

V 2 = 479.4ms” 1 = V(2 x 1.145 x 10 3 x AT aCT ) 

A Pact = 100.36 (nozzle) 

Turbine AT = 2 x 100.36 = 200.72 °C, T (n = 1000 K. 

200.72 


, /TS = o,83 = 1Q00 _ ^ 

Pus _ 3.5 x 10 3 / 1000 

P? 


= 758. 18 K 


p 3 \758.19 
1.157 x 10 s 


p 3 = 1.157 x 10 5 NnT 


P3 = 


= 0.532 kg m 
287 x 758.19 6 

Power generated = 120 kW, therefore 


-3 
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V A 

(=^3) 



(c) 



Figure A. 20 Solution 9.4. 


— — — = 1.145 X 10 3 x 200.8 x m 
U.o 

m = 0.653 kgs' 1 

Area of flow = -(1.05 2 - 0.02625 3 ) = 0.0081 m 2 . therefore 
4 

0.653 = 0.0081 x 0.532 x V' A3 , V A3 = 151. 2 ms~‘ 

Thus at the exducer the maximum diameter velocity triangle is given by 
Fig. A. 20(c): 

fi 3 = 28.8° 

Wi = 313.7 m s -1 

a = V(1.33 x 287 x 758.19) = 537.96ms- 1 
relative M n = 0.583 

At the exducer the hub diameter n H = 68.73 and (Fig. A. 20(d)) /? 3 = 
65.56°. 

Second part of solution 

If 0% reaction then ST nozi[c = turbine AT and, following the solution 
through 

p 3 = 0.829 

m = 2 x 0.653 = 1.306 kgs -1 
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V A3 = 194.49 ms' 1 

a is now 579ms' 1 and W 3 = 336.74ms' 1 , which gives 

relative M n = 0.582 

9.5 In the pumping mode Q = 74 m 3 s' 1 in the draft tube, which gives V 2 = 
2. 16 m s _l . V 1 = 4.87 m s _1 and V 7 ,,, = 3.52 ms -1 . The friction loss is 


4 x 0.005 x 78 x 3.52 2 


1.76 J kg -1 


2 x 5.5 

Applying the energy equation from draft tube exit to its inlet (planes 2 
to 1) 

10 5 2.16- p, 4.87 2 n , ^ 

_ + _ + 29f ,,£L + -_ + 0+ 1.76 

thus pi = 0.373 x 10 6 Nm“ 2 . 

In turbining mode Q = 76.5 m 3 s' 1 , V x = 5.03ms' 1 and V 2 = 
2.23ms' 1 . Draft tube loss is 

f5.03 2 2.23 2 \ 


0.8 


= 8.13 J kg' 


Thus the lowest draft tube pressure is in pumping mode and, since the 
inlet head is 318 m 

Pin,,., = 318g x 10 3 = 3.1196 x 10 6 N m' 2 
Applying the energy equation from draft tube inlet to its exit, 
p i , 5.03 2 
10 3 2 


10" 2.23“ 

0 = ^ + -y- + 29g + 8.13 


Pi = 0.383 x 10 6 Nm - 
0.373 x 10 6 


= 0.1197 


3.1196 x 10 6 

Applying modelling concepts (Chapter 2), if the full size machine has 
subscript "p‘. and model has subscript ‘m\ then 


Q P = 74 . 


H p = 318, 
Qp Qn 


N p = 333, 

Up 


% 


88.5% 


V / >' 


N m D™ 


N 2 D 2 N 2 D 2 

1 ' p iV p 1 ’ m iy m 


74 


333 x 10 : ' 
318 


3000 xl 3 ’ 

H, 


Q m = 0.667 m s 


3c-' 


333’ x 10 2 3000 2 x l 2 

Using the Hutton formula (Table 2.1) 


5 = 258 m 


1 - > 1 p 
1 - >lm 


= 0.3 4- 0.7 


/ /?em V’- _ 
\Re p/ 


0.3 + 0.7 


( 3000 _1_ 

333 X 10 2 


0.2 
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Since //,, = 88.5, rj m = 84.3% and therefore 

74 x 10 3 x 258# 

model power = 0 ~g 43 ~ 2.22 MW 


9.6 Total energy drop over unit is 


10 3 \ 10 5 

100g + Jip) ~ Jo 5 _ °' 5g = 99 - 5?Jkg ' 


Hydraulic power is 


50 x 10 6 


= 10 3 x Q x 99.5 g 


Q = 56.92 gm 3 s -1 

The radial velocity at entry to runner is 


7r x 6 x 0.3 


= 10.07 m s" 


The velocity at entry to draft tube is 


- x 4 2 

4 


The velocity at exit from draft tube is 


= 4.53ms' 


- 2.011 ms' 


Looking at the draft tube (DT) 

^ + i^ + 0 ^ + ^ + 2g + DTIOSS 

P 2 p 2 


DT loss = 0.83 




P in = 0.1175 x 10 6 Nm 2 (absolute) 
Considering now the guide vane path, from Fig. A. 21(a): 
u = 29.3 x § = 87.9ms” 1 
10.07 = V, sin 25, V t = 23.83 ms” 1 
therefore guide vane loss is 

73 S3 2 

0.06 x = 17.03 J kg” 1 
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Figure A. 21 Solution 9.5. 

Applying the energy equation through the guide vanes 

100g = h? + + 0 + 1703 

giving the pressure at inlet to the runner of 0.68 x 10 6 Nm 2 . The inlet 
velocity triangle is (Fig. A. 21(b)) a 2 = 21.69° and p 2 = 31.02°. 

Construction as shown by Fig. 7.7, gives Fig. A. 21(c); the nozzle 
angle to triangle is 


5 777 

a' = tan -1 ' = 20.49 

14.11 
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